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A model  solar-thermal  concentrator-collector  that  is  l/6th  scale  (in  cross 
section)  relative  to  a system  considered  to  be  of  practical  dimensions  for 
power  generation  has  been  constructed  and  tested.  It  consists  of  an  Archimedes 
(i.e. , Fresnel)  mirror-concentrator,  5 cm  ID  x 1.5  m long  glass  receiver  pipe, 
0.3  cm  thick  graphite  absorber,  and  gaseous  (Ar)  heat  transfer  medium. 
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moves  east-west  within  two  1.67  m dia.  tracks  located  near  its  ends.  The 
glass  receiver  pipe  holds  the  thin  material  that  absorbs  the  concentrated 
solar  radiation  focused  on  it  and  transfers  it  as  thermal  energy  to  the 
flowing  gas,  which  exits  at  high  temperature.  This  receiver  pipe  remains 
fixed.  The  insulation  (2.5  cm  thick)  surrounding  the  pipe,  except  for  a 
narrow  (2  cm  wide)  longitudinal  aperture,  rotates  in  such  a manner  that 
alignment  between  the  aperture  and  the  tracking  concentrator  is  always 
maintained.  The  collector  is  tilted  at  41  (the  latitude  of  State  College, 
PA)  from  the  horizontal  so  that  its  axis  is  parallel  to  the  earth's 
rotational  axis  when  the  device  is  oriented  north-south.  The  small  size 
of  the  aperture  and  the  internalized  flux  absorption  and  heat  transfer 
within  the  cavity  receiver  reduce  thermal  and  radiative  losses  substantially. 
An  alternative  0.09  cm  thick  copper  absorber  with  a selective  oxide  surface 
^has  also  been  employed.  /sj  f-- 


On  clear  days  (direct  solar  flux  ca.  800  stagnation  (l.e.,  no-flow 

condition)  temperatures  of  ca.  370°C  have  been  reached  at  the  center  of  the 
experimental  device,  where  end  effects  are  minimal.  Experiments  have  been 
performed  with  gas  flow  rates  up  to  ca.  100  llters/min  (200  RE  2000), 
resulting  in  a bulk  fluid  temperature  rise  in  the  axial  (flow)  direction  of 
up  to  285^.  These  experiments  have  been  utilized  in  a study  of  heat 
transfer  processes  within  the  collector.  A computer  model  of  this  system 
which  takes  into  account  most  of  the  influential  variables  (e.g.,  combined 
forced-free  convection,  establishment  of  velocity  and  temperature  profiles, 
and  axial  heat  conduction)  has  been  developed.  Its  behavior  agrees  well 
with  observations.  Also,  after  ca.  200  h of  high-temperature  operation, 
materials  degradation  was  measured  and  found  to  be  very  low.  ^ 

The  good  agreement  between  the  analytical  and  experimental  results  permitted 
the  formulation  of  a second  computer  model,  of  the  full-scale  device,  with 
confidence.  The  effect  of  various  parameters  such  as  Insulation  thickness, 
aperture  width,  and  antlreflectlve  coatings  has  been  investigated  in  an 
optimazation  sutdy  of  the  proposed  full-scale  collector.  It  was  concluded 
that  the  aperture  width,  aperture  emittance,  and  mirror  reflectance  exert 
the  strongest  influence  on  the  exit  temperature  and  collector  efficiency. 
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ABSTRACT 

A model  solar-thermal  concentrator-collector  that  is  l/6th 
scale  (in  cross  section)  relative  to  a system  considered  to  be  of 
practical  dimensions  for  power  generation  has  been  constructed 
and  tested.  It  consists  of  an  Archimedes  (l.e.,  Fresnel)  mirror- 
concentrator,  5 cm  I.D.  X 1.5  m long  glass  receiver  pipe,  0.3  cm 
thick  graphite  absorber,  and  gaseous  (Ar)  heat  transfer  medium. 

The  mirror-concentrator  utilizes  40  suitably  tilted  reflective 
strips,  each  1.6  cm  wide  by  1.5  m long;  the  strips  are  mounted  In  a 
parallel  planar  array  which  moves  east-west  within  two  1.67  m dia. 
tracks  located  near  its  ends.  The  glass  receiver  pipe  holds  the 
thin  material  that  absorbs  the  concentrated  solar  radiation  focused 
on  It  and  transfers  It  as  thermal  energy  to  the  flowing  gas,  which 
exits  at  high  temperature.  This  receiver  pipe  remains  fixed.  The 
insulation  (2.5  cm  thick)  surrounding  the  pipe,  except  for  a narrow 
(2  cm  wide)  longitudinal  aperture,  rotates  in  such  a manner  that 
alignment  between  the  aperture  and  the  tracking  concentrator  Is 
always  maintained.  The  collector  Is  tilted  at  41°  (the  latitude  of 
State  College,  PA)  from  the  horizontal  so  that  its  axis  is  parallel 
to  the  earth's  rotational  axis  when  the  device  is  oriented  north-south. 
The  small  size  of  the  aperture  and  the  internalized  flux  absorption 
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and  heat  transfer  within  the  cavity  receiver  reduce  thermal  and 
radiative  losses  substantially.  An  alternative  0.09  cm  thick 
copper  absorber  with  a selective  oxide  surface  has  also  been 
employ>.a. 

_2 

On  clear  days  (direct  solar  flux  ca.  800  W m ),  stagnation 
(i.e.,  no-flow  condition)  temperatures  of  ca.  370®C  have  been 
reached  at  the  center  of  the  experimental  device,  where  end  effects 
are  minimal.  Experiments  have  been  performed  with  gas  flow  rates 
up  to  ca.  100  liters  min  ^ (200  < Re  < 2000),  resulting  in  a bulk 
fluid  temperature  rise  in  the  axial  (flow)  direction  of  up  to 
285'’C.  These  experiments  have  been  utilized  in  a study  of  heat 
transfer  processes  within  the  collector.  A computer  model  of  this 
system  which  takes  into  account  most  of  the  influential  variables 
(e.g.,  combined  forced-free  convection,  establishment  of  velocity 
and  temperature  profiles,  and  axial  heat  conduction)  has  been 
developed.  Its  behavior  agrees  well  with  observations.  Also,  after 
ca.  200  h of  high- temperature  operation,  materials  degradation  was 
measured  and  found  to  be  very  low. 

The  good  agreement  between  the  analytical  and  experimental 
results  permitted  the  formulation  of  a second  computer  model,  of  the 
full-scale  device,  with  confidence.  The  effect  of  various  parameters 
such  as  insulation  thickness,  aperture  width,  and  antireflectlve 
coatings  has  been  investigated  in  an  optimization  study  of  the 
proposed  full-scale  collector.  It  was  concluded  that  the  aperture 
width,  aperture  emlttance,  and  mirror  reflectance  exert  the  strongest 
influence  on  the  exit  temperature  and  collector  efficiency. 
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CHAPTER  I 


SOLAR  THERMAL  PROCESSES 


A.  Introduction 


Solar  thermal  power  generation  is  one  method  currently 
being  investigated  as  a potential  means  of  converting  solar  energy 
to  electricity.  This  involves  a number  of  steps:  collection  of 
the  solar  flux  as  thermal  energy;  transfer  of  this  energy  to  a 
suitable  fluid;  conversion  to  mechanical  energy  via  a heat  engine; 
and,  finally,  coupling  of  the  engine  to  an  electric  generator. 

Except  for  differences  in  energy  sources,  the  process  is  essentially 
the  same  as  that  employed  in  large  nuclear  and  fossil-fuel  electric 
generating  stations.  Thus,  research  centers  on  the  solar-to-heated 
working  fluid  steps,  these  having  a decisive  influence  on  overall 
system  efficiency. 

The  heat  engine  (usually  a Rankine  or  Brayton  cycle)  has  the 
Carnot  efficiency  as  a theoretical  upper  bound.  Because  for  a 
specified  heat  rejection  temperature  the  Carnot  efficiency  is  a 
function  solely  of  the  operating  fluid  temperature  at  the  inlet  of 
the  heat  engine,  a high  exit  fluid  temperature  from  the  solar  col- 
lector is  a prerequisite  for  a realistic  design.  The  required  high 
temperature  can  only  be  achieved  by  a focusing  (i.e.,  concentrating) 
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collector.  Focusing  collectors  utilize  optical  systems  - reflectors 
or  refractors  - to  increase  the  (normally  low)  intensity  of  solar 
radiation  impinging  on  the  absorber.  Higher  energy  flux  on  that 
surface  results  in  a smaller  surface  area  for  a given  total  amount 
of  energy  transfer,  and  correspondingly  reduced  thermal  losses,  thus 
permitting  higher  temperatures  to  be  reached. 

Unfortunately,  focusing  collectors  have  a number  of  disad- 
vantages when  compared  with  the  simpler  non-focusing  flat-plate 
collectors: 

(1)  optical  transmission  losses  become  significant  because 
of  the  imperfect  optics  of  a real  concentrator; 

(2)  generally,  only  the  direct  component  of  solar  radiation 
can  be  utilized; 

(3)  the  added  cost  and  complexity  of  concentrator  and 
tracking  systems  must  be  borne. 

It  is  obvious  that,  for  a given  incident  flux,  the  maximum 
useful  output  of  a collector  will  be  obtained  when  the  sum  of  the 
thermal  and  optical  losses  is  a minimum.  There  is  a trade-off  to 
be  made  here,  for  at  high  concentration  ratios  (i.e.,  small  receiver) 
much  of  the  energy  will  miss  the  receiver  altogether,  while  at  low 
concentrations  the  large  receiver  surface  will  lose  much  of  the 
energy.  A large  number  of  designs  have  been  proposed  that  attempt 
to  improve  collector  performance.  These  designs  can  usually  be 
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placed  in  one  of  two  broad  categories;  long  linear  collectors  (also 
termed  "distributed  collection  systems")  of  the  solar-farm  type,  or 
tower  receivers  with  flat  heliostats. 


B.  The  Solar  Farm 


Long  linear  collector  designs  covering  large  land  areas  have 
been  extensively  treated  by  the  Meinels  (1),  and  by  Duffle  and  Beck- 
man (2).  The  absorber  pipe  is  often  placed  within  an  evacuated  glass 
envelope,  and  the  solar  flux  is  most  frequently  concentrated  by  either 
a parabolic  trough  reflector  (3)  or  the  less  common  Fresnel  refractor 
(4) , which  has  been  investigated  by  Nelson  ^ al.  (5) . A selective 
surface  is  proposed  for  the  absorber  exterior,  and  part  of  the  enve- 
lope is  silvered  to  be  highly  reflective  (see  Figure  1)  (4,6).  Heat 
at  high  temperature  is  extracted  via  a fluid  flowing  through  the 
absorber  pipe,  exiting  at  ca.  500°C;  the  fluid  can  be  sodium  (4) 
or  steam  (7).  Parasitic  fluid  pumping  losses,  for  a receiver  up 
to  1 km  in  length,  have  been  shown  to  be  negligibly  small  in  compar- 
ison with  the  power  produced  (7) . 

The  theoretical  analysis  of  such  a collector's  performance 
has  frequently  been  extremely  simplified  (1,3,4).  A more  complex 
mathematical  formulation  of  heat  exchange  within  the  receiver  and 
to  the  ambient  environment,  as  well  as  the  effects  of  several  fluid 
types,  the  inlet  temperature,  and  the  flow  rate,  is  given  by  Lee  and 
Schimmel  (6) . For  turbulent  flow,  the  percentage  of  energy  extracted 
by  the  fluid  from  the  incoming  solar  flux  is  estimated  at  40-50%; 
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Figure  1.  Cylindrical  parabolic  collector. 
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selective  coatings  on  the  absorber  were  found  to  be  relatively  unim- 
portant if  the  surrounding  envelope  is  silvered  (6) . A fairly 
comprehensive  analysis  is  presented  by  Edenburn  (8) , who  consid- 
ered both  an  evacuated  and  a non-evacuated  annular  receiver,  with 
Therminol  66  as  the  working  fluid.  Experimentally  determined  and 
calculated  efficiencies  of  heat  transfer  (not  defined)  are  stated 
to  be  in  excellent  agreement,  and  range  between  35-62%,  with  fluid 
inlet  temperatures  spanning  475-578®K  (8) . 

Recent  experiments  with  water  (9)  and  an  evacuated  envelope 

about  the  receiver  pipe  resulted  in  peak  efficiencies  at  SOC’C 

(defined  as  net  heat  delivered/perpendicular  solar  flux)  of  50%  and 

30%,  respectively,  for  a selectively  coated  and  a black  painted 

tube.  Cobble  and  Smith  (10)  analyzed  in  detail  a projected  5 kW 

solar  electric  generating  station  with  air  as  the  working  fluid  in 

an  open  Brayton  cycle  and  predict  an  overall  maximum  efficiency 

(cycle  power/flux  on  concentrator)  of  about  6%,  with  an  optimum 

z (z  = , over  the  cycle)  of  approximately  2.9.  No 

envelope  or  insulation  was  considered  around  the  absorber  pipe. 

This  conversion  efficiency  is  only  about  one-half  the  projected  12% 

calculated  by  Meinel  and  Meinel  for  a steam  turbine  cycle  (7) . The 

2 

latter  authors  propose  a 1.61  km  "farm,"  with  a collector  fill 
factor  of  50%,  operating  over  a AT  = lOO'C,  and  producing  a peak 
electrical  output  of  90  MWe  at  the  quoted  efficiency. 

Optimization  studies  have  been  done  by  many.  For  example, 
Singh  and  Cheema  (11)  studied  the  performance  of  a cylindrical- 
parabolic  collector,  while  Howell  and  Bannerot  (12)  investigated 
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maximizing  cycle  work  output,  and  Turner  (13)  considered  minimizing 
the  piping  and  insulation  requirements.  The  entire  solar  farm 
concept  has  recently  been  the  subject  of  very  comprehensive  research 
efforts  by  the  University  of  Minnesota  and  Honeywell  (14) ; however, 
further  evolutionary  development  of  the  concept  appears  to  have  been 
brought  nearly  to  a halt  in  favor  of  the  solar  tower  approach. 

C.  The  Solar  Tower 

The  solar  tower  concept  is  currently  under  intensive  develop- 
ment by  DOE  (the  federal  Department  of  Energy)  with  progressively 
larger  test  facilities  planned.  This  concept  differs  from  others 
in  that  it  replaces  most  of  the  plumbing  required  for  transport  of 
heat  to  the  power  plant  by  optical  transmission  of  the  solar  flux 
to  a central  receiver.  Very  high  concentration  ratios,  on  the  order 
of  1000,  are  achievable;  thus,  the  solar  flux  can  theoretically  be 
utilized  at  the  highest  temperatures  tolerated  by  the  materials 
employed,  with  relatively  low  energy  losses  from  re-radiation  and 
convection. 

A typical  solar  tower  power  plant  is  shown  in  Figure  2. 
Heliostats,  comprised  of  two-axis  steerable  flat  reflectors  and 
ancillary  equipment,  are  situated  at  ground  level  about  a tower. 

The  solar  flux  incident  on  the  reflectors  is  continuously  directed 
at  a receiver  on  top  of  the  tower  where  it  is  converted  to  heat  in 
an  energy  transfer  fluid.  The  heat  is  transported  to  the  ground  and 
used  to  operate  a heat  engine  (15,16).  For  a nominal  100  MWe  power 


Figure  2.  Cross-sectional  view  of  solar  tower  system. 


plant,  the  overall  thermal  collection  efficiency  attainable  was 

estimated  at  about  60%  (17)  early  in  the  design  process.  The 

system  would  require  eight  426  x 426  m modules  each  containing 
2 

1840  - 37.2  m heliostats  and  a 122  m high  tower  with  a cavity  steam 
generator  receiver  of  approximately  94%  conversion  efficiency  (heat 
into  steam/energy  into  cavity)  (18).  The  land  area  requirement  here 
is  about  the  same  as  for  a linear  receiver  farm  of  similar  output. 

Although  seemingly  simpler,  the  solar  tower  design  has  some 
disadvantages  peculiar  to  itself.  These  are: 

(1)  onerous  tracking  requirements; 

(2)  flux  and  temperature  peaks  at  potentially  unacceptable 
levels  on  the  receiver; 

(3)  additional  costs  for  the  support  tower  and  heliostat 
mounts/ tracking  mechanisms. 

To  achieve  the  desired  concentration  factor  of  ca.  1000, 
the  total  optical  aberration  must  be  limited  to  2-4  milliradians 
(mr) ; this  requires  optical  flatness  of  the  reflector  surface  to 
ca.  2 mr,  and  tracking  accuracy  to  ca.  1 mr  (19).  Flux  distribu- 
tions on  the  receiver  surface  have  been  studied  by  Lipps  (20) , and 
by  Riaz  and  Gurr  (21),  among  others;  Blake  et^  al_.  (18)  also  consid- 
ered the  temperature  distribution.  If  all  the  heliostats  are  aimed 
at  the  central  zone  of  a receiver,  unacceptably  high  Gaussian-like 
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flux  peaks  result  (19)  with  the  currently  favored  design  of  a 
spirally-wound  hemispherical  boiler  and  a cylindrically-shaped 
receiver.  The  peak  flux  can  obviously  be  reduced  to  acceptable 
levels  by  displacing  the  aim  points  of  some  heliostats,  although  this 
does  increase  the  receiver  size  and  the  amount  of  bypassed  flux. 

Also,  the  nominal  design  solar  heat  flux  for  the  tower  receiver  is 
approximately  an  order  of  magnitude  greater  than  the  fluxes  encount- 
ered in  conventional  steam  boilers.  Thus,  design  operating  conditions 
and  fabrication  techniques  for  the  receiver  will  necessarily  depart 
from  conventional  practice.  The  generation  of  superheated  steam  by 
means  of  a thin-walled  receiver  of  tubular  construction,  using 
technology  garnered  from  aerospace  research  on  rocket  thrust  chambers, 
has  been  suggested  by  Sobin  ^ (22).  This  type  of  receiver,  with 

either  single-pass  or  helical-spiral  tube  arrangements,  could  readily 
* handle  the  anticipated  fluxes,  and  would  be  compact  as  well.  An 

alternate  design  uses  a liquid  metal  intermediate  fluid  which 
generates  saturated  steam  (23).  Liquid  metals,  such  as  sodium,  are 
known  to  possess  excellent  heat  transfer  properties. 

Yet  another  approach  foresees  a receiver  constructed  of  open- 
ended  ceramic  cavities  (24).  Heat  transfer  is  by  impingement- jet 
techniques;  the  fluid  is  made  to  flow  through  a porous  baffle  plate 
^ and  to  impinge  on  the  heated  receiver  surface.  Heated  air  at  ca. 

1000°C  is  supplied  to  a regenerative  Brayton  cycle  operating  at  a 
pressure  ratio  of  four.  The  ceramic  materials  are  resistant  to  the 
high  temperatures  encountered,  permitting  operation  in  an  oxidizing 
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atmosphere  at  temperatures  to  ca.  1650°C.  A quartz  window  over  the 
flux-inlet  side  of  the  cavity  offers  no  improvement  in  theoretical 
efficiency,  when  compared  to  a windowless  cavity  (25). 

Tower  costs  are  a function  of  both  height  (varying  approxi- 
mately as  height  squared)  and  receiver  weight.  A tapered  free- 
standing concrete  tower  is  the  preferred  choice  for  the  heavier  steam 
generating  receivers;  the  lighter  liquid  metal  or  air  cycle  receivers 
can  use  a more  economical  guyed  steel  construction  (19) . 

Heliostats  must  be  individually  controlled  about  two  axes  to 
maintain  the  solar  flux  on  the  receiver,  adding  appreciably  to  the 
cost  and  complexity  of  this  design.  Furthermore,  the  mean  mirror 
utilization  factor  here  is  only  ca.  0.70  (19).  This  includes  the 
averaged  effects  of  foreshortening  or  non-perpendicularity  of  mirrors 
to  the  solar  flux  (ca.  0.75),  shading  and  blocking  between  adjacent 
mirrors  (0.85-0.95),  and  bypassing  the  receiver  (0. 8-1.0).  Easton 
et  al.  (26)  developed  a procedure  for  optimizing  the  energy  collec- 
tion portion  of  a central  receiver  power  plant.  Still,  the  low 
utilization  factor  for  the  heliostats  largely  offsets  the  material/ 
cost  savings  that  this  design  otherwise  offers  relative  to  a solar 
farm,  since  the  mirror  utilization  factor  can  be  significantly 
higher  in  the  latter. 

Wen  and  Wu  (27)  present  a theoretical  performance  comparison 
between  the  central  tower  receiver  system  (solar  tower)  and  a distrib- 
uted system  (solar  farm  - linear  receiver) . These  systems  were 
analyzed  at  a rated  capacity  of  30  MWth,  delivered  in  the  form  of 
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superheated  steam  at  SSS’C  and  68  atm.  On  an  annual  basis  the  boiler 

section  (which  provides  22.7  Ml^th  at  rated  capacity)  of  a N-S  trough 

collector  has  an  energy  performance  factor  of  ca.  0.32  (where 
2 

1858  kWth-hr/m  -yr  = 1) , a two-axis  trough  shows  an  improvement 
to  ca.  0.51.  A tower  with  a conical  absorber  achieves  ca.  0.63, 

If  while  the  considerably  more  expensive  (26)  tilted  cavity  receiver 

» ‘ 

manages  ca.  0.85.  These  values  are  at  reflector  surface  slope 
error  = 0.20“  and  a pointing  error  of  0.16“,  and  decline  with 
increasing  errors.  But  the  land  utilization  factor,  defined  as 
the  ratio  of  the  collector  area  to  the  land  area,  varies  consider- 
ably with  system  employed.  The  central  receiver  system  has  a land 
utilization  factor  < 30%.  Higher  factors  exist  for  the  other 
types:  36%  for  the  two-axis  linear  collectors,  and  71%  for  a 
single-axis  trough  collector  oriented  N-S  (27) . These  would  tend 
to  even  out  land  requirements  for  a given  capacity;  reflector 
manufacturing  costs,  including  the  costs  associated  with  tracking 
mechanisms,  control,  maintenance,  and  auxiliary  components  such  as 
piping  and  insulation  (linear  systems)  and  tower  construction 
(central  receiver  systems)  must  also  be  considered  in  a comprehen- 
sive total  cost  evaluation  program. 

A preliminary  and  rather  simplified  cost  analysis  and  compar- 
ison of  a distributed  sys 'em  with  the  central  receiver  is  presented 
by  Caputo  and  Truscello  (28) . The  electrical  energy  cost  is  claimed 
to  be  48  mills/kWe-hr  in  the  case  of  the  central  receiver,  and 
approximately  double  this  for  linear  systems.  It  is  not  clear  if 
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all  the  aforementioned  factors  influencing  costs  of  these  systems 
were  taken  into  account  here,  but  the  inclusion  of  the  linear 
collector  capital  cost  at  twice  the  central  receiver  collector 
cost  (28)  makes  this  doubtful.  Clearly,  although  much  work  has 
recently  been  accomplished  on  solar  thermal  power  systems,  it 
appears  that  additional  work  is  necessary  before  one  system  can  be 
shown  conclusively  superior  to  others. 

Sandia  Laboratory,  Albuquerque,  NM,  is  currently  overseeing 
construction  and  testing  of  a 91  m tower,  with  a field  of  320  - 
1.2  X 1.2  m heliostats.  It  is  expected  to  yield  steam  at  1010°C, 
for  a maximum  output  of  ca.  5 MWth,  thus  making  it  the  largest 
solar  facility  in  the  world  (29).  Its  primary  purpose  is  to  test 
boiler  designs  and  materials,  mirtors,  and  other  components  for  a 
10  MW  facility  that  has  recently  been  approved  for  construction. 

The  site  for  the  larger  plant,  which  will  employ  conventional  steam 
generating  equipment,  is  on  130  acres  at  Barstow,  CA.  The  10  MW 
plant  is  expected  to  be  operational  in  1980  or  1981. 

D.  Linear  Radiation-Cavity  Collector 

Research  at  The  Pennsylvania  State  University  is  aimed  at 
developing  a concentrator-collector  of  the  solar  farm  type  that  is 
built  of  readily  available  materials,  is  simple  in  operation,  and 
will  heat  a gas  to  fairly  high  temperatures  (525-575°C)  so  it  can 
move  directly  into  a gas  turbine  in  a closed-cycle  Brayton  scheme. 
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To  this  end,  advantage  is  taken  of  the  ability  of  a small  aperture 
in  the  wall  of  a large  cavity  to  act  essentially  as  a blackbody; 
i.e.,  practically  all  the  radiation  incident  upon  it  is  absorbed 
(30) . Palmer  (31)  and  Palmer  and  Kuo  (32)  exploited  this  feature 
in  a proposed  modified  solar  farm  concept,  and  also  designed  an 
Archimedes  (i.e.,  Fresnel)  mirror-concentrator  to  replace  the  more 
complex  parabolic  trough  (see  Figure  3) . The  mirror  rotates  about 
the  collector  pipe,  which  remains  stationary  and  contains  a thin 
slab  of  commercial  graphite  as  flux  absorber.  In  the  original 
proposal  (32)  the  individual  mirror  strips  were  rotated  simulta- 
neously at  varying  rates,  concomitant  with  rotation  of  the  entire 
mirror  array.  In  the  present  work  a somewhat  simpler  tracking 
mechanism  is  employed. 

The  solar  farm  concept  examined  here  was  estimated  to  provide 

2 

an  output  of  25  MWe  from  an  area  smaller  than  1 km  (including  the 
power  plant) , under  what  seemed  to  be  realistic  assumptions  regarding 
losses  and  optical  features  (32).  The  total  conversion  efficiency, 
conservatively  calculated  relative  to  the  solar  flux  at  the  top  of 
the  atmosphere,  was  found  to  be  in  the  range  from  about  4. 5-5. 9%. 
Relative  to  the  flux  entering  the  collector,  the  efficiency  was 
estimated  to  be  from  12-14.5%. 

E.  The  Purpose  and  Method  of  the  Present  Study 

The  writer  believes  that  in  the  near  future  our  society  must 
develop  and  exploit  power  production  methods  such  as  those  described 
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above.  However,  choosing  among  the  numerous  methods  proposed  is 
difficult.  The  criterion  to  be  used  is  primarily  one  of  cost, 
which  is  inseparably  linked  with  overall  system  size,  complexity, 
reliability,  and  efficiency.  Idealized  mathematical  models  of 
proposed  designs  presented  in  the  past  and  discussed  above,  with 
simplistic  experimental  support  (if  any) , preclude  a rational  selec- 
tion process.  Thus,  experiments  must  be  performed  and  realistic 
models  constructed  to  determine  rather  accurately  the  potential 
efficiency  of  a full-scale  version  of  a particular  power-production 
method. 

To  this  end  a somewhat  simplified  first  model  of  the  linear 
cavity  receiver  design  described  in  Section  D,  above,  was  built  to 
approximately  1:11  scale  in  cross  section.  Experience  gained  with 
this  first  model  (described  in  Chapter  II)  led  to  the  design  and 
construction  of  a larger  (1:6  scale),  more  realistic,  and  more 
versatile  second  model,  representing  the  full-scale  device  much  more 
closely.  Detailed  studies  were  performed  (described  in  Chapter  III) 
on  the  temperature  rise  of  the  gas  as  a function  of  its  flow  rate 
as  it  traverses  the  receiver  pipe,  and  on  the  asymptotically 
limiting  gas  temperature  (i.e.,  stagnation  temperature  (1))  achiev- 
able at  a given  solar  flux  level.  Attendant  computer  modeling  of 
the  experimental  device  yielded  results  in  close  agreement  with 
empirical  (see  Chapter  IV) , and  provided  confidence  in  developing 
a computer  model  of  the  full-scale  device  (described  in  Chapter  V) . 
This  research  indicates  that  the  basic  concept  studied  here  has  proven 


CHAPTER  II 

FIRST  EXPERIMENTAL  MODEL  (STATIC  GAS) 

A.  Introduction 

An  experimental  model  of  the  full-scale  design  (see  Chapter 
V)  was  constructed  at  The  Pennsylvania  State  University  early  in 
1975.  The  cross-sectional  scale  factor  was  approximately  1:11. 

The  receiver  aperture,  which  was  not  insulated  in  order  to  allow 
the  solar  flux  to  enter  the  pipe,  was  proportionately  much  larger: 
ca.  0.6  of  the  receiver  pipe  diameter  versus  0.2  for  the  full-scale 
device  (see  Figures  4-6) . Insulation  also  was  much  thinner  (Table  1) , 
so  that  conduction  heat  losses  were  relatively  large.  The  short 
length  of  the  model  introduced  end  effects  - axial  temperature 
gradients  - which  channeled  heat  away  from  the  center,  lowering  the 
maximum  temperature  there. 

B.  Construction 

The  glass  receiver  tube  was  sealed  (rounded)  at  one  end  and 
capped  at  the  other  end  with  a Transite  plug.  This  had  two  small 
holes  drilled  through  it  (ca.  0.25  cm  diameter)  for  thermocouple 
insertion.  Graphite  absorber  material  was  placed  in  the  center  of 


igure  5.  View  into  aperture 


Figure  6.  Main  features  of  1/11  scale  experimental  model. 


Table  1.  Materials  description  and  properties. 


Receiver 

Absorber 

^“S^25'’C 

Pipe 


^100®C 

^^S^25°C 


Outside  dia.  x wall 
thickness 


Commercial  graphite,  0.34  cm  tk. 
0.87  (see  Chapter  IV) 

Corning  clear  1720  aluminosilicate 
glass 

0.90 

0.91  (assumed  same  as  1723  glass: 
1720  glass  not  listed) 


3.0  cm  X ca.  0.1  cm 


Aperture  L x W 
Material 

^100°C 

Insulation 

Thickness 

’^lOO-C 

Reflector 

Type 

Material 
Number  of  strips 
Strip  dimensions 

^^S^25°C 

Miscellaneous 

Thermocouple 

Temperature  readout 
instrument 

Heated  medium 


57  cm  X 1.9  cm  (at  the  glass) 

Johns-Manville  Transite  (i.e., 
asbestos  board) 

0.93  (see  Chapter  IV) 

Carborundum  Fiber f rax  (bulk  fiber) 
1.4  cm  (2.8  cm  in  some  tests) 
ca.  3.8  W cm  cm  ^ °C  ^ 


Archimedes 

3-M  back  aluminized  plastic 
22 

57  cm  X 1.6  cm 

0.67  (see  Table  2 and  Figure  7) 


22  ga.  Chromel-Alumel 

Hoskins  pyrometer,  type  AH 
Atmospheric  air 


Table  2.  Reflectance  of  plastic  mirror  strip.  Measured  at  25°C  with  a Beckman  DK-2A  spectrore- 


Figure  7.  Reflectance  of  plastic  mirror  strip  (measured). 


the  cube.  A thin-walled  outer  fiberglass  tube  with  a longitudinal 
aperture  supported  the  glass  tube  and  its  surrounding  insulation; 
the  ends  were  well  Insulated  also  to  minimize  axial  heat  losses. 

The  mirror  was  made  using  reflective  tape  with  an  adhesive 
backing,  mounted  on  0.3  cm  thick  Plexiglas  support  strips.  Screws 
and  locknuts  at  the  ends  fastened  the  strips  to  an  aluminum  angle 
frame;  this  frame  also  supported  the  receiver.  A second  metal 
frame  served  as  a base  for  the  device. 

Each  mirror  strip  had  to  be  individually  aligned  to  correctly 
reflect  the  solar  beam  onto  the  receiver  aperture.  The  alignment 
procedure  was  as  follows: 

(1)  The  receiver  tube  was  placed  in  its  holder  over  the 
mirror  array,  with  the  aperture  facing  the  central  axis  of  the  array. 

(2)  The  mirror  strips  were  rough-aligned  by  eye. 

(3)  The  receiver  was  removed  and  a string  placed  where  the 
central  axis  of  the  aperture  (in  the  plane  of  the  glass  window 
surface)  had  previously  been  located. 

(4)  Two  horizontal  strings  attached  to  heavy  structures  were 
placed  ca.  0.5  m apart  vertically  and  aligned  in  f vertical  plane  by 
means  of  a plumb  bob. 

(5)  The  device  was  placed  under  the  two  strings  in  such  a 
manner  that  these  were  centered  over  an  individual  strip  (to  be 
aligned).  The  array  was  locked  in  a horizontal  position. 
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(6)  A sighting  was  then  taken  across  the  string  which  had 
replaced  the  aperture  axis,  toward  the  strip  of  interest. 

(7)  The  strip  was  rotated  until  all  three  strings  appeared 

to  coincide  along  its  longitudinal  axis.  After  this  was  accomplished 
the  locknuts  holding  the  strip  were  tightened. 

Some  misalingment  was  present  as  a result  of  mirror  strip 
flexure  and  waviness  of  the  mirror  surfaces.  Hence,  the  actual 
width  of  the  beam  somewhat  exceeded  the  1.9  cm  window  aperture 
width:  the  measured  beam  width  in  the  aperture  plane  was  approximately 
2.5  cm,  a pointing  error  of  about  32%.  However,  not  all  strips  con- 
tributed equally  to  this  error;  it  was  chiefly  produced  by  those 
strips  farthest  from  the  array  center.  Therefore,  the  true  pointing 
error  flux  loss  probably  was  significantly  lower  than  the  above  value. 

C.  Experimental  Method 

The  small  experimental  collector  was  placed  on  the  roof  of  the 
Mineral  Industries  Building  of  The  Pennsylvania  State  University, 
where  it  was  tested  approximately  a dozen  times  in  the  period 
February-May , 1975  (Table  3).  Bright,  sunny  days  were  chosen  for 
the  test  runs.  Tests  generally  lasted  1-2  hours  near  solar  noon.  The 
data  collected  consisted  of  temperature  measurements  of  the  still  air 
within  the  collector.  Orientation  of  the  axis  of  the  device  was 
north-south  (N-S) . An  electric  clock  motor  was  employed  to  rotate 
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6 = declination  (l.e.,  angular  position  of  the  sun  at  solar  noon  with  respect  to  the  equator)  north 
positive.  Can  be  calculated  by  the  approximate  equation  of  Cooper  (36),  where  6 = 23.5  sin 
[360(284  + n)/365];  n = day  of  the  year. 
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the  entire  collector  at  15'’/h  in  order  to  track  the  sun.  However, 
this  motor  often  did  not  perform  adequately,  especially  early  in 
the  year  when  the  collector  tilt  from  the  horizontal  (of  the  longi- 
tudinal axis,  to  correct  for  solar  declination)  was  greatest.  It  was 
found  easier  to  adjust  the  collector  manually  then,  and  thus  maintain 
it  normal  to  the  incident  solar  rays,  the  preferred  orientation  of 
the  collector  axis  being  then  an  east-west  (E-W)  one.  This  orienta- 
tion required  fewer  adjustments  than  the  N-S,  over  a given  time 
interval  near  solar  noon.  The  thermocouple  was  normally  placed  in 
the  upper  semicircle  of  the  receiver,  opposite  the  aperture  (see 
Figure  6) . 

The  Pennsylvania  State  University  Meteorological  Observatory 
continuously  measures  the  total  solar  insolation  on  a horizontal 
surface,  with  an  Eppley  pyranometer.  For  use  in  the  experiments 
described  herein,  it  was  necessary  to  correct  the  insolation  data  for 
the  solar  angle  from  the  vertical,  and  then  to  determine  the  direct 
component  theoretically  (see  Chapter  IV  for  a full  discussion  of  the 
procedure) . 

D.  Experimental  Results  (see  Table  3) 

Maximum  (i.e.,  asymptotic)  temperatures  well  in  excess  of 
200‘’C  were  readily  attained  by  the  still  air  at  the  center  of  the 
receiver.  A plot  of  ths,  se  temperatures  versus  the  direct  solar  flux 
incident  on  the  concentrator  (for  a "standard  configuration"  without 


30 


extra  insulation  or  convection  covers)  is  shown  in  Figure  8.  Some 
scatter  is  evident,  but  the  trend  toward  higher  asymptotic  tempera- 
tures at  higher  flux  levels  is  apparent. 

While  normally  measurements  were  made  of  the  temperature  in 
the  upper  half  of  the  receiver,  placing  the  thermocouple  in  the 
lower  half  produced  a gain  of  ca.  10®C.  Most  of  this  rise  was  no 
doubt  caused  by  the  solar  flux  now  incident  on  the  thermocouple  Itself. 
Other  temperature  measurements  at  several  axial  locations  along  the 
receiver  indicated  a temperature  "droop"  of  ca.  yO'C,  as  the  ends 
were  approached. 

Figure  9 illustrates  the  time-dependent  temperature  variation 
and  the  fact  that  the  asymptotic  temperature  was  reached  in  approxi- 
mately 50  minutes  from  the  initial  exposure  of  the  device  to  the  sun. 
Although  for  the  test  run  illustrated  a double  thickness  of  insulation 
was  used  on  the  receiver,  the  results  shown  have  a general  validity. 

In  every  test  run,  irrespective  of  actual  conditions  of  solar  flux, 
insulation  thickness,  etc.,  asymptotic  temperatures  were  attained  in 
about  50  minutes. 

Varying  a number  of  experimental  parameters  had  the  following 
effects : 

(1)  Doubling  the  receiver  insulation  thickness,  from  1.4  cm 
to  2.8  cm,  resulted  in  a temperature  increase  of  about  10“C.  While 
the  added  insulation  obviously  reduced  conduction  losses  somewhat,  it 
also  decreased  the  available  flux  by  about  7%,  since  the  extra  insula- 
tion cast  a larger  shadow  on  the  mirror-concentrator. 
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(2)  Sealing  of  the  receiver-insulation  aperture  (see  Figure 
6)  with  a clear  polymethyl  methacrylate  (Plexiglas)  cover  strip,  ca. 
0.3  cm  thick,  reduced  convection  losses  and  increased  the  asymptotic 
temperature  by  O-IO^C.  The  effectiveness  of  the  cover  in  reducing 
convective  losses  appeared  to  depend  on  the  wind  velocity.  Since  the 
Plexiglas  material  grew  quite  warm  and  sagged  considerably  during 
testing,  it  seems  that  a significant  amount  of  infrared  flux  was 
being  absorbed.  The  I.R.  spectrum  of  the  material  indicated  high 
absorptance  in  a thickness  of  0.15  cm,  for  A 5 3 pm  (33).  When  a thin 
Saran  film  was  placed  over  the  aperture  instead,  severe  fogging  of 
the  film  occurred,  and  its  use  was  discontinued  after  a brief  trial. 

E.  Calculated  Results 

-2  -1 

A direct  perpendicular  solar  flux  of  800  W m s is  somewhat 
arbitrarily  selected  as  a typical  experimental  condition  from  Figure 
8.  Thermal  losses  are  assumed  to  consist  only  of  radial  ones,  and 
any  axial  losses  (end  effects)  are  ignored  here.  The  calculation  is 
based  on  a 1 cm  length  of  receiver  (at  its  axial  center)  purely  for 
ease  in  computation. 


The  total  concentrated  solar  flux,  per  centimeter  of  receiver 
length  is: 

• (Pjj)  (22  strips)  (1.6  cm/strip)  (1  cm  length) 
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At  the  steady-state  (i.e.,  asymptotic)  temperature  the  incoming 
flux  is  just  balanced  by  the  outgoing  thermal  losses.  An  energy 
balance  is  written: 


qg(l  - y)(1  - F3  • Pg)(l  - Pp)  = + e2F2A2a  • 


+ k • A_(T  - T )/(insul.tk.)  + h • A, (T  - T ) , 
3 g 00'  ' ' 4 g ® 


where  the  individual  terms  on  the  R.H.S.  of  the  equation  represent, 
respectively: 


(1)  radiative  loss  from  the  glass  window,  assumed  to  be  at 
the  temperature  of  the  gas; 


(2)  radiative  loss  from  the  aperture  sides  (see  Chapter  IV 
for  a discussion),  at  an  assumed  average  wall  temperature; 

(3)  conductive  heat  loss  through  the  insulation  (radial 

only) , where  the  temperature  of  the  insulation  inside  surface  is 

taken  = T ; 

g 

(4)  convective  heat  loss  from  the  aperture  (glass  + sides) , 
assumed  to  be  at  the  same  temperatures  as  in  Items  (1)  and  (2) 
above , 


and  the  term  on  the  L.H.S.  is  the  flux  actually  available  to  the 
receiver.  y»  the  pointing  error,  is  assumed  to  be  equal  to  0.5  of 
the  excess  beam  spread  at  the  aperture.  Therefore,  the  pointing 


error  is  set  equal  to  0.16  for  the  present  calculation  (see 
earlier  discussion).  Before  the  above  equation  can  be  solved, 
the  various  view  factors  (F's),  the  free  convection  coefficient 
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(h) , and  the  insulation  area  (A^)  must  be  evaluated. 


1 . View  Factors  (F's) 


The  relevant  receiver  dimensions  are  shown  in  Figure  10.  If 
one  makes  the  assumption  of  two-dimensional  radiative  interchange, 
the  "crossed-string”  method  (see  Chapter  IV  for  a full  discussion) 
can  be  used  to  evaluate  the  required  view  factors.  For  a long 
enclosure  (i.e.,  long  relative  to  width), 


Vc-D  • H , 


where  is  the  transverse  length  (width,  actually)  of  side  A,  etc. 
as  shown  in  the  accompanying  cross  section  of  a long  enclosure,  and 

Va-d 


Figure  10.  Receiver  dimensions. 
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• V-K-Ap.rture  " 2-')  " <2-3  + l-3)l/2 

or  F,  * 1.4/1. 9 * 0.74  . 

SlmlLrly,  1.3  • - (1.3  + 3.5  - 2.7)/2  - Fj  - 0.81,  and 

2-3  • "Oraphlta-Apottura  ' '<2.8  + 2.8)  - (1.6  + 1. 6)1/2  - F, 

= 0.43  . 

(This  slightly  underestimates  reflection  loss  from  the  graphite 
absorber;  that  loss  is  entirely  diffuse  and  three-dimensional). 

2.  Free  Convection  Coefficient  (h) 

h “ f (Gr,Pr)  in  free  convection.  Equations  are  available 
for  calculating  h.  First,  however,  Gr  must  be  calculated  to 
determine  if  the  external  flow  past  the  aperture  is  laminar  or 
turbulent  (at  an  assumed  wind  velocity  = 0) . The  relevant  Grashof 

3 

number  is  defined  here  as  Gr  * (g  6 L AT)/2,  where  AT  = T - T , and 

w ^ 

L » z,  the  axial  length.  Evaluation  of  the  gas  (air)  properties  at 

an  average  temperature  T^  = (T^  + T^)/2  corrects,  with  very  little 

error,  for  the  large  AT  (37).  B is  taken  to  be  equal  to  1/T  , 

00 

although  strictly  speaking  this  is  true  only  for  an  ideal  gas; 

however,  it  is  a reasonable  approximation  over  the  range 

0.5  < < 3.0  (38).  Assuming  that  the  receiver  (i.e.,  Pyrex  + 

graphite  + air)  is  isothermal  radially,  T = T = 235®C  at  a flux  of 

g w 
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800  W m ^ (see  plot  in  Figure  8).  Thus,  if  T = 0°C,  then  T /T  = 

(273  + 235)/273  * 1.86  - 8 = 1/273  = 3.66  x 10~^  °K~^,  and  = 

(235  + 0)/2  = 118°C  = 391°K.  At  this  T^,  v air  = 24.9  x lO"^  s“^, 

and  k . = 3.3  x lO"^  W m~^  °K~^,  and  Pr  = 0.69  (39).  Therefore, 

air  *' 

Gr  = (9.8  m s"^)(3.7  x 10“^  ‘’k”^)(28.5  cm  x 10“^  m cm"^)^ 

(235  - 0°K)/(24.9  x lO"^  s"^)^ 

—8 

= 3.18  X 10  and  the  flow  is  laminar  (40). 

For  an  inclined  surface,  the  local  Nusselt  number  Nu  becomes 

z 

Nu  = 0.508[Pr^/(0.952  + Pr)]°‘^^  [ (g  6 AT  cosa  z^)/v^]°‘^^  (40), 
z 

where  a = the  angle  from  the  vertical  and  is  taken  to  be  approximately 
50°.  Therefore,  at  the  axial  center  of  the  device. 


Nu  = 0.508[0. 69^/(0.952  + 0.69) (9.8  m s"^) 
z 


(3.66  X 10“^  °C“^)  (235°C) (cos  50°)  (28.5  cm  x 10~^  m cm“^)^/ 


„ ,.-6  2 -1,2,0.25 

(24.9  X 10  ms)] 
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Nu  = 4A,5,  and  since  Nu  = h • z/k  -*-11  =»  Nu  • k/z,  and  so 

z z z z z 

= (44.5) (3.3  X 10“^  W ®K~^)/(28.5  x lO"^  m)  = 5.15  W °k"^. 
3.  Insulation  Mean  Area 

For  the  insulation,  R /R  = (1.5  cm  + 1.4  cm)/1.5  cm  = 1.93; 

o i 

when  this  ratio  is  less  than  2,  the  mean  insulation  area  can  be 
approximated  with  A = (A^  + A^)/2  (Reference  40),  and  here 
A^  = tt(5.8  + 3.0)cm^/2  = 13.8  cm^,  for  360®  coverage.  Because  the 
aperture  occupies  an  approximate  75®  sector,  the  correct  area  now 
becomes  A^  = 13.8  cm^  (360®  - 75®)/360®  = 10.9  cm^. 

Upon  substituting  the  above  values  into  the  energy  equation, 
this  now  becomes 

1.89(1  - 0.16)11  ~ (0.43)(0.13)] (1  - 0.09)  W ^ 

0.9(0.74)  (5.67  x lO"®  W ®k“‘^) 

(235  + 273®K)^(1.9  x 10“'^  m^)  + (0.93)  (0.81) 

(5.67  X lO"^  W m"^  °K“^) (235/2  + 273®K)^(2.6  x 10"^m^) 

+ (3.8  X 10“'^  W cm“^  ®C“^)(10.9  cm^)  (235®C)/1.4  cm 
+ (5.15  W m‘^  ®k"^)[(1.9  cm^)(235®C)  + 

(2.6  cm^)  (235/2  ®C)  ] (lO”'^  m^  cm"") 

or 
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^ 1.36  W = (0.478  + 0.258  + 0.695  + 0.387)  W 

i 1.82  W . \ 

Clearly,  the  assumptions  of  radial  isothermicity , and  of 
skin  temperature  = are  overly  conservative.  It  will  now  be 
assumed  that: 

(1)  insulation  exterior  temperature,  T^) 

= + 50°C.  This  was  not  measured,  but  the  outside  of  the  insula- 

tion felt  warm  (i.e.,  the  temperature  there  was  > body  temperature, 
ca.  37°C). 

(2)  Tp  = Tg  over  the  insulated  portion  of  the  pipe;  at  the 
window,  Tp  = Tg  - 50'’C  (not  measured)  . 

Now  using  T = 235°C,  as  before,  and  combining  some  of  the 
S 

constants,  the  energy  equation  becomes: 

1.36  i (7.17  X 10"^^) (235  - 50  + 273)^  + (1.11  x lO"^^) 

[(235  - 50  + 50)/2  + 273]'^  + (2.96  x 10“^)  (235  - 50) 

+ (5.15  X 10"^)^{(1.9)(235  - 50  - 0)  + (2'.6)[(235  - 50  + 50)/2  -0  ]}. 

a:  h varies  little  with  temperature  and  has  not  been  recalculated. 
Its  accuracy  is  at  any  rate  probably  no  better  than  + 20%  (40) . 
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1.36  W = 0.32  + 0.26  + 0.55  + 0.34  - 1,35  W = 1.46  W . | 

i 

i 

The  result  is  acceptable,  considering  the  assumed  radial 
temperature  distributions,  pointing  error,  and  the  neglect  of  axial 
losses.  The  device  could  be  more  rigorously  modeled,  but  since  it 
was  only  the  precursor  of  the  desired  larger,  flowing-gas  experimental 
model  this  was  not  undertaken.  The  first  model  served  its  purpose 
in  verifying  the  flux-collection  potential  of  the  cavity  receiver 
design.  Furthermore,  the  chief  heat  loss  terms  had  been  Identified 
and  adequately  characterized. 

F.  Conclusions 

This  first  experimental  model  showed  that  moderately  high 
temperatures  were  achievable,  even  with  a very  small  version  of  the 
proposed  design.  The  combination  of  experiments  and  calculations 
indicated  that: 

(1)  The  dominant  loss  was  through  the  insulation,  so  this 
would  have  to  be  thicker  in  the  next  model. 

(2)  The  radiation  from  the  glass  window  was  large  and  might 
become  the  largest  loss  term  at  higher  temperatures.  It  would  be 
best  to  keep  the  aperture  as  narrow  as  possible. 

(3)  The  convective  loss  from  the  aperture  was  also  large  and 
presented  the  greatest  uncertainty  in  calculations.  Wind  speed  and 


direction,  aperture  shape,  collector  tilt,  and  the  amount  of 
receiver  rotation  affect  the  magnitude  of  the  loss  by  this  mode. 

It  is  difficult  to  quantify  the  influence  exerted  by  these 
parameters. 

On  the  basis  of  the  above,  it  was  decided  that  the  next  model 
should  be  of  approximately  1:6  scale,  with  a receiver  ca.  6 cm  in 
diameter,  an  aperture  no  wider  than  2 cm  at  the  glass  pipe,  ca. 

2.5  cm  of  insulation,  and  a reflector  of  40  x 1.6  cm-wide  mirror 
strips.  Very  simplified  calculations  with  these  parameters  indicated 
that  an  asymptotic  temperature  of  ca.  400°C  could  be  expected.  For 
a receiver  ca.  1.5  m long,  one  should  attain  a sufficiently  steep 
axial  temperature  gradient  for  detailed  study  of  heat  transfer  and 
collector  thermal  efficiency.  Studies  should  be  carried  out  at 
steady-state,  at  various  laminar  flow  rates. 
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CHAPTER  III 


SECOND  EXPERIMENTAL  MODEL  (FLOWING  GAS) 


A.  Introduction 


Work  with  the  small  model  discussed  in  Chapter  II  provided 
guidance  in  design  of  the  second  experimental  device.  It  was  clear 
that  the  model  should  be  substantially  larger  in  cross  section  in 
order  to  more  closely  simulate  a full-scale  system  and  in  order  to 
reduce  heat  loss  by  conduction  through  the  insulation.  Restrictions 
on  available  space  for  installation,  finances,  and  facilities  for 
fabrication  suggested  that  a reasonable  size  for  the  receiver  pipe 
would  be  about  5 cm  I.D.  with  a length  of  about  two  meters.  The 
system  should  be  designed  for  both  flowing-gas  and  static-gas 
conditions.  Also,  provision  for  gas  preheating  should  be  incorpor- 
ated into  the  system.  The  Archimedes  mirror-concentrator  should  be 
more  rigidly  constructed  and  should  have  a reliable  sun-tracking 
mechanism.  Measurements  on  the  system  should  at  least  include 
incident  solar  flux,  flow  rates,  and  the  temperature  distributions 
along  the  length  of  the  graphite  absorber,  the  gas,  and  the  external 
surface  of  the  glass.  Information  on  optical  properties  of  all 
materials,  mirror  reflectance,  and  thermal  conductivity  of  the 
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insulation  should  be  available  either  from  the  manufacturer,  the 
literature,  or  by  measurement. 

The  various  measurements  and  information  on  properties  should 
then  permit  development  of  a reliable  mathematical  model  for  describing 
the  performance  of  the  device,  and  subsequently,  of  the  full-scale 
design.  The  mathematical  model  would  serve  to  extend  the  necessar- 
ily limited  scope  of  experimental  study  and  allow  general  conclusions 
to  be  drawn. 

Clearly,  the  most  tractable  formulation  of  the  mathematical 
model  would  be  for  the  case  of  turbulent  fluid  flow  in  the  receiver. 

In  that  case,  forced  convective  heat  transfer  would  prevail  over  all 
other  modes  of  internal  heat  transfer,  and  radial  temperature  grad- 
ients would  be  sufficiently  small  to  be  of  no  consequence.  However, 

4 

the  difficulty  in  achieving  fully  turbulent  flow  (Re  > 10  (41)) 

within  such  a small  device  precluded  this  approach.  Very  high  flow 
rates,  with  a large  and  expensive  pump,  would  have  been  required. 
Furthermore,  turbulent  flow  would  produce  a rather  flat  axial  bulk 
temperature  profile,  undesirable  for  the  elucidation  of  the  effects 
of  various  experimental  parameters  to  which  the  profile  should  be 
sensitive. 

Hence,  although  it  was  anticipated  that  analytical  modeling 
would  be  more  difficult  for  the  laminar  flow  case  than  for  fully 
developed  turbulence,  it  was  decided  to  perform  experiments  over 
the  laminar  flow  range  only,  the  pump  and  recirculating  lines  being 


sized  accordingly. 


B.  Construction 


Approximately  one  year  (July,  1975-August,  1976)  was  required 
for  construction  of  the  model;  the  elapsed  time  covers  the  stages 
of  preliminary  design  to  completion  of  shake-down  runs  and  modifi- 
cations. In  constructing  the  model,  costs  and  other  constraints  led 
to  the  selection  of  readily  available  items  and  materials  over  more 
exotic  articles. 

The  model  has  been  designed  to  one-sixth  scale  (in  cross 
section)  relative  to  a system  considered  to  be  of  practical  dimensions 
for  power  generation  (32) . The  solar-thermal  concentrator-collector 
uses  an  Archimedes  (i.e.,  Fresnel)  mirror  concentrator,  5 cm  I.D.  x 
1.5  m long  glass  pipe  receiver,  0.34  cm  thick  graphite  absorber, 
and  gaseous  (Ar)  heat  transfer  medium.  An  alternative  0.09  cm 
thick  copper  absorber  with  a selective  oxide  surface  has  also  been 
employed.  The  absorbers  were  constructed  slightly  narrower  than 
the  pipe  I.D.,  and  simply  rested  within  it;  axial  support  was 
provided  by  a metal  plate,  of  the  same  thickness  as  the  absorber, 
which  extended  the  length  of  the  lower  (i.e.,  inlet)  end  piece. 

Tabs  on  the  plate  kept  the  absorber  immobile. 

The  mirror-concentrator  consists  of  40  suitably  tilted  reflec- 
tive strips,  each  1.6  cm  wide  by  1.5  m long;  the  strips  are  mounted 
with  a flexible  adhesive  on  metal  support  angles  in  a parallel 
planar  array  that  travels  E-W  within  two  1.67  m diameter  tracks 
located  near  its  ends.  The  glass  receiver  pipe  holds  the  thin 
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material  that  absorbs  the  concentrated  radiation  focused  on  it  and 
transfers  it  as  thermal  energy  to  the  flowing  gas,  which  exits  at 
high  temperature.  The  insulation  (2.5  cm  thick)  surrounding  the 
pipe,  except  for  a narrow  (2  cm  wide)  longitudinal  aperture,  is 
housed  within  a thin-walled  metal  tube  which  is  attached  to  the 
mirror  array  with  angles.  Thus,  the  insulation  rotates  in  such  a 
manner  as  to  always  maintain  alignment  between  the  aperture  and  the 
tracking  concentrator.  This  last  feature  is  one  of  the  major  strong 
points  of  the  design,  since  it  obviates  the  need  to  rotate  the 
receiver  pipe.  Also,  the  small  size  of  the  aperture  and  internalized 
flux  absorption  and  heat  transfer  within  the  cavity  receiver  reduce 
thermal  and  radiative  losses  substantially.  A 0.25  mil  (Mylar) 
film  over  the  aperture  suppresses  convective  losses  (see  Figures 
11-18  and  Table  4). 

The  concentrator-collector  assembly  is  tilted  41  degrees  (the 
latitude  of  State  College,  PA)  from  the  horizontal  so  that  the  axis 
of  the  collector  is  parallel  to  the  earth's  rotational  axis,  when 
the  device  is  oriented  N-S.  Mirror  array  tracking  (and  rotation 
of  the  receiver  insulation)  is  by  means  of  a very  accurate  device, 
a high-torque  stepping  motor  which  takes  400  steps/revolution.  This 
motor  is  connected  through  a chain  drive  to  the  array;  a 30:1  reduc- 
tion worm/worm  gear  combination  multiplies  torque  and  inhibits  array 
motion  with  the  motor  stopped.  Control  (no  feedback  loop)  is  by  a 
continuously  variable  solid-state  timing  circuit. 
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Figure  11.  1/6  Scale  experimental  model  - side  elevation 
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P’lgiire  14.  Drive  detail 
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Figure  15.  Flux  bulb. 
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Figure  17.  Kquipinent  sclicMnatic 
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Table  4.  Materials  description  and  properties. 


Receiver 
Glass  tube 


^25°C 

25°C 
^ 250“C 
Flanges 


Flange  inserts 


Gaskets  (standard) 


Corning  type  7740  clear 

borosilicate  (Pyrex) , with 
conical  ends 

50  mm  I.D.  x 5 mm  wall  x 1829  mm 
long 

#72-7502 

1.09  W m~^  °k“^ 


Aluminum  alloy,  two  required 
50  mm 

Corning  #72-9062 

Hard,  bonded  asbestos,  two  required 
50  mm 

Corning  #72-9058 

Asbestos  type  A,  two  required 
50  mm 

Corning  #72-9130 


Gaskets  (high  temp.) 
Thickness 


Silicone  rubber 
0.3  cm 


Insulation 


Amount  used 
^ 100°C  (35,43) 
(Ps)25''C 

Insulation  support  tube 


Carborundum  ceramic  fiber 
(Fiberfrax)  blanket 

2.5  cm  X 96.1  kg  m ^ 

1 2 
ca.  1 m 

3.8  W cm  m“^  °C~^ 
ca.  1.00^ 


Aluminum 

11.4  cm  O.D.  X 0.09  cm  wall 


Aperture  L x W 
Sides 


152.5  cm  X 2 cm 

Johns-Manville  Transite  (i.e., 
asbestos  board) 
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Table  4,  Continued. 

Receiver 

Thickness 

^200°C 
End  pieces 
Size 

Absorber  support 
(inlet  only) 


0.32  cm 
0.94^ 

Copper  (outlet);  st.  stl.  (inlet) 

ca.  5 cm  O.D.  x 0.15  cm  wall  x 
46  cm  long 

0.34  cm  thick  stainless  steel 
with  tabs 


Convective  cover 
Thickness 

25°C 

Absorber  (standard)  (45) 

(Ps)25'’c 

^“S^25°C 

^250'’C 

H50°C 

Absorber  (selective) 
Substrate  material 
Thickness  ' 
Selective  surface 
Thickness 

CPs')25°C 

(ctg)25‘’C 

“IR 

Pre-heat  (inlet) 

Catalog  no. 

Capacity 

Size 

Control 


DuPont  Mylar 
0.25  mil 
0.85^ 

Airco-Speer  mold-stock  graphite, 
commercial  787G  Grade 

0.13^  (nearly  100%  diffuse) 

1 - 0.13  = 0.87 

0.87^ 

173  W m ^ ®C  ^ (with  grain) ^ 

Commercial  copper 
0.09  cm 
CU2O  + CuO 
Not  known 
0.11® 

1 - 0.11  = 0.89 

0.20 

Heating  tape  (electrical  resistance) 
Fisher  Scientific  #ll-463-51D 
768  W 

2.5  cm  X 244  cm 
V'ariac 
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Table  4.  Continued. 

Receiver 

Heat  transfer  fluid 
Purity 

Thermocouples  (external) 
Material 
Gauge 

Thermocouples  (internal) 
Material 
Sheath 
Type 

Compression  fitting 
Thermocouple  switching 

Temperature  readout  instrument 
(dial) 

Type 

Accuracy 

Internal  resistance 

Temperature  readout  instrument 
(digital) 

Model  No. 

Accuracy  (rated) 

Resolut ion- repeat ability 

Lead  resistance 

Reflector  - concentrator 
Mirror 
Size 
Number 

^^S^25°C 


Argon 

Standard  commercial  grade 
Type  K 

Chromel-Alumel 

28 

Type  K 

Chromel-Alumel 

Stainless  steel,  0.16  cm  O.D. 

Omega  Engineering  Cat.  No. 
CASS-116G  and  CASS-116E 

Omega  Cat.  No.  BRLK-116 

2 rotary,  1 toggle,  std. 
commercial 

Hoskins 

AH 

Note  (b) 

50  ohm 

Precision  Digital  Corporation 

7701-K 
+ 2°C^ 

1°C 

To  1000  ohm  for  < 1°C  error  • 

PPG  back-silvered  float  glass 
1.6  cm  X 153  cm  x 0.32  cm  tk. 
40 

0.80^,  direct  only 
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Table  4.  Continued. 

Reflector  - concentrator 

Supports  2.5  cm  x 1.3  cm  x 0.3  cm  tk.  A1 

angles 

Mounting  Screws  and  locknuts 

Adhesive  Silicone  sealant 

Array  frame  Aluminum  angle 

Size  5 cm  X 5 cm  x 0.8  cm  tk. 

Wheels  5 cm  dia.  Bassick  /*9389  casters 


Track 

Size 

Dia.  (rolled  to) 


Aluminum  tee,  two  required 
10  cm  X 7.6  cm  X 1.0  cm  tk. 
1.67  m 


Main  frame 
Size 

Length  required 
Wheels 


Slotted  steel  angle 

5.7  cm  X 3.8  cm  x 0.16  cm  tk. 

ca.  60  m 

7.5  cm  dia.  Bassick  #3619  casters 


Tracking 

Motor 

Type 

Torque 


Motor  control 
Type 

Drivetrain 

Sprockets 

Size 

Shaft 


Superior  Electric  Company  stepping 
motor 

#M111-FD327,  capable  of  200  or  400 
steps/revolution 

28.8  kp  cm  at  slow  stepping  speed 
and  for  400  steps/rev  (i.e., 
0.9°/step) 

Superior  Electric  Company,  solid 
state 

#STM1800DV  translator  module;  also 
RC  timing  pulse  circuit 


Boston  Gear  No.  80,  steel,  two 
required 

2.5  cm  pitch,  12  tooth 

1.9  cm  dia.  x 2 m long  aluminum, 
with  1.6  cm  dia.  ends 
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Table  4,  Continued. 


Trackinfi 


Bearings 


Chain 


Gears 


Worm  gear 


Fluid  flow  system 


Capacity 

Max.  press,  (gauge) 


Motor  size 


Return  line 


Flow  measurement 


Manometer  fluids 


Absolute  pressure 


Nozzle  AP 


Filter 


Dessicant 


Flux  measuring  device 


Absorber 


Boston  Gear  Cat.  No.  06904-3H 
ball  brg.  pillow  blocks,  1.6  cm 
bore,  5 required 


Boston  Cat.  No.  C2042D,  2.5  cm  pitch 


Roller,  with  oversize  Delrin 
rollers 


Boston  Cat.  No.  G1061  cast  iron,  30 
tooth,  R.H.  thread 


Boston  Cat.  No.  GH1066  soft  steel 


Graphite  vane,  centrifugal 


Gast  model  #0522-V103-G18D  oil-less 


ca.  100  liters/minute 


ca.  1 atm. 


1/4  HP 


ca.  1 cm  dia.  copper 


ca.  2.5  cm  dia.  constricting  to 
ca.  0.5  cm  dia.  nozzle 


n-butyl  phthalate  (p  of. 
1.042  g cm-3) 


Whatman  Gamma-12  unit  with  R40 
filter  tube 


1 mm  thick  Pyrex 


10  cm  copper  plate 
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Table  4.  Continued. 


Flux  measuring  device 


(dg)  25‘’C 
^25°C 


Support  rod 
^25®C 


Measured  here 
^Discussed  in  text 


0.39^ 

0.03 

0.208  dia.  x 10.3  cm  long  copper 
3.89  W cm~^  '’C'^ 


« 
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The  overall  insulated  length  of  the  receiver  pipe  is  approxi- 
mately 2.7  m.  This  length  includes  a ca.  0.6  m long  piece  at  each 
end  which  serves  to  minimize  end  effects  and  has  360°  of  insulation 
(i.e.,  no  window  aperture).  The  window  aperture  occupies  153  cm  of 
length  in  the  mid-section.  Of  this,  only  the  91  cm  long  portion 
beginning  at  32  cm  and  extending  to  123  cm  (the  distances  refer 
to  positions  measured  in  the  same  direction  as  the  flow,  from  the 
beginning  of  the  aperture  - refer  to  Figure  17)  is  normally  studied. 
The  mirror  length  is  also  153  cm;  hence,  a shadow  at  either  the  inlet 
or  outlet  end  of  the  aperture  is  present,  its  length  a function  of 
the  solar  altitude  (i.e.,  day  of  the  year).  With  the  device  tilted 
at  a fixed  angle  of  41°,  this  shadow  is  longest  at  the  solstices, 
measuring  ca.  25  cm  at  those  times.  The  32  cm-123  cm  central 
section  remains  fully  illuminated  during  an  experiment  at  any  time 
of  the  year. 

The  pump  used  to  circulate  the  fluid  cannot  operate  at 
temperatures  much  in  excess  of  ambient;  thus  the  fluid,  after 
traversing  the  exit  flow-smoothing  end  piece  (which  is  fully  insu- 
lated) , leaves  at  high  temperature  and  is  made  to  pass  through  an 
uninsulated  copper  return  line,  to  re-enter  the  pump  at  near-ambient 
temperature.  This  restriction  on  receiver  inlet  temperature  is 
removed  by  means  of  an  electrical  preheat  section  (770  W heating 
tape  wrapped  around  the  metal  inlet  end  piece) , which  allows  simula- 
tion of  a longer  receiver,  under  gas  flow  conditions. 

The  experimental  device  was  located  on  the  roof  of  the 
Mineral  Industries  Building  at  The  Pennsylvania  State  University. 
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Because  it  was  not  fully  weatherproofed,  a shed  was  specifically 
constructed  for  it  (part  of  which  is  visible  in  Figures  11  and  12) , 
in  which  it  was  housed  when  not  being  tested.  On  a test  day,  the 
device  was  wheeled  out  of  the  shed  and  oriented  with  its  longitudinal 
(i.e.,  flow)  axis  N-S. 

In  actual  fabrication/construction  of  the  device,  three  main 
constraints  have  been  adhered  to,  since  they  are  crucial  to 
performance: 

(1)  The  maximum  tolerable  concentrator  array  pointing  error 
has  been  set  at  10%  of  the  aperture  width  of  2 cm,  or  0.2  cm.  With 
the  array  55  cm  below  the  aperture,  the  corresponding  angular  devia- 
tion 0 is 


6 = tan"^(0.2  cm/55  cm)  = 0.21' 

(this  assumes  perfect  mirror  optics)  and  the  permissible  error  in 
array  tilt,  for  ca.  100  cm  array  width  between  the  supporting  wheels, 
is  then 


y/100  = tan(0.21')  y = 0.36  cm  . 

Therefore,  the  tracks  have  been  constructed  with  a tolerance  of 
+ 0.18  cm  on  the  83.8  cm  track  radius  in  order  to  maintain  the 
pointing  error  within  the  specified  limit. 
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(2)  Infiltration  of  air  into  the  receiver  must  be  prevented; 
dilution  of  the  heat  transfer  fluid  would  affect  its  heat  capacity 
and  transfer  characteristics.  Furthermore,  the  possibility  of  oxida- 
tion damage  to  receiver  internals  at  the  high  temperatures  encountered 
during  experiments  must  be  avoided.  Thus,  the  system  is  hermetically 
sealed  from  atmospheric  air  and  argon  is  employed  as  the  working 
fluid,  with  the  choice  of  fluid  based  on  its  inert  nature,  reason- 
able cost,  and  ease  of  preventing  leaks. 

(3)  Temperature  measurements  of  the  fluid  must  be  representa- 
tive of  its  true  energy  level.  For  this  purpose,  a number  of 
thermocouples,  installed  through  compression  fittings  at  both  metal 
end  pieces  of  the  receiver,  extend  to  five  pre-selected  measuring 
positions  inside  the  glass  receiver  (see  Figure  19  and  Table  5) . 

Some  difficulty  has  been  encountered  in  meeting  the  second 
and  third  criteria.  Gas  leakage  at  the  pump  cannot  be  eliminated, 
so  a plastic  envelope  with  provision  for  feeding  in  a trickle  of 
excess  argon  has  been  constructed  around  it.  This  arrangement 
functions  well,  at  some  additional  cost  in  argon  usage.  The  hard 
graphite  vanes  of  the  pump,  in  rubbing  on  the  inside  of  the  pump 
chamber,  gradually  grind  away.  The  resulting  very  fine  particulate 
matter  is  readily  entrained  into  the  flowing  gas  stream.  Although 
the  pump  is  provided  with  a felt  filter  (50  pm),  this  has  proven 
ineffectual  in  preventing  particulate  carry-over  into  the  receiver, 
where  it  coats  internals.  A 2 pm  filter  has  been  installed  in  the 


Figure  19.  Thermocouple  location. 


Table  5.  Thermocouple  description. 


No. 

Type^ 

Position 

Temperature  Measured 

T1 

E 

Inside  flux  bulb 

Copper  plate 

T2 

E 

Outside  flux  bulb 

Heat  sink 

T3 

E 

2^^ 

Insulation  exterior 

T4 

E 

2 

Pyrex  exterior 

T5 

E 

2 

T^,  of  gas 

T6 

G 

2 

Absorber 

T7 

E 

74 

Pyrex  exterior 

T8 

E 

151 

Insulation  exterior 

T9 

E 

151 

Pyrex  exterior 

TIO 

E 

151 

T^,  of  gas 

Til 

G 

151 

Absorber 

T12 

E 

Nozzle 

T^,  of  gas 

T13 

E 

Thermometer 

Ambient  at  collector 

T14 

E 

32 

T^,  of  gas 

T15 

G 

32 

Absorber 

T16 

E 

74 

T^,  of  gas 

T17 

E 

123 

T^,  of  gas 

T18 

G 

123 

Absorber 

T19 

E 

74 

Pyrex  at  aperture 

= exposed;  G = grounded 

°axial,  from  beginning  of  aperture  at  inlet  end,  cm  (see  Figure  19) 
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line  (pressure  drop  across  it  is  acceptable) , and  captures  essen- 
tially 100%  of  the  entrained  particulates.  Also,  the  binder  in  the 
asbestos  gasket  at  the  inlet  metal/glass  joint  tends  to  vaporize 
at  higher  temperatures  (e.g.,  when  the  heating  tape  is  used  to 
preheat  the  incoming  flow)  and  has  been  replaced  with  silicone 
rubber.  This  can  sustain  the  highest  temperatures  (probably  > 400°C 
at  the  joint)  without  damage. 

The  preliminary  design  of  the  apparatus  had  only  four  thermo- 
couples (two  for  the  fluid,  two  for  the  graphite)  placed  inside  the 
receiver.  These,  being  located  very  near  the  receiver  ends,  gave 
temperature  readings  that  were  strongly  influenced  by  end  effects. 
After  a few  experiments,  five  additional  internal  thermocouples, 
extending  farther  into  the  receiver,  were  installed.  Even  with 
these,  there  is  some  question  as  to  the  representativeness  of  the 
measurements  obtained.  A full  treatment  of  this  complex  topic  is 
given  in  Chapter  IV.  All  thermocouples  have  been  fabricated  and 
installed  to  comply  with  common  practice  (49,50). 

To  assure  impingement  of  the  concentrated  solar  flux  on  the 
graphite  absorber,  the  absorber  was  initially  constructed  x-shaped 
in  cross  section.  Thus,  the  receiver  pipe  was  divided  into  quad- 
rants. Problems  were  encountered  with  unequal  heating  in  individual 
quadrants  and  large  transverse  heat  conduction  losses.  A single 
flat  absorber  was  subsequently  installed,  bisecting  the  receiver 
into  two  semicircles.  It  functions  very  satisfactorily. 

A few  experiments  have  been  made  with  a selectively  surfaced 
absorber  replacing  the  graphite.  The  selective  surface,  an  oxide 
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semiconducting  layer,  was  formed  on  a brightly  polished  commercial 


Cu  sheet.  A number  of  oxidizing  solutions  of  different  concentra- 


tion, plus  varying  exposure  (immersion)  times,  were  tried  on  small 


Cu  samples  in  order  to  optimize  the  selective  properties.  Adhering 


to  the  recommendations  of  Christie  (46)  resulted  in  what  appeared  to 


be,  visually,  the  best  surface:  a slightly  glossy  gray-blue  surface. 


This  surface,  which  had  the  highest  a„  measured  among  samples,  was 


produced  by  immersing  a 5 cm  x 120  cm  copper  strip  (two  such  were 


used  in  the  receiver  at  one  time)  in  2500  ml  of  solution  containing 


300  g NaOH  and  150  g NaC102,  maintained  at  60°C  and  continuously 


stirred,  for  10  minutes. 


C.  Material  Properties 


The  properties  of  the  relevant  materials  employed  in  the 


collector  have  a decided  influence  on  its  performance.  Knowledge  of 


them  is  a prerequisite  for  an  understanding  of  the  optical  and  thermal 


processes  in  the  collector.  These  properties  have  been  either  measured 


here  or  taken  from  the  literature. 


1.  Optical  Properties 


a . Solar  range  (ca.  300-2500nm) 


On  a clear  day,  at  the  bottom  of  the  atmosphere,  nearly  99% 


of  the  sun's  energy  spectrum  lies  between  300  and  2500  nm  (51,52) 


Optical  properties  in  this  range  can  be  readily  measured  with  a 


Beckman  DK-2A  spectroref lectometer . This  instrument  does  not  measure 
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the  absolute  transmittance  (x)  or  reflectance  (p) , but  compares  the 
sample  with  an  air/MgO  or  MgO  reference  standard.  Because  the  MgO 
reflectance  is  very  high,  on  the  order  of  0.97  (48,53,54),  it  is 
common  to  assume  that  it  represents  a perfect  reflector,  and  use 
an  MgO  versus  MgO  calibration  curve  as  representing  total  (i.e., 

100%)  reflectance  (53) . The  measured  sample  properties  are  then 
corrected  for  any  small  departures  of  this  curve  from  the  100%  line 
on  the  graph.  Both  total  and  diffuse  x and  p can  be  measured,  but 
only  at  room  temperature.  The  variation  of  these  properties  with 
temperature  will  be  discussed  in  the  next  section. 

Solar  irradiance  (defined  as  the  solar  flux  on  a unit  area 

-2 

of  surface,  at  the  angle  of  incidence  of  the  radiation,  in  W m ) 
varies  with  wavelength,  and  thus  a further  correction  of  properties 
is  necessary  to  obtain  a value  applicable  in  solar  experiments. 
Several  irradiance  curves  have  been  determined;  older  works  (55,56) 
have  u^'ually  been  based  on  the  data  of  Moon  (57).  Thekaekara  (51,52) 
has  more  recently  done  extensive  experimental  and  theoretical  work 
in  this  area.  Solar  irradiance  spectra  generated  by  a computer 
model  of  '"he  atmosphere  are  available  from  NASA  (52),  and  are  used 
here. 

The  irradiance  is  a function  of  the  air  mass  (AM) , atmospheric 
ozone  and  precipitabie  water  vapor  content,  the  turbidity  coefficient 
6,  and  the  wavelength  exponent  a.  A U.S.  standard  atmosphere  with 
20  mm  of  precipitabie  H^O  and  3.4  mm  of  0^  is  generally  assumed  to 
be  present  (51,52).  The  air  mass  is  the  ratio  of  the  mass  of  air 
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in  a column  with  unit  cross  section  along  the  path  of  the  solar 
rays,  to  the  mass  in  a vertical  column  extending  to  the  top  of  the 
atmosphere.  For  moderate  z (i.e.,  where  no  account  has  to  be  taken 
of  the  curvature  of  the  solar  ray  due  to  refraction  in  an  increas- 
ingly denser  atmosphere) , 


AM  = sec  z , (51) 

where  z is  the  zenith  angle  (the  angle  between  the  vertical  at  a 
point  on  the  earth  and  the  line  joining  that  point  to  the  sun) . 
Here,  given  the  41°  latitude  of  State  College,  PA,  and  the  fact 
that  most  experiments  have  been  performed  near  the  equinoxes,  an 
AM  = 1.5  is  assumed  throughout. 

For  atmospheric  turbidity  effects,  Thekaekara  (51)  uses 
Angstrom's  equation.  A.  K.  AngstrSra  (58)  developed  an  equation  for 
the  optical  thickness  due  to  aerosols  in  terms  of  a and  S,  as 

t^(A)  = SA"“  , 

where  the  optical  thickness  is  defined  as  the  number  of  mean  free 
paths  of  a radiation  photon.  The  exponent  a is  influenced  by  the 
atmospheric  particle  size  distribution.  Extensive  observations 
indicate  that  a good  average  value  is  a = 1.3  + 0.2  (58).  The 
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coefficient  B is  a measure  of  the  turbidity,  or  haze  content,  of 
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the  atmosphere  (59) . It  varies  widely,  from  several  hundredths 
(clear  atmosphere)  to  several  tenths  (low  atmospheric  transparency) 
(60).  Values  of  a = 1.3  and  6 = 0.02,  signifying  a relatively  clear 
atmosphere  (51),  are  selected  as  appropriate  parameters.  The  inte- 
grated direct  perpendicular  solar  irradiance  resulting  from  the 

_9 

values  selected  is  870  vV  m “ (52) . An  examination  of  the  measured 
irradiances  during  experiments  performed  here  shows  that  this  is  a 
reasonably  good  average  value  for  a clear  day  in  State  College,  PA. 
Therefore,  this  irradiance  spectrum  is  used  in  the  present  work 
(see  Table  6) . It  might  be  noted  here  that  the  solar  spectral 
energy  partition  is  not  very  sensitive  to  departures  in  a and  S 
(and  thus  the  integrated  irradiance)  from  the  chosen  values  (51) , 
and  it  is  this  spectral  partition  which  is  of  interest  in  determining 
averaged  optical  material  properties  in  the  solar  spectral  range 
(see  Table  6) . 

b.  Thermal  infrared  (TIR) 

An  instrument  for  measuring  optical  properties  at  X > 2.5  ym 
was  not  available  at  The  Pennsylvania  State  University.  With  one 
exception  (the  graphite  absorber)  the  values  to  be  discussed  are 
taken  from  the  literature.  Quite  often  the  published  values  are 
not  very  useful  because  either  full  information  on  sample  and  source 
temperatures  is  not  presented,  or  the  sample  temperature  is  much 
different  from  the  anticipated  application  temperature. 

For  radiative  heat  exchange,  two  parameters  are  needed  to 
specify  the  desired  material  property.  One  is  the  temperature  of 
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Continued 


Corning  data  (42) 
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the  source  (or  the  wavelength  corresponding  to  the  maximum 
emissive/absorptive  flux  of  a blackbody  at  that  temperature,  as 
determined  by  Wien's  displacement  law).  Since  radiative  properties 
of  a given  material  can  vary  appreciably  with  temperature  of  the 
material  (not  to  be  confused  with  the  wavelength  dependence 
mentioned  above) , the  temperature  of  the  material  in  question  is 
also  important. 

A consistent  nomenclature  will  be  adopted  here;  thus 


(e_), 

a b 


or 


a b 


a b 


will  be  taken  to  signify  a source  temperature  or  wavelength  "a" 

(or  the  symbols  S and  IR  will  be  used  to  denote  a property  averaged 
over  the  solar  or  infrared  wavelengths,  respectively),  while  "b" 
will  correspond  to  the  temperature  of  the  material  in  question  at 
which  the  property  was  actually  measured.  If  b = a,  only  one 
subscript  will  be  shown.  Also,  unless  specified  otherwise,  direc- 
tional dependence  will  be  assumed  negligible  and  the  properties 
measured  perpendicular  to  a surface  will  be  used,  when  available. 

c.  Data 

Neuroth  (61)  reports  little  change  in  the  infrared  trans- 
mittance (X  > 2.5  ym)  of  borosilicate  glass  at  high  temperature, 
while  Markin  and  Sobolev  (62)  indicate  that  the  infrared  reflec- 
tances of  02^3  Si02  (the  main  constituents  of  Pyrex  glass) 
remain  unchanged  to  1000'’C.  Touloukian  (53)  shows  a near-constant 
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transmittance  of  Pyrex  to  800®K;  thus,  the  value  measured  here  at 

25°C  will  be  used  as  applicable  to  the  much  higher  temperatures  of 

the  receiver  pipe.  The  manufacturer  measured  r „ = 0.89  in  one 

25®C 

instance  (34),  but  recommends  e = 0.94,  without  specifying  any 

temperature  or  wavelengths,  in  another  booklet  (63).  The  sample 

thickness  is  not  stated  in  either  case. 

Touloukian  (53)  reports  p =0.84  for  a 1.3  cm-thick 

250'’C 

sample,  with  almost  no  temperature  dependence  of  e to  900°K. 

Pirani  (64)  measured  Pyrex  emittance  in  the  250-800°C  range  and 
gives  = 0.94  (no  sample  thickness  stated,  but  6.9  mm  quartz 

standard  used),  again  little  affected  by  higher  temperatures. 
Extrapolation  of  the  manufacturer's  data  on  transmittance  yields 
e = 1.00  for  a 5 mm  thick  sample  (63).  Because  of  the 

uncertainty  involved,  and  since  the  aperture  actually  radiates 
from  a width  somewhat  greater  than  its  nominal  width  of  2 cm  (the 
aperture  sides  do  not  quite  touch  the  glass  pipe,  thus  allowing  a 
larger  area  to  radiate  than  is  available  for  the  entering  flux) , a 
conservative  ^250°C  ~ will  be  assumed.  The  250°C  temperature 

above  has  been  selected  as  an  average  value,  representative  of  the 
entire  collector  pipe  length,  for  all  experiments. 

The  reflectance  can  be  calculated  from  the  well-known  modified 
Fresnel  equation  (65) 

p = [(n  - l)/(n  +1)]^  , 


where  n is  Che  index  of  refraction.  Pyrex  has  n = 1.474  (63), 

which  can  be  considered  a constant  to  750°C  (66),  and  hence 

^^S^  250‘’C  ” 0.037/interface,  or  ca.  0.07  for  both  surfaces.  Thus, 


S"250°C 


=l-T-p=l-  0.85  - 0.07  = 0.08.  Only  the  reflected 


flux,  0.07  of  the  concentrated  solar  flux,  is  taken  as  lost;  Che 
absorbed  component  (0.08)  is  included  in  the  energy  balance  on  the 
receiver  (see  Chapter  IV) . 

The  reflectance  of  the  insulation  material  is  very  high,  about 
equal  to  that  of  the  MgO  reference  standard,  and  nearly  100%  diffuse. 

It  is  assumed  to  be  temperature-independent  over  the  range  of  interest, 
i.e.,  (P2)  250°C  ~ characteristic  temperature  is  again 

taken  as  250®C.  Singham  (44)  lists  the  integrated  hemispherical 
emittance  of  asbestos  paper  (the  e of  asbestos  board  being  nearly 
identical  at  other  temperatures,  but  not  shown  for  temperatures  of 


interest  here)  as  e 


20Q°C 


0.94,  and  = 0.93. 


The  graphite  absorber  material  presents  more  difficulties. 

The  measured  (°‘3)25°C  0.87  is  in  good  agreement  with  the  value 

^'^6000'^K^25°C  * *^‘85  measured  by  Sieber  (67)  for  a graphite  material, 

but  is  lower  than  the  (a  ) = 0.91-0.92  of  "7087"  graphite  in 

^ 25  C 

Touloukian  (53).  With  a few  exceptions,  most  graphites  exhibit 
optical  properties  that  are  nearly  independent  of  wavelength  and  of 
sample  temperature  to  at  least  1000°C  (53,68,69).  Data  on  the  tempera- 
ture dependence  of  for  the  787G  graphite  used  here  are  apparently 
not  available.  It  will  therefore  be  assumed  that  25®C  ” 0.87 

is  valid  at  250'’C,  Che  assumed  typical  operating  temperature.  Because 
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the  view  factor  from  the  graphite  to  the  aperture  (i.e.,  the  out- 
side) is  so  small,  any  inaccuracy  in  the  absorptance,  and  hence  the 
reflectance,  will  have  little  bearing  on  calculated  results. 

The  infrared  emittance  of  Speer  "7087"  graphite  has  been 
measured  (53,68)  (Speer  does  not  manufacture  a 7087  graphite;  787G 

was  probably  intended  (70)).  The  reported  value,  e o-  = 0.85-0.9 

23U  L 

(53,68)  is  also  valid  for  some  other  graphite  types  (68,69).  A 

sample  of  the  absorber  was  sent  to  the  W.  B.  Meinel  Laboratory 

(71),  which  tested  emittance  to  sample  temperatures  of  200°C,  and 

over  2 < X < 15  um.  The  results,  which  indicate  strongly  selective 

properties  (i.e.,  high  a and  low  , conflict  both  with  extrapo- 

b IR 

lated  trends  of  the  measurements  made  here  and  also  with  the 
literature  cited  (except  Sieber  (67),  and  even  there  the  selectivity 
is  much  less) . 

The  measured  (a^)  250^  = 0.89  of  the  selective  Cu  absorber  is 
in  good  agreement  with  the  measurements  of  Christie  (A6)  and  with 
Hottel  and  Unger  (72).  The  one  available  measurement  of  emittance, 
“65°C  ~ 0.20,  is  by  Christie.  Temperature  effects  are  not  known, 

although  some  of  Singham's  (44)  data  suggest  that  they  may  be  negli- 
gibly small. 

Tien  e^  al.  (73)  have  measured  the  t of  1/4-mil  Mylar;  it 

IR 

remains  quite  high  except  for  some  strong  resonances  in  the  8-10  pm 
region.  Therefore,  the  Mylar  absorbs  little  energy  and  probably 
operates  here  at  temperatures  only  slightly  above  ambient;  hence, 
its  measured  (t^) 250c  “ 0.85  should  apply.  For  the  reflector  and 

flux  bulb,  ambient  operation  is  assumed. 


Thermal  Properties  (Conductivity) 


The  required  thermal  properties  of  the  materials  used  either 
are  taken  from  the  literature  or  are  available  from  the  manufacturers. 
Although  the  conductivity  is  a function  of  temperature,  it  is  so  con- 
sidered only  for  the  insulation  and  for  the  argon  gas  (the  viscosity 
and  density  of  the  gas  are  also  taken  as  variable) ; the  density  of 
n-butyl  phthalate  in  the  flow  measuring  manometer  is  also  considered 
a temperature-dependent  property.  The  conductivities  of  graphite,  Cu, 
and  Pyrex  (and  the  Hg  density  in  the  pressure  manometer)  are  taken  as 
constants  here  because  resulting  errors  in  the  mathematical  analysis 
will  be  slight,  and  the  analysis  (see  Chapter  IV)  will  become  more 
tractable  thereby. 

As  in  the  case  of  optical  properties,  difficulties  were 
encountered  in  characterizing  the  graphite.  The  manufacturer  (45) 
was  unable  to  provide  date  on  the  thermal  conductivity,  but  did 
supply  information  on  the  resistivity  (which  is  known  to  vary 
inversely  with  the  thermal  conductivity) , and  the  coefficient  of 
thermal  expansion  (CTE) . 787G  graphite  is  anisotropic;  its  proper- 

ties vary  because  of  the  extrusion  process  used  in  its  manufacture. 

-4 

The  electrical  resistivity  is  7.1  x 10  ohm-cm  measured  parallel  to 
the  grain  (i.e.,  extrusion),  and  the  CTE  is  2.4  x 10  ^ °C  ^ (parallel) 
and  4.3  X 10  ^ °C  ^ (perpendicular  to  the  extrusion  axis),  yielding  a 
CTE  anisotropy  ratio  of  2.4/4. 3 = 0.56.  The  Kirk-Othmer  Encyclopedia 
(74)  gives  data  on  many  graphite  types:  an  AGOT  graphite  has  a resis- 
tivity (parallel  to  the  extrusion  axis)  of  7.35  x 10  ^ ohm-cm  and  a 
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GTE  ratio  of  0.57,  closely  matching  similar  properties  of  the  787G 
graphite.  Its  thermal  conductivity  is  173  W m ^ °C  Touloukian 
(53)  gives  data  on  the  conductivity  of  grade  "7087"  (i.e.,  assumed 
to  be  787G,  as  discussed  earlier)  Speer  graphite,  but  only  over  the 
temperature  range  800--1900°K.  Extrapolating  this  to  300°K  gives 
k = 210  W m ^ °C  decreasing  rapidly  to  136  W m ^ °C  ^ by  the 
first  data  point  at  800°K.  Because  of  the  uncertainty  in  the 
conductivity,  and  graphite  temperatures  above  300°K  in  the  present 
work,  a value  of  k = 173  W m ^ °C  ^ is  selected  as  typical. 

D.  Operation 

1 . Preliminary  Alignment  and  Calibration 

Before  experiments  could  be  performed  with  the  device,  the 
following  tasks  had  to  be  performed; 

(1)  Calibration  of  the  flow-rate-determining  nozzle.  This 
had  been  constructed  with  a diameter  ratio  (constriction  diameter/ 
upstream  diameter)  well  under  the  maximum  of  0.5  recommended  (76) 
for  accurate  metering  at  low  Re.  The  nozzle  was  also  made  suffi- 
ciently long  for  flow  perturbations  to  be  minimal  (76,77).  It  was 
calibrated  against  a wet  test  meter  (Precision  Scientific)  using 
the  actual  pump  of  the  experimental  device.  Various  flow  rates,  to 
a maximum  deliverable  108.5  liters/minute  of  air  at  ambient  con- 
ditions, were  used.  The  flow  was  then  corrected  to  provide  a 


J 


calibration  for  argon  flow  at  standard  temperature  and  pressure 
(STP,  25°C  and  760  torr  here)  (see  Table  7 and  Figure  20). 
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(2)  For  accurate  E-W  tracking  at  a constant  15°/h,  the 
receiver  axis  has  to  be  parallel  to  the  earth's  axis  of  rotation 
(i.e.,  N-S) . The  required  alignment  was  accomplished  as  follows: 
the  time  of  solar  noon  was  determined  from  a solar  ephemeris  tabu- 
lation (these  are  available  from  several  sources,  e.g.,  K and  E 
Company) . The  shadow  cast  by  a vertical  structure  is  at  that 
instant  in  the  true  N-S  plane.  The  shadow  position  was  marked 
(for  minimum  error,  this  procedure  was  repeated  a number  of  times, 
on  widely  spaced  days)  and  the  receiver  axis  was  aligned  along  the 
markings . 


(3)  Each  mirror  strip  of  the  Archimedes  concentrator  was 
focused  onto  the  receiver  aperture.  The  method  was  straightforward: 
all  the  strips  except  one  were  first  covered  with  opaque  material. 
The  device  was  oriented  N-S,  and  the  mirror  arr  .y  positioned  (motor- 
driven)  perpendicular  to  the  plane  defined  by  the  receiver  axis  and 
the  sun.  The  single  uncovered  strip  was  quickly  rotated  until  its 
reflected  beam  fell  onto  the  aperture,  meanwhile  tracking  the  sun 
with  the  step  motor.  An  accurate  tracking  speed  (i.e.,  pulse/ 
time  by  the  timing  circuit  to  the  stepping  motor,  ca.  2.5  pulse/s 
at  a motor  step  increment  of  0.9°/pulse)  was  then  determined  by 
trial-and-error  over  long  time  periods.  No  provision  was  made  for 
feedback  of  tracking  error  to  the  timing  circuit,  so  this  had  to  be 
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Table  7.  Nozzle  calibration. 


4P,  cm  H^O® 

b 

Flow  rate  , air, 
liter  min  ^ 

Equivalent  4P 

■>'  ‘=16"22‘>4' 

Flow  rate*^,  Ar 
liter  min 

14.5 

108.45 

13.9 

97.31 

11.6 

96.56 

11.1 

86.64 

9.0 

86.65 

8.6 

77.73 

6.8 

76.17 

6.5 

68.32 

5.0 

67.11 

4.0 

60.20 

4.7 

65.70 

4.5 

58.93 

3.8 

60.31 

3.7 

54.10 

3.3 

56.35 

3.2 

50.55 

2.5 

49.84 

2.4 

44.71 

1.6 

41.34 

1.5 

37.08 

0.6 

32.56 

0.6 

29.21 

0.4 

20.10 

0.4 

18.03 

Measured,  see  Figure  20 

'^Measured  average  over  5 min.,  at  ambient  conditions  (T  «•  24.5  *C, 

P » 729  torr,  P water  vapor  in  wet  test  meter  * 23.2  torr,  rel. 
humidity  * 87%) 

c 

n-butyl  phthalate,  cm,  at  25®C 

*^i.e.,  at  760  torr,  25®C,  determined  by:  (flow  rate,  air) 

[(729  + 23. 2)/760] (0.905) (298/297.5) , where  Reuther  (75)  recommends 
a factor  of  0.905  to  correct  for  air/argon  differences.  Ideal  gas 
behavior  is  assumed. 


done  quite  carefully.  After  what  was  considered  acceptable  tracking 
accuracy  had  been  achieved  (i.e.,  no  manual  corrections  were  neces- 
sary for  a period  of  2-3  h) , the  remaining  mirror  strips  were 
aligned  singly  by  uncovering  each  separately,  focusing  its  beam, 
tightening  the  locknuts  fixing  it  to  the  array  frame,  and  then 
covering  it  again  while  continuously  tracking  Che  sun.  As  an 
additional  precaution,  the  mirror  strip  which  had  been  aligned  first 
was  kept  uncovered  throughout,  and  its  focus  monitored.  The  focusing 
was  performed  several  times,  at  different  array  positions  on  the  track, 
as  the  optimum  alignment  for  a particular  position  was  not  the  best 
for  all  positions;  a compromise  focus  was  the  result,  with  an  esti- 
mated 5%  (i.e.,  2 mirror  strlps/40)  loss  of  the  flux  from  misalignment 
(see  Figure  13) . It  is  assumed  that  all  of  the  flux  missing  the 
window  is  lost,  even  though  it  impinges  on  the  aperture  sides. 

Once  the  mirrors  were  properly  aligned,  no  further  corrections  were 
required  thereafter  as  long  as  the  receiver  position  did  not  alter, 
as  after  changing  absorbers  - which  minutely  displaced  the  receiver 
axis. 


(4)  Taking  advantage  of  the  selective  radiative  properties  of 
2 

copper,  a thin  10  cm  plate  was  mounted  within  an  evacuated  Pyrex 
glass  bulb.  When  this  device  is  exposed  to  sunlight,  the  plate  warms. 
The  only  important  thermal  loss  occurs  by  conduction  along  a copper 
support  rod,  hence  the  incident  direct  solar  flux  can  be  calculated 
from  the  difference  between  the  measured  plate  and  ambient  tempera- 
tures. The  device  was  mounted  on  the  tracking  mirror  array,  and  was 
masked  so  the  plate  only  "saw"  a small  angle  of  the  sky  (see  Figure  15) . 
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2.  Experimental  Procedure 

The  concentrator,  which  focuses  only  the  direct  component  of 
the  solar  flux  onto  the  receiver  aperture,  limits  operation  to  cloud- 
less days.  Tracking  was  initiated  ca.  3 hours  prior  to  solar  noon. 

The  receiver  required  approximately  90  minutes  to  attain  its  asymp- 
totic temperature.  The  first  set  of  readings  was  taken  at  the  end 
of  this  time  period.  Flow  was  started  thereafter  with  attainment 
of  steady-state  conditions  within  30-60  minutes  (depending  on  the  flow 
rate  selected).  Appropriate  temperature,  pressure,  and  solar  flux 
measurements  were  then  made  every  30  minutes  until  the  end  of  the 
experiment.  In  greater  detail: 

(1)  The  collector  was  wheeled  out  of  its  shed  and  its  axis 
aligned  N-S.  The  entire  device  was  light  enough  so  that  this  was 
easily  performed  by  one  person. 

(2)  The  receiver  was  filled  with  argon  (or  still  contained 
argon  from  a recent  past  experiment) . 

(3)  The  reflector  was  cleaned  with  a commercial  cleaner 
(Wlndex)  and  acetone,  in  that  order. 

(4)  The  mirror  array  was  positioned  on  the  tracks  such  chat 
the  reflected  and  concentrated  solar  beam  was  focused  onto  the 
aperture.  A fast-forward  provision  in  the  timing  circuit  allowed 
this  to  be  done  quite  rapidly;  Che  stepping  motor,  although  small, 
had  sufficient  torque  (ca.  29  kp  cm  max.  in  its  0.9*’/step  mode) 
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for  this  maneuver.  The  track  diameter  and  array  width  combine  to 
restrict  permissible  array  travel  to  ca.  45®  off-horlzontal,  or 
to  ca.  6 hours  total  operation  when  centered  about  solar  noon. 

The  normal  tracking  speed  (pulse  rate)  was  set,  and  tracking  began. 
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(5)  After  attainment  of  steady-state  stagnation  conditions 
(Figure  21),  temperature  readings  were  taken  and  recorded.  This 
was  done  to  check  thermocouple  response  (the  switches  could  affect 
this.  If  not  cleaned  frequently).  Also,  the  asymptotic  temperatures 
gave  some  Indication  of  the  solar  flux  level;  If  a wide  discrepancy 
between  the  two  solar  flux  measurements  was  observed,  one  could 
have  an  Indication  as  to  which  was  the  more  accurate  measurement 

by  comparing  these  measurements  with  previous  experimental  flux 
and  asymptotic  temperature  measurements. 

(6)  The  system  was  re-charged  with  Ar  for  ca.  2 minutes  at 
a gauge  pressure  of  ca.  100  torr.  The  pump  was  run  at  a low  flow 
rate  (AP  * 2 cm  phthalate)  for  ca.  15  s during  this  Interval. 

(7)  The  pump  was  started  and  the  flow  rate  set.  A trickle 
of  Ar  was  fed  Into  the  envelope  around  It,  keeping  that  gently 
Inflated. 

(8)  If  used,  pre-heat  was  set  and  maintained  at  some  chosen 
level  on  the  Varlac.  T5  and  T6  monitor  the  preheat  and  attainment 
of  steady-state  conditions  with  It  operational. 
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Figure  21.  Time-dependent  temperature  response  (measured) 
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(9)  Readings  of  temperatures,  pressures,  and  other  parameters 
were  taken  every  30  minutes  thereafter  (Table  8) . Although  steady- 
state  conditions  were  achieved  within,  at  most,  an  hour  after 
initiating  flow,  the  flow  rate  was  generally  maintained  the  same 
throughout  the  day  as  a check  on  reading  accuracy,  etc.  Occasionally, 
two  or  even  three  different  flow  rates  have  been  employed  on  the  same 
day  to  take  advantage  of  the  infrequent  occurrences  at  State  College, 
PA,  of  bright,  clear  skies. 

E.  Experimental  Results 

1.  Flow  Range  Restrictions 

One  complication  in  studying  a small  thermodynamic  model  is 
that  various  effects,  such  as  end  losses,  transverse  temperature 
gradients,  etc.,  which  may  play  a negligible  role  in  a full-scale 
device,  often  exert  a strong  influence  on  the  experimental  results. 

In  the  present  study,  it  is  found  that  at  flows  much  greater  than 
40  liters  min  ^ (ca.  Re  » 800,  or  ca.  40  cm  s ^ mean  velocity  in 
the  receiver,  STP)  flow  disturbances  near  the  receiver  outlet  cause 
anomalous  temperature  profiles  in  the  downstream  half  of  the  receiver. 
The  effect  worsens  with  increasing  flow  rate,  although  examination  of 
the  data  (see  Table  8)  reveals  that  the  heat  transfer  rate  is  not 
affected.  The  disturbances  are  most  likely  caused  by  the  too-short 
exit  end  piece,  and  the  abrupt  constriction  at  the  small-diameter 
copper  return  line.  This  causal  hypothesis  can  be  readily  tested 
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Table  8.  Experimental  results  (flowing  gas). 


* 

Date 

19. IX. 

76 

23. IX 

.76 

24.  IX 

.76 

29. IX. 76^ 

Time 

1:45 

2:30 

2:00 

3:00 

1:45 

2:45 

3:30 

Angle,  cm 

78.5 

91.5 

89 

114.5 

84 

102 

131 

AP  , , .cm 

2.8 

11.7 

1.9 

10.9 

0.8 

5.4 

3.1 

pntn. 

A^Hg* 

-1.6 

5.1 

-2.9 

5.1 

-1.6 

1.2 

-0.7 

P , cm 

72.3 

72.3 

72.5 

72.5 

72.8 

72.7 

72.3 

T , “C 

28 

27 

22.5 

23 

21 

20 

20 

-1 

Flux,  ly  min 

1.16 

1.06 

1.06 

0.96 

1.12 

1.02 

0.81 

Shadow,  cm 

0 

0 

1.5 

1.5 

1.5 

1.5 

4 

-1 

Wind,  m s 

0 

0 

4 

4 

1 

1 

2 

Tl,  “C 

— 

— 

44.5 

45.5 

49 

47 

40.5 

T3 

47 

42 

36 

37 

50 

43 

37 

T4 

93 

68 

70 

60 

96 

77 

153 

T5 

55 

43 

45 

37 

68 

43 

229 

T6 

102 

74 

92 

66 

127 

84 

226 

T14 

163 

83 

170 

87 

235 

127 

257 

T15 

233 

172 

231 

173 

296 

232 

300 

T7 

220 

168 

210 

175 

283 

231 

245 

T19 

128 

113 

122 

118 

208 

147 

135 

T16 

238 

130 

233 

138 

319 

215 

278 

T17 

270 

150 

271 

158 

340 

246 

290 

T18 

298 

228 

300 

228 

350 

295 

316 

T8 

60 

55 

46 

46 

63 

66 

48 

T9 

216 

172 

211 

177 

264 

237 

224 

TIO 

270 

180 

267 

189 

315 

272 

286 

Til 

272 

228 

273 

238 

309 

292 

286 

T12 

45 

53 

35 

46 

31 

46 

31 

(T16)g®,  “C 

331 

342 

372 

352 

i. 
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Table  8.  Continued. 


* 

Date 

3.x. 

76® 

4.x. 

76® 

11. X. 

76 

12.x. 

76 

Time 

1:15 

2:45 

1:50 

3:00 

2:15 

3:15 

1:15 

2:30 

Angle,  cm 

91.5 

114.5 

99 

119.5 

114.5 

140 

101.5 

120.5 

‘‘■phth.-  “ 

11.4 

8.5 

1.2 

6.1 

1.5 

3.1 

12.3 

5.9 

/^Pjjg.  cm 

5.6 

3.8 

-3.1 

2.2 

3.3 

4.9 

15.0 

8.5 

1*00  > 

72.4 

72.4 

72.8 

72.8 

73.4 

73.4 

73.2 

73.1 

T , ®C 

23 

24 

23 

22 

14 

15 

16 

16.5 

-1 

Flux,  ly  min 

1.02 

0.92 

0.97 

0.85 

0.95 

0.79 

i.OO 

0.90 

Shadow,  cm 

5 

5 

5 

5 

8.5 

8.5 

9 

9 

-1 

Wind,  m s 

2.5 

3 

2 

2.5 

1.5 

1.5 

3 

3 

Tl,  "C 

47 

45 

47 

45 

42 

40 

42.5 

42 

T3 

43 

49 

57 

43 

29 

29 

26 

28 

T4 

131 

136 

162 

152 

35 

40 

36 

40 

T5 

166 

182 

259 

197 

28 

27 

35 

31 

T6 

176 

187 

240 

206 

65 

55 

45 

51 

T14 

196 

208 

270 

224 

95 

71 

60 

64 

T15 

257 

260 

296 

276 

206 

145 

145 

158 

T7 

224 

210 

259 

239 

216 

184 

160 

170 

T19 

145 

135 

172 

123 

115 

115 

106 

85 

T16 

224 

235 

299 

251 

219 

168 

111 

138 

T17 

228 

240 

313 

260 

271 

208 

125 

164 

T18 

283 

282 

326 

294 

300 

250 

169 

206 

T8 

58 

62 

74 

56 

59 

56 

48 

45 

T9 

226 

207 

250 

223 

251 

227 

193 

201 

TIO 

246 

245 

302 

271 

285 

237 

141 

192 

Til 

278 

265 

292 

282 

298 

289 

226 

248 

T12 

48 

51 

36 

45 

32 

37 

54 

39 

(T16)g®,  “C 

- 

331 

361 

370 

91 


Table  8.  Continued. 


Date*  15.x. 76'^  23.X.76  14. XI. 76  19.XI.76‘^ 


Time 

1:30 

2:45 

1:30 

2:30 

1:00 

12:30 

1:51 

Angle,  cm 

105 

131 

118 

140 

164 

161 

203 

^Vh.’ 

3.0 

12.1 

1.4 

6.0 

12.7 

0.25 

0.9 

5.3 

14.8 

3.2 

8.6 

15.4 

0.9 

2.2 

P».  cm 

72.1 

71.9 

73.1 

73.1 

72.5 

71.9 

71.' 

Too'  “C 

24 

25 

12 

13 

6 

14 

14 

-1 

Flux,  ly  min 

0.91 

0.79 

0.90 

0.80 

0.74 

0.67 

0.5: 

Shadow,  cm 

10 

10 

12.5 

12.5 

20.5 

21.5 

21.: 

Wind,  m s ^ 

8 

8 

3 

3.5 

5.5 

8.5 

8.5 

Tl,  “C 

46 

45.5. 

34 

33.5 

28 

32 

30 

T3 

31 

30 

19 

22 

9 

17 

16 

T4 

39 

34 

27 

25 

14 

25 

19 

T5 

31 

34 

19 

23 

17 

22 

17 

T6 

47 

43 

36 

32 

19 

33 

26 

T14 

67 

52 

56 

42 

26 

71 

35 

T15 

157 

119 

139 

115 

77 

128 

90 

T7 

180 

126 

182 

125 

102 

205 

127 

T19 

173 

145 

170 

140 

123 

174 

130 

T16 

157 

102 

172 

103 

66 

212 

133 

T17 

202 

107 

228 

126 

82 

272 

193 

T18 

240 

147 

263 

166 

122 

292 

223 

T8 

49 

38 

43 

40 

22 

32 

30 

T9 

223 

154 

226 

170 

143 

206 

164 

TIO 

227 

126 

251 

149 

98 

267 

219 

Til 

265 

200 

275 

222 

193 

259 

231 

T12 

36 

45 

23 

37 

28 

15 

18 

(T16)  ®,  "C 
s 

353 

333 

348 

ca.  300 

1 
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Table  8.  Continued. 


Date 

14. XII. 76 

11.11.77 

2. III. 77 

9. III. 77 

11. III. 77 

Time 

12.:45 

1:15 

1:30 

1:15 

1:00 

Angle,  cm 

193 

135 

109 

96.5 

91.5 

^^hth.’ 

2.0 

12.6 

5.9 

6.0 

4.4 

3.9 

15.9 

9.0 

9.0 

7.4 

Pc.> 

73.1 

72.9 

73.4 

73.2 

73.3 

T^,  cm 

1 

14 

5 

20 

23 

Flux,  ly  min  ^ 

0.66 

0.84 

1.03 

1.04 

1.01 

Shadow,  cm 

25.5 

15 

9 

5 

4.5 

Wind,  m s ^ 

2.5 

6 

6.5 

5 

2.5 

Tl,  “C  ’ 

22 

35 

29 

45 

45 

T3 

7 

19 

15 

28 

36 

T4 

9 

21 

18 

40 

47 

T5 

5 

24 

13 

29 

33 

T6 

10 

28 

28 

52 

61 

T14 

17 

35 

42 

63 

76 

T15 

72 

101 

102 

155 

174 

T7 

136 

111 

141 

160 

173 

T19 

138 

133 

156 

150 

184 

T16 

109 

73 

116 

134 

158 

T17 

175 

98 

152 

161 

187 

T18 

207 

135 

133 

214 

231 

T8 

25 

30 

34 

41 

57 

T9 

200 

146 

185 

195 

215 

TIO 

210 

111 

176 

186 

214 

Til 

249 

189 

239 

242 

261 

T12 

10 

37 

23 

36 

46 

(T16)  °C 

— 

340 

379 

367 

366 

L 
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i 

I 


1^ 

Table  8.  Continued 


f 

K 

A 

Date 

1.IV.77 

12. IV. 77^ 

15. IV. 77 

17. IV. 77^ 

8.V.77^ 

t 

\ 

Time 

12:45 

1:45 

2:15 

1:45 

1:35 

^ngle , cm 

67.5 

71 

81 

68 

42 

> 

3.0 

1.55 

2.5 

1,45 

1.5 

i 

> 

5.8 

3.8 

4.4 

3.0 

3.3 

P , cm 

73.3 

73.5 

73.0 

73.1 

72.3 

T , cm 

14 

30.5 

25 

26 

22 

f 

♦ 

% 

Flux,  ly  min  ^ 

1.22 

1.09 

1.15 

1.16 

1.38 

Shadow,  cm 

0(4) 

0(7.5) 

0(9) 

0(10) 

0(16) 

t 

Wind,  m s ^ 

1.5 

7.5 

4 

2 

7.5 

Tl,  “C 

39 

50.5 

50 

49 

43 

T3 

30 

54 

45 

64 

28 

T4 

61 

152 

75 

226 

95 

T5 

30 

206 

42 

333 

99 

T6 

82 

219 

99 

321 

149 

1 

T14 

94 

240 

114 

343 

169 

; 4 

T15 

196 

287 

213 

364 

250 

T7 

191 

222 

190 

256 

167 

T19 

200 

213 

189 

248 

195 

■ • 

T16 

187 

272 

206 

356 

238 

T17 

228 

288 

241 

355 

265 

T18 

271 

315 

283 

372 

287 

T8 

48 

59 

59 

65 

34 

T9 

226 

210 

196 

228 

177 

1 

TIO 

253 

281 

250 

321 

248 

Til 

268 

256 

242 

278 

215 

T12 

36 

46 

39 

43 

30 

(T16)g®,  “C 

388 

364 

390 

378 

358 

L 
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Table  8.  Continued. 

Date*  14.V.77^  26.V.77‘^ 


Time 

1:45 

2:45 

1:30 

2:30 

3:30 

Angle,  cm 

39.5 

53.5 

33 

43 

66 

‘'’phth.-  ™ 

2.9 

2.85 

0.65 

1.0 

1.5 

SPhj.  •=• 

5.3 

5.2 

2.0 

2.8 

3.3 

72.6 

72.6 

72.9 

72.9 

72.9 

^00  ’ 

30 

33 

33 

34 

32 

Flux,  ly  min 

1.36 

1.27 

1.34 

1.28 

1.19 

Shadow,  cm 

0(18) 

0(18) 

0(20) 

0(20) 

0(20) 

Wind,  m s ^ 

4 

4.5 

4 

4 

4 

Tl,  “C 

53 

55 

53 

56 

53 

T3 

49 

58 

49 

51 

49 

T4 

105 

149 

80 

77 

74 

T5 

95 

182 

58 

51 

47 

T6 

150 

209 

141 

133 

122 

T14 

152 

221 

173 

164 

151 

T15 

250 

287 

264 

248 

230 

T7 

216 

238 

242 

232 

216 

T19 

214 

225 

218 

210 

188 

T16 

235 

263 

278 

256 

220 

T17 

245 

272 

305 

288 

262 

T18 

276 

297 

326 

315 

289 

T8 

51 

56 

56 

59 

56 

T9 

183 

189 

189 

182 

174 

TIO 

241 

262 

246 

253 

244 

Til 

218 

231 

2:8 

226 

221 

T12 

46 

46 

38 

43 

42 

(T16)  "C 

s 

375 

384 

Table  8.  Continued. 


* 

Date 

3. VI. 

4. VI. 

,77^ 

Time 

1:30 

2:30 

3:30 

1:30 

2:30 

3:30 

Angle,  cm 

30.5 

42 

63.5 

30.5 

42 

63.5 

2.0 

2.5 

3.0 

1.0 

1.5 

2.0 

4.1 

4.7 

5.4 

2.6 

3.3 

4.0 

73.2 

73.2 

73.2 

73.2 

73.2 

73.2 

T^,  cm 

25 

26 

27 

33 

34 

36 

Flux,  ly  min 

1.41 

1.36 

1.22 

1.36 

1.31 

1.19 

Shadow,  cm 

0(23) 

0(23) 

0(23) 

0(23) 

0(23) 

0(23) 

i! 

Wind,  m s ^ 

4 

4 

4 

3.5 

3.5 

3.5 

'V 

Tl,  "C 

48 

49 

50 

55 

56 

57 

i«  * 

T3 

39 

40 

42 

47 

44 

44 

■- 

T4 

66 

61 

64 

74 

72 

70 

T5 

38 

36 

39 

52 

46 

45 

T6 

107 

103 

98 

129 

117 

109 

. 

T14 

127 

113 

100 

165 

144 

126 

■’ 

T15 

219 

206 

193 

234 

217 

200 

1 

T7 

227 

203 

193 

247 

224 

212 

■ * 

T19 

198 

187 

173 

222 

207 

193 

T16 

210 

194 

181 

256 

228 

205 

T17 

247 

230 

212 

284 

260 

237 

T18 

274 

258 

244 

306 

289 

261 

/i 

T8 

45 

46 

48 

55 

54 

54 

t 

T9 

172 

160 

153 

181 

170 

162 

> 

TIO 

233 

221 

209 

242 

237 

224 

Til 

213 

207 

198 

216 

215 

207 

T12 

38 

40 

41 

40 

43 

44 

(T16)  “C 

359^ 

333^ 

I 


The  terms  and  symbols  in  this  column  of  Table  8 have  the  following 
significance: 


date  The  date  of  the  experiment. 

time  The  time  the  data  was  recorded. 

angle  Horizontal  length  of  shadow  cast  by  a vertical  91.5  cm 
high  rod. 

AP^hth  Pressure  drop  across  the  flow  measuring  constriction.  The 
manometer  fluid  was  n-butyl  phthalate. 

AP„  Fluid  static  pressure  at  the  flow  measuring  constriction. 

The  manometer  fluid  was  mercury;  one  side  of  the  manometer 
was  open  to  the  atmosphere. 

P^  The  ambient  atmospheric  pressure,  in  centimeters  of  mercury. 

T^  The  ambient  temperature. 

flux  The  total  solar  flux  on  a horizontal  surface  in  langley 

min”l  (langley  • gram  calorie  cm“2) , Measured  by  The 
Pennsylvania  State  University  Meteorological  Station. 

shadow  Length  of  aperture  shadow,  at  the  receiver  inlet  end. 


wind  Wind  velocity.  Measured  by  The  Pennsylvania  State  University 

Meteorological  Station. 

T1  - T18  Temperatures  measured  by  various  thermocouples  (see  Figure 
19  and  Table  5) . 


Preheat  used 

^Mylar  cover  installed  over  aperture;  all  subsequent  experiments  had 
this  feature 

New  digital  pyrometer  used  here  and  thereafter 

Selective  absorber  used  here  instead  of  the  graphite 

No  shadow  at  inlet  hereafter;  in  parentheses  is  the  length  of  aperture 
shadow,  now  at  the  outlet  end 

^Low  reading;  probable  cause  is  absorber  buckling  in  the  vicinity  of 
the  thermocouple,  bringing  it  close  to  the  relatively  colder  glass  tube 

a 

Steady-state  stagnation  temperature,  T16,  measured  just  prior  to  the 
start  of  flow 
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by  installation  of  a longer  exit  end  piece,  and/or  a larger- 
diameter  return  line.  However,  the  size  of  the  storage  shed, 
which  was  barely  large  enough  for  the  collector  as  is,  prevented 
such  a test. 

This  perturbed  condition  limits  the  upper  flow  range  avail- 
able for  experiment  and  analysis.  No  attempt  has  been  made  to 
formulate  the  perturbed  behavior  in  mathematical  terms,  because 
it  appears  to  be  too  complex  to  treat.  Very  low  flow  rates  are 
also  ruled  out,  since  below  Re  ■ Pr  » 300  axial  conduction  in  the 
gas  may  no  longer  be  negligible  (78) ; it  is  certainly  a factor  at 
Re  • Pr  ^ 100  (79). 

2.  Temperature  Reading  Accuracy 

Temperatures  were  at  first  measured  with  an  old  Instrument, 
a Hoskins  pyrometer.  These  readings  were  found  to  be  inaccurate 
because  the  instrument: 

(1)  Was  not  temperature-compensated.  A correction  was 
made  by  attaching  a thermocouple  to  a calibrated  thermometer  in 
the  ambient  environment.  This  became,  in  effect,  the  reference 
temperature. 

(2)  Above  100®C  readings  were  higher  than  actual,  the 
difference  (even  as  a percentage)  Increasing  with  higher  tempera- 
tures. The  error  was  considerable:  for  example,  at  511°C  the 
Instrument  read 
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(3)  The  instrument  had  a low  internal  resistivity.  The 
electrical  resistivity  of  28  ga.  Chromel-Alumel  wires  (the  thick- 
ness used  here)  is  about  12.5  ohm/double  m (i.e.,  for  both  wires; 
the  value  is  a sum  of  individual  resistivities) . Wire  lengths  of 
up  to  5 m were  used  and  the  readings  from  the  longer  wires  were 
thus  significantly  lowered. 

A digital  readout  instrument  has  therefore  been  used  for  most 
measurements.  It  was  calibrated  at  0®C,  lOO^C,  and  511®C;  its 
internal  cold  junction  reference  was  also  calibrated,  as  prescribed. 
The  high  internal  impedance  of  this  instrument  permits  accurate 
measurement  even  with  the  longest  wires  used  here.  As  a check  on 
accuracy,  a voltage  corresponding  to  689®C  was  fed  into  the  longest 
lead  (which  went  through  two  switches) . The  instrument  error  then 
was  on  the  order  of  1®C.  When  all  the  experiments  had  been  completed, 
this  check  was  repeated,  with  similar  results.  Also,  an  ice  bath 
was  employed  on  several  occasions  (as  well  as  the  previously- 
mentioned  thermocouple/thermometer  combination  for  ambient  tempera- 
ture comparison,  which  was  used  throughout)  to  check  zero  drift  and 
internal  compensation  accuracy.  The  readings  were  generally  within 
1°C  of  actual,  with  an  occasional  maximum  error  of  3®C.  Considering 
the  prevailing  operating  conditions  - outside  environment,  with 
varying  sun,  wind,  and  temperature  - this  was  very  good  performance. 
Repeatability  was  excellent,  always  within  1®C. 


The  old  Instrument  was  calibrated  against  the  new  one, 
thermocouple-by-thermocouple,  and  at  various  temperatures.  The 
data  acquired  by  means  of  the  old  Instrument  were  then  corrected, 
and  are  so  shown  here  (Table  8) . 

3.  Material  Degradation 

At  the  end  of  the  experiments,  all  materials  whose  optical 
properties  had  been  previously  measured  (Table  6)  were  measured 
again  to  determine  the  effect  of  operation  at  high  fluxes  and/or 
high  temperatures.  After  ca.  200  hours  of  operation,  the  graphite 
absorber  exhibited  no  measurable  change  In  Its  properties.  The 
total  reflectance  of  the  mirrors,  as  well  as  the  total  transmit- 
tance of  the  Mylar  cover  also  remained  unchanged.  The  diffuse 
properties  of  the  mirror  and  the  Mylar  did  show  some  slight  degrada- 
tion: the  diffuse  reflectance  of  the  mirror  Increased  by  ca.  3 
percentage  points,  and  the  diffuse  component  of  transmittance  In 
the  Mylar  Increased  by  about  4 percentage  points.  These  values 
are  over  the  wavelength  range  400  < X < 1600  nm,  which  contains 
ca.  90%  of  the  solar  energy.  The  Mylar  cover  was  replaced  several 
times  (e.g.,  to  clean  the  Pyrex  window),  and  this  particular  sample 
had  been  subjected  to  a concentrated  flux  for  ca.  72  hours. 

Because  degradation  was  so  minuscule,  and  Its  effect  on 
performance  difficult  to  quantify  (also,  the  degradation  rate  may 
not  have  been  constant) , It  was  not  Incorporated  Into  the  calcula- 
tions of  Chapter  IV.  The  Pyrex  was  not  remeasured  since  no  acceptable 
method  had  been  found  to  expose  a small  sample  to  the  environment 
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"seen"  by  the  glass  receiver  pipe,  while  the  copper  absorber  surface 
was  so  badly  degraded  (after  only  ca.  20  h of  operation)  that  measure- 
ments would  have  been  pointless. 

F.  Conclusions 

The  experimental  device  achieves  high  temperatures  and  performs 
essentially  as  expected.  Stagnation  temperatures  show  that  the 
32-to-123  cm  axial  section  of  the  aperture  exhibits  a quite  flat 
temperature  profile  near  the  equinox  (Figure  22) . Off-equinox 
operation,  when  significant  shadowing  occurs  over  one  aperture  end 
section,  indicates  some  end-effect  penetration  into  the  32-123  cm 
section  (Figure  23).  The  influence  is  not  great,  and  would  probably 
be  even  smaller  under  gas  flow  conditions.  Its  effect  on  the  axial 
gas  bulk  temperature  profile  is  considered  negligible;  i.e.,  the 
32-cm  gas  bulk  temperature  (the  starting  temperature  for  calcula- 
tions) already  incorporates  any  such  effect,  if  present,  while  the 
123-cm  thermocouple  position,  being  upstream  of  any  shadow  at  the 
outlet,  would  be  negligibly  affected  thereby.  Examination  of  the 
experimental  results  (Table  8)  reveals  large  variations  in  tempera- 
tures with  flow  rate;  ambient  temperature,  aperture  shadowing  (at 
the  inlet),  and  the  solar  flux  also  exert  an  influence,  A full 
discussion  of  these  is  reserved  for  Chapter  IV. 


AXIAL  POSITION,  cm 


I Figure  22.  Steady-state  stagnation  temperatures  of  receiver 

* measured  on  23. III. 77,  12:00  p.m.  The  direct 

' solar  flux  was  ca.  870  W m-2^  and  the  ambient 

temperature  was  7'’C. 
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AXIAL  POSITION,  cm. 

Figure  23.  Steady-state  stagnation  temperatures  of  receiver 
measured  on  14.V.77,  12:15  p.m.  The  direct  solar 
flux  was  ca.  800  W m"^,  and  the  ambient  temperature 
was  25'’C. 


CHAPTER  IV 


MATHEMATICAL  MODEL  (FLOWING  GAS) 


A.  Introduction 

One  would  like  to  be  able  to  accurately  predict  the  performance 
of  this  collector  design  In  all  possible  configurations.  Prior  to 
such  an  attempt,  a mathematical  model  of  the  device  must  be  devel- 
oped that  adequately  accounts  for  Its  observed  behavior.  In 
practice,  the  mathematical  model  was  formulated  In  successive 
stages  of  complexity  and  refinement.  William  of  Ockham's  principle 
(l.e.,  the  simplest  model  to  account  for  observed  phenomena  is  the 
most  likely  to  represent  the  true  situation)  was  observed. 

As  in  Chapter  II,  the  calculation  method  for  asymptotic 
conditions  is  straightforward.  This  temperature  occurs  when  the 
incoming  solar  flux  is  just  balanced  by  the  various  thermal  losses 
(l.e.,  there  is  no  net  heat  gain  by  the  gas  which  implies  either  gas 
exiting  from  an  infinitely  long  collector  or  a static  gas).  Two 
principal  loss  terms  were  considered:  radial  heat  conduction 
through  the  insulation  surrounding  the  glass  receiver,  and  re- 
radiation  by  the  aperture,  both  by  the  glass  and  the  asbestos  sides. 

A more  detailed  mathematical  formulation  is  required  for 
handling  gas  flows  at  various  rates.  Internal  heat  transfer  modes 
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and  their  magnitudes  must  be  considered.  The  laminar  flow  regime, 
tilt,  and  short  length  of  the  model  create  a situation  far  more 
intricate  than  that  within  the  long  full-scale  device  with  its 
highly  turbulent  flow  and  small  temperature  gradients. 

B.  Stagnation  (V  = 0)  Calculations 

A number  of  subsidiary  calculations  and  operations  must  be 
performed  before  formulating  an  expression  for  heat  transfer 
within  the  device.  These  are  discussed  below. 

1.  Solar  Flux  Computation 

The  solar  flux  available  to  the  device  is  determined  as 
follows : 

a.  Flux  bulb 

The  temperature  difference  across  this  instrument  (described 
in  Chapter  III)  is  proportional  to  the  direct  solar  flux  at  the 
reflector.  Since  the  bulb  is  masked  so  that  the  absorber  plate 
within  it  only  "sees"  the  sun,  and  is  at  the  same  angle  as  the 
reflector,  no  corrections  for  the  diffuse  flux  and  tilt  are  neces- 
sary. Thus,  using  the  one-dimensional  heat  conduction  equation  for 
the  copper  support  rod  (of  the  exposed  plate) : 


(A^  k AT)/[A2(ag)25< 


C ^1^'^S^25' 


where  the  first  area  is  the  cross-sectional  area  of  the  support 

rod,  and  the  second  is  the  area  of  the  plate;  the  absorptance  is 

for  the  plate,  while  the  transmittance  is  for  the  Pyrex  bulb; 

AT  = T1  - T . 

00 

Substituting  in  the  appropriate  values  (see  Table  4), 

q = (Tr/4)  (0.208  cm) ^ (lO^cm^m"^)  (3 . 89  W cm"^  '’c"^)(AT)/ 

(10  cm^) (0.39) (10.3  cm) (0.91) 

=36.4  (AT)  W m“^  . 

The  selective  radiative  properties  of  copper  have  been  used  to 

advantage  here.  The  very  low  near-ambient  plate  tempera- 

JLK 

tures  permit  neglecting  the  radiative  loss  from  the  plate, 
b . Eppley  pyranometer 

Since  it  is  essential  to  define  the  flux  as  accurately  as 
possible,  an  alternative  approach  is  taken.  The  Pennsylvania  State 
University  operates  a U.S.  Weather  Bureau  Station  and  measures  the 
total  steradian)  horizontal  flux  with  an  Eppley  pyranometer. 

The  accuracy  of  the  particular  instrument  used  is  currently  estimated 
to  be  of  the  order  of  95%  (80) , which  is  fairly  typical  (51) . With 
the  horizontal  flux  reading  one  can  calculate  the  perpendicular 
flux  from  the  solar  altitude  (defined  here  as  the  angle  of  the  sun 
above  the  horizon,  i.e.,  90®  minus  the  zenith  angle),  and  then  the 
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direct  component  from  an  empirical  correlation.  A further  overall 
correction  for  non-perpendicularity  of  the  array  to  the  solar  flux 
is  also  made. 

More  specifically: 


perpendicular  total  flux  = (horizontal  total  flux) /sin  9 

where  the  angle  9 is  the  pertinent  solar  elevation  in  the  vertical. 
Here,  a vertical  rod  and  its  shadow  are  used  to  find  9,  i.e.: 

9 =»  tan  ^(91.5  cm,  vertical  rod  ht.)/ 

(measured  rod  horizontal  shadow  length,  cm)  . 

The  flux  available  at  the  mirror  then  becomes: 

available  direct  flux  = 

(perpendicular  total  flux)  • (1  - diffuse  component)  • 

(angular  correction  for  non-perpendicular  collector)  . 

The  flux  of  diffuse  sky  radiation  reaching  the  ground  may  be 
computed  by  solving  the  equation  of  radiative  transfer  for  a turbid 
atmosphere.  However,  the  problem  is  a demanding  one  (81)  and  theoret- 
ical predictions  of  the  diffuse  component  are  inadequate  at  this  time 
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(60,82).  The  relationship  between  diffuse  and  total  radiation  is 
best  determined  from  empirical  correlations.  Liu  and  Jordan  (82) 
developed  such  a correlation  for  Blue  Hill,  MA.  This  correlation 
has  been  widely  used.  Since  it  is  based  on  data  from  only  a single 
site,  and  employs  the  older  value  of  the  extraterrestrial  solar 
constant  (2.00  ly  min  which  is  ca.  3%  higher  than  the  more 
accurate  current  measurement  (51)),  it  results  in  an  underprediction 
of  the  diffuse  component,  when  compared  with  other  sites  with 
similar  climatological/geographical  features  (83).  Some  error  also 
results  (82)  from  the  use  of  a monthly  average  to  account  for  daily 
variation  in  the  extraterrestrial  solar  constant. 

More  recently,  Ruth  and  Chant  (83)  have  suggested  forming  a 
correlation  based  on  local  conditions;  the  requisite  data  are  not 
available  for  State  College,  PA.  Their  standard  correlation, 
developed  from  four  Canadian  sites,  can  be  used  with  only  small 
error  at  other  locations.  Additional  slight  error  is  introduced 
when  applying  the  correlation,  which  is  based  on  daily  fluxes,  to 
hourly  fluxes  near  noon  here.  But  Orgill  and  Hollands  (84)  have 
shown,  in  their  study  of  hourly  dif fuse/total  solar  flux  radiation 
for  Toronto,  that  the  difference  between  daily  (83)  and  hourly 
ratios  is  insignificant. 

Tuller  (85)  found  a linear  relationship  adequate  for  a similar 
Canadian  data  base,  but  only  over  a restricted  range  of  values.  To 
cover  the  full  range  of  conditions  encountered  during  the  experiments 
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at  The  Pennsylvania  State  University,  a computer  program  that 
calculates  a polynomial  expression  to  fit  data  by  the  method  of 
least  squares  was  employed.  Ruth  and  Chant's  standard  correlation 
(83)  was  fitted  to: 

Kp  » 3.01  - 6.85(K^)  + 4.19(K^)^  , 

with  a reduced  chi-squared  value  of  0.83  x 10~^,  for  0.65  $ ^ 0.83, 

where  the  nomenclature  is  as  in  the  references  cited  above,  i.e.: 

H = total  insolation  on  a horizontal  surface  (measured  here) , 
Hq  =■  extraterrestrial  insolation  on  a horizontal  surface, 

Hjj  = diffuse  insolation  on  a horizontal  surface, 

S ■ 


and 


S - "'"o- 

_2 

Hq  varies  somewhat  (+  3%)  from  its  annual  average  of  1353  W m (51) 
because  of  the  annual  variation  in  the  sun-earth  distance.  Although 
the  plot  of  the  flux  over  a year  appears  as  a sinusoid  (2) , it  will 
be  represented  here  by  a linear  relationship,  to  sufficient  accuracy 
for  computing  K^. 
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2.  Natural  Convection  Within  Enclosed  Aperture 


During  most  experiments  the  receiver  aperture  was  covered 
with  a transparent  Mylar  film,  for  two  reasons:  to  reduce  convec- 
tive losses  from  the  hot  aperture  surfaces,  and  to  permit  a more 
accurate  mathematical  formulation.  Although  incorporation  of  a 
natural  convection  term  in  the  heat  loss  equation  for  the  first 
experimental  model  did  not  present  any  insuperable  difficulties, 
it  was  a source  of  large  error.  The  influence  of  the  wind  on  a 
complex,  tilted  shape  is  difficult  to  quantify  even  if  the  wind 
direction  and  velocity  are  accurately  measured. 

With  an  enclosed  aperture,  there  still  exists  the  possibil- 
ity of  enhanced  heat  loss  across  it  from  natural  convection  cell 
formation  within  it.  However,  the  presence  of  the  hot  Pyrex  at 
the  top  of  the  enclosure  tends  to  produce  a stable  system,  i.e., 
no  appreciable  air  motion. 

Catton  et^  (86)  solved  the  two-dimensional  equations  of 
motion  governing  natural  convection.  A decrease  in  heat  transport 
at  the  top  of  an  enclosure  was  observed,  when  it  was  heated  from 
above.  The  likely  cause  is  density  stratification,  creating  a 
stagnant,  constant  temperature  layer  adjacent  to  the  heated  top 
wall.  While  equations  were  developed  for  a finite  Prandtl  number 
rluid,  a solution  was  determined  only  for  a large  Pr  (i.e.,  high 
kinematic  viscosity)  fluid.  Examination  of  the  results  for  a slot 
of  aspect  ratio  (length/height)  of  20,  at  an  angle  of  tilt  = 60“ 

4 

shows  that  the  Nusselt  number  = 1 to  a Rayleigh  number  of  ca.  10  , 
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Ra  = Pr  g 6 AT/v^  , 


where  the  properties  of  the  air  trapped  within  the  aperture  are 

evaluated  at  T,,,  (37),  which  is  equal  to  (T,  + T )/2  here; 

tixm  ti  c 

S ■ 1/T  * 1/T  (38),  as  in  Chapter  II. 

Substituting  in  the  appropriate  values. 


- 143‘’C  ° (0.7)  (9.8  m s“^)  (143  + 273  ‘’K)“^(0.025  m)^ 


(261  - 25“C)/(26.9  x lO"^  m^s"^)^ 


8.4  X 10 


Rai  , 7gO(,  = (0.7)  (9.8)  (78  + 273)"^(0. 025)^(78  - 25)/(21.8  x 10“^)^ 


- 3.41  X 10  . 


The  effect  of  the  Transite  sides  has  been  ignored,  and  the  height  of 
the  aperture  (L  in  the  definition  of  Ra)  has  been  conservatively  taken 


as  2.5  cm. 


The  two  calculated  Ra's,  when  compared  with  the  above  analyti- 
cal and  experimental  results,  indicate  that  probably  1 < Nu  < 2 
(Nu  = 1 -►  pure  conduction).  Thus,  natural  convection  effects  are 
small,  and  given  the  somewhat  higher  conductivity  of  the  insulation 


A very  efficient  insulation  was  placed  around  the  Pyrex 
receiver.  Its  conductivity  varies  appreciably  with  temperature, 
and  a least-squares  polynomial  expression  was  fitted  to  the  data 
(35)  on  its  properties,  i.e.; 

k,  W cm  m"^  = 2.87  - 4.15  x 10“^(f)  + 1.76  x 10"^(f)^  , 

where  T is  the  mean  insulation  temperature,  ®K.  This  expression 
yields  a reduced  chi-squared  value  of  0.124  x 10 

Since  R^/R^  « (2.5  + 0.5  + 2.5)/(2.5  + 0.5)  * 1.83  is  less 
than  2,  an  arithmetic  average  insulation  area  can  be  an  acceptable 
approximation  for  use  in  computing  the  heat  transfer.  For  greater 
accuracy,  the  logarithmic  mean  area,  devined  by  Kreith  (40)  as: 

I ~ {k  - A.)/ln(A^/A,) 
o 1 0 1 

* 2Tr(5.5  - 3.0  cm)  (1  cm  length) /ln(5 . 5 cm/3.0  cm) 

2 

- 25.9  cm  /cm  receiver  length 


[• 
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will  be  used  here.  The  radial  temperature  drop  across  the  0.5 
cm-thick  Pyrex  is  very  small  and  is  ignored. 

j 

i 

4.  Radiation  Losses 

The  so-called  view  factor  characterizes  the  influence  of 
geometry  on  radiative  heat  exchange  between  two  surfaces.  The 
diffuse  view  factor  ^ from  finite  surface  A to  finite  surface  B 
is  defined  by  5zisik  (90)  as: 

F^_g  » (Radiative  energy  leaving  surface  A that  strikes  B directly)/ 

(Radiative  energy  leaving  surface  A in  all  j 

directions  in  the  hemispherical  space) 

For  a very  long  enclosure  where  each  surface  is  an  isothermal, 
diffuse  reflector  and  diffuse  emitter,  the  radiative  exchange  along 
the  longitudinal  axis  is  negligible  and  the  exchange  becomes  two- 
dimensional.  The  stagnation  condition  approximates  such  a situation 
quite  well  (Figures  22  and  23) . Under  flowing  gas  conditions  there 
is  an  axial  temperature  gradient,  but  it  is  too  small  (and  the  abso- 
lute temperatures  involved  are  too  low)  for  appreciable  axial  radiative 
exchange  to  be  present.  The  low  aspect  ratio  of  the  enclosed  aperture 
also  precludes  significant  axial  radiation  (see  following  discussion). 

Therefore,  the  approximation  of  essentially  two-dimensional  radiative 
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losses  from  the  aperture  surfaces  will  be  extended  to  all  cases. 
Hottel  and  Sarofim  (91)  present  a simple  method  for  calculating 
the  two-dimensional  view  factors  for  an  enclosure,  and  suggest  its 
use  in  approximating  a real  engineering  system. 

This  "crossed-string"  method  is  applied  to  a long  four-sided 
enclosure  as  follows  (90,91).  Imaginary  strings  are  stretched 
between  the  corners  (transverse  section)  and  the  view  factor  from 
surface  A to  B becomes: 

g » [(sum  of  lengths  of  2 crossed  strings)  - 
(sum  of  lengths  of  2 uncrossed  strings) ]/2  x length  of  side  A 

This  is  illustrated  on  Figure  24,  which  also  shows  the  important 
receiver  dimensions.  Thus, 

Vb  ■ "4  * S>  - • ’•A 


or 


Vb  * 5.0)]/2(2.5)  = 0.24  = Fg_^  , 


by  reciprocity.  Similarly,  F * F = 0.68,  F « 0.1, 

r “A 

Fp_^  = 0.80,  F^_p  » 0.08,  and  Fg_p  = 0.28. 
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5 . End  Effects  and  One-Dimensional  Approximation 

The  calculation  is  performed  only  at  the  center  (axially)  of 
the  device  where  the  temperature  profile  is  flat,  and  the  cooler 
ends  are  both  distant.  The  effects  of  the  axial  dimension  are 
therefore  negligible.  On  the  further  approrimation  of  peripherally 
isothermal  conditions  (at  the  Pyrex  exterior,  by  use  of  an  average 
temperature) , the  problem  reduces  to  one  dimension,  with  the  view 
factor  effectively  incorporating  two  dimensions  into  one. 


6 . Energy  Balance  | 

j 

The  calculation  is  again  (as  in  Chapter  II)  based  on  a 1-cm  ] 

length  of  receiver.  The  total  concentrated  flux  per  cm  of  receiver 

length  is:  i 

( 

qg  = (P^) (Wj^) (no.  mirrors) (1  - pointing  error) (1  - Pp) (1  - • F^)  | 


^^Mylar 


)(meas.  direct  perpendicular  flux  at  mirror)  . 


i 


This  effectively  delimits  one  control  volume  surface  as  the  Pyrex 
aperture.  The  tilt  of  individual  strips  is  not  taken  into  account, 
since  the  maximum  tilt  is  only  18°.  Substituting  in  the  proper  values. 


qg  = (0.80) (1.59  cm)(40)(l  cm  length) (10  ^ m^  cm  ^) (1  - 0.05) 


(1  - 0.07) (1  - 0.13  • 0.28) (0.85) (flux,  W m"^) 
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= (3.67  X 10  ) (direct  perpendicular  flux),  W,  per  centimeter  of 

length,  which  is  available  to  the  Pyrex  receiver.  At  steady-state 
this  is  counterbalanced  by  the  thermal  losses,  where 


qg  = a(e  F A T^)p  + a(E  F A A(Tp  - T^) /insulation  tk.  , 


and  where  the  first  RH  term  represents  the  raaiative  heat  loss  from 

the  glass  window,  the  second  term  is  the  radiative  loss  from  the 

aperture  sides,  and  the  last  term  is  the  loss  by  conduction  through 

the  insulation.  It  is  assumed  that  T = (T„  + T )/2  = (T19  + f )/2, 

w P e e 

where  = (T3  + T8)/2,  provides  an  adequate  measure  of  the 
temperature  of  the  aperture  walls.  This  (linear)  assumption 
probably  somewhat  underestimates  the  actual  loss  (because  of  the 
exponential),  but  this  is  mitigated  by  the  poorer  view  factor  of 
the  hotter  (i.e.,  nearer  the  glass)  parts  of  the  walls.  Radiative 
exchange  in  the  form  of  wall  reradiation  to  the  glass  is  ignored. 

The  view  factor  is  very  small  (0.08),  and  a conservatively  high  Ep 
for  the  glass  was  selected  (1.00),  partly  in  order  to  compensate 
for  any  absorbed  reradiation.  Tp,  the  average  peripheral  glass 
temperature,  is  taken  as  (T16  + T7)/2;  i.e.,  the  glass  temperature 
is  assumed  to  vary  in  linear  fashion  from  its  hottest  points,  where 
the  graphite  rests  against  it  (taken  as  ca.  = T16,  at  9 and  3 
o'clock)  to  the  coolest  points  (T19,  which  = Tp,  and  T7)  at  6 and 
12  o'clock.  The  graphite  temperature  at  the  axial  center  of  the 


device  was  not  measured,  but  is  probably  only  slightly  higher  than 
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T16,  the  gas  temperature  there  (see  Figures  22  and  23).  Given  its 
high  thermal  conductivity,  the  graphite  ought  to  be  isothermal,  in 
the  transverse  direction. 

The  measured  stagnation  temperatures  on  March  23,  1977,  at 
12:00  p.m. , were:  T7  - 297‘’C,  T19  = 261,  T16  = 388,  T3  =•  25,  and 
T8  * 30.  That  day  was  very  windy,  which  lowered  the  T19  (i.e., 
glass  window)  reading,  the  thermocouple  being  relatively  exposed. 
Thus  T19,  which  is  usually  “ T7,  is  here  nearly  U0°C  lower,  and 
will  not  be  used  in  the  calculation.  Substituting  the  appropriate 
values , 

(872  W m~^)(3.67  x lO"^)  ^ (5.67  x 10“®  W m“^  “K~^) (1.00) (0.80) 
(297  + 273'’K)^(2  cm^)  (lO"^  m^  cm“^) 

+ (5.67  X 10"®  W m"^  “k"^)  (0.94)  (0.68)  (5  cm^)  (lO"'^  m^  cm"^) 
[(297  + 28)/2  + 273‘’K]^ 

+ (4.66  W cm  m"^  '’c"^)(2.6  x 10~®  m^)[(388  + 297)/2  - 
(30  + 25)/2]/2.5  cm^. 


I 
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3.2  W 3.14  W,  which  is  quite  close. 
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If  one  attempts  to  calculate  the  convective  loss  from  the 
receiver  insulation  to  the  ambient  environment,  a large  error 
results.  For  example,  with  the  above  temperatures  and  a wind 
velocity  = ca.  8 m s ^ (measured,  at  the  Meteorological  Station), 
the  calculated  convective  heat  loss  from  the  external  skin  is 
2.47  W.  A well-known  empirical  equation  (8)  was  used  in  obtaining 
this  value;  a cross-flow  condition  was  assumed. 

Since  the  true  skin  loss  is  seriously  overpredicted,  this 
loss  has  not  been  calculated  in  subsequent  analyses.  That  loss 
should,  of  course,  be  equal  to  the  conduction  loss  through  the 
receiver  insulation;  the  latter  is  readily  and  accurately  deter- 
mined. Because  skin  temperatures  are  generally  only  ca.  25 °C 
above  ambient  (see  Table  8) , it  is  evident  that  the  convective 
loss  from  the  skin  is  not  a limiting  factor  on  receiver  performance. 

Selecting  another  date  at  random  (but  near  the  equinox,  for 
minimal  aperture  shadowing/end  loss),  as  March  9,  1977  at  12:00  p.m., 
with  T7  - 281,  T19  = 263,  T16  = 367,  T3  = 34,  and  T8  = 40,  and 
again  using  T7  exclusively  results  in 

3.17  W = 2.86  W . 

The  losses  (RHS)  are  a bit  low  for  the  flux,  buc  the  flux  reading 

discrepancy  between  the  two  measurements  was  somewhat  high 
_2 

(90  W m ).  On  another  day,  on  which  the  two  flux  mtasurements 
were  nearly  equal,  as  on  April  1,  1977,  at  11:45  a.m.,  with 


T7  - 288,  T19  = 281,  T16  » 388,  T3  - 44,  and  T8  - 45,  one  calculates 


The  above  calculations,  using  the  actually  measured  stagna- 
tion temperatures,  show  that  the  main  heat  loss  mechanisms  have 
been  correctly  characterized.  Agreement  between  the  available 
solar  flux  and  the  sum  of  the  losses  is  good,  always  within  10%. 
Considering  the  assumptions/approximations  which  have  been  made, 
and  experimental  error,  this  agreement  is  adequate  for  the  purpose 
of  calculating  some  of  the  losses  under  gas  flow  conditions. 

C.  Flowing  Gas  Calculations  (Axial  Bulk  Temperature  Profile) 

1.  Introduction 

I 

At  steady-state,  the  three  parameters  that  almost  wholly 

determine  the  axial  temperature  profile  of  the  fluid  (i.e.,  T,  ) 

b 

along  the  receiver  are  gas  inlet  temperature,  gas  flow-rate,  and 
available  solar  flux.  The  mathematical  model  that  faithfully 
represents  the  experimental  device  should  succedd  in  predicting  the 
temperature  profile  from  these  initial  conditions. 

Two  complications  are  now  encountered: 

(1)  The  axial  temperature  gradient  present  with  gas  flows 
provides  a thermal  driving  force,  so  that  axial  conduction  losses 
can  no  longer  be  Ignored. 

(2)  An  appreciable  thermal  resistance  is  offered  by  the  laminar 
flow  of  the  gas,  resulting  in  a significant  transverse  temperature 


gradient . 


These  mechanisms  must  be  elucidated  and  incorporated  into 
the  mathematical  model.  Because  this  flowing-gas  model  considers 
much  of  the  receiver's  (relatively  short)  total  length,  end 
effects/losses  can  exert  an  influence  and  must  be  thoroughly 
examined.  Thus,  the  spectrum  of  heat  transfer  mechanisms  that 
will  now  be  investigated  in  detail  includes: 

(1)  axial  conduction; 

(2)  radiative  end  loss; 

(3)  effects  of  shape  (semi-circular  receiver)  and  tilt; 

(4)  development  of  velocity  and  temperature  profiles; 

(5)  forced  convection; 

(6)  natural  convection. 

Before  proceeding  further,  some  general  remarks  on  the 
approach  taken  here  are  in  order.  A closed-form  solution  to  the 
problem  under  consideration  is  clearly  not  available.  The  full  set 
of  mass,  momentum,  and  energy  conservation  equations,  under  the 
conditions  of  developing  flow  and  with  buoyancy  effects  included, 
has  not  been  solved  at  present,  even  by  employing  numerical  tech- 
niques and  adopting  constant  boundary  conditions  and  other 
simplifying  assumptions  (92) . Hence,  it  appeared  pointless  to 
attempt  this  approach,  particularly  in  view  of  the  coupled  boundary 
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conditions  and  end  effects  in  the  experimental  device.  The  solution 
presented  here  has  been  constructed  by  utilizing  investigations  of 
particular  and  separated  aspects  of  the  total  problem.  As  these 
interact  to  some  extent,  assumptions  had  to  be  made  that  would 
take  this  interaction  into  account  where  necessary. 

2.  Axial  Heat  Conduction 

The  absorber  and  the  glass  receiver  are  relatively  thick,  and 
the  absorber  has  a high  thermal  conductivity.  Axial  heat  conduction 
by  these  two  parts  constitutes  a significant  loss  at  the  inlet, 
especially  at  low  flows  (high  axial  temperature  gradient) . It  can 
be  computed  by  an  iterative  technique  from  the  temperature  profile, 
which  is  first  calculated  without  inclusion  of  this  loss.  The 
conductivities  are  taken  as  constant;  the  absorber  conductivity 
is  not  accurately  known,  while  the  Pyrex  has  a relatively  low 
conductivity. 

3.  Radiation  End  Loss 

Both  the  graphite  absorber  and  the  glass  receiver  have  high 
emittances  at  the  temperatures  attained  by  the  collector.  There  is 
therefore  some  heat  radiated  from  the  receiver  test  section  to  the 
cooler  end  pieces,  and  the  magnitude  of  the  loss  must  be  estimated. 

Usiskin  and  Siegel  (93)  solved  the  applicable  integro- 
differential  energy  equation  by  three  methods,  but  only  for  the 
limiting  case  of  pure  radiation  in  a tube.  Siegel  and  Perlmutter 
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(94,95)  used  part  of  this  analysis  in  approximating  the  geometrical 
shape  factors  with  a separable  kernel,  while  considering  combined 
convective  and  radiant  heat  transfer  in  a tube.  This  approach 
changes  the  energy  equation  to  a second-order  ordinary  differential 
equation,  which  can  be  readily  solved  by  numerical  methods.  Both 
tabulated  and  graphical  solutions  are  provided  (94,95),  for  a 
variety  of  conditions. 

Thorsen  (96) , and  Thorsen  and  Kanchanagom  (97)  performed 

experiments  with  air  in  narrow,  long  (L/D  = 50)  tubes  at  laminar 

flows  and  high  heating  rates  (e.g.,  T /T  . , * 5 at  the  entrance 

w gas  inlet 

and  ca.  9.5  near  the  exit,  while  T . ji..  = 4.5;  the 

temperatures  are  in  °F,  with  an  inlet  temperature  of  500°?.  Axial 

wall  heat  conduction  and  temperature  profile  development  were  studied 

simultaneously.  The  results  are  close  to  the  theoretical  predictions 

of  Siegel  and  Perlmutter  (94,95):  internal  wall  radiation  exchange 

with  the  tube  ends  is  limited  to  ca.  the  first  and  last  20%  of  tube 

length  (L/D  = 50),  with  the  greater  loss  near  the  exit.  The  influence 

* 

on  the  Nusselt  number  over  the  first  80%  of  the  tube  length  is  slight 
(96,97).  The  high  Nu  resulting  from  flow  development  near  the 
entrance  reduces  the  wall  temperature,  and  decreases  the  radiative 
end  loss  there. 

Now  an  experiment  with  high  inlet  and  outlet  temperatures  and 
a high  solar  flux  is  selected,  as  on  March  19,  1977,  1:30  p.m. 

(Table  8),  where  T14  = 120'’C,  T17  = 281,  and  the  flux  was  872  W m ^ 
with  a flow  rate  of  32.3  liter  min  Here  the  temperatures  and 
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other  parameters  at  the  32  and  123-cm  axial  positions  are  chosen 
as  representative  of  the  receiver  inlet  (0-cm  position)  and  outlet 
(152-cm  position),  respectively,  and  are  used  in  the  calculations 
that  follow.  The  reason  for  this  approximation  is  that  the  32  and 
123-cm  parameters  are  well-known  from  computer  results.  This  approx- 
imation is  conservative  at  the  receiver  inlet,  because  the  temperature 
at  the  32-cm  position  is  considerably  higher  than  that  at  0-cm;  at 
the  receiver  outlet,  the  123-cm  and  153-cm  temperatures  differ  only 
by  a small  amount,  so  the  approximation  is  quite  valid  there  (Table  8). 
The  influence  of  four  parameters,  as  defined  in  (95),  must  be 
examined : 

a.  Stanton  number 

The  Stanton  number,  S,  is  defined  as  4 Nu/Re  • Pr.  Here, 

Re  * 554,  Ra  = 354  -*■  Re  • Ra  = 1.96  x 10^  (taken  from  computer 

results;  see  Appendix  II)  and  (Nu  ) , -o/ (Nu,)„o  = 0.27 (Re  • Ra)^'^^ 

c 45  r U 

(see  discussion  on  combined  convection,  below) . Thus 

(Nu^)^go/ (Nu^)qo  = 1.90,  indicating  a 90%  enhancement  of  heat 

transfer  by  natural  convection  effects.  Also,  z"*"  = z/D,  Re_  Pr  = 

h D 

32  cm/(2.86  cm) (554) (0.67)  = 0.030,  resulting  in  a (Nu^)qo  = 5.2 
(see  discussion  on  developing  flow,  below),  versus  4.36  without 
development  effects,  in  a circular  tube.  Therefore  the  actual 
Nusselt  number  =*  (5.2)  (1.90)  (4.09/4.36)  = 9.3,  where  the  fraction 
corrects  for  a lower  Nu  (i.e.,  4.09)  in  a semicircular  sector  tube. 

The  Stanton  number  can  now  be  computed: 


4 


S = 4(9. 3)/(554) (0.67)  = 0.10 
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When  this  value  is  compared  with  the  results  of  Perlmutter  and 
Siegel  (95),  it  is  evident  that  radiant  losses  are  low.  Although 
only  the  range  of  0.005  ^ S ^ 0.02  is  shown  ((95),  Figure  4(a)),  the 
radiative  end  loss  is  relatively  insensitive  to  variation  in  the 


Stanton  number  over  this  range.  The  loss  is  negligible  in  the  60% 
central  length  of  tube;  extrapolation  to  S = 0.10  should  yield 


similar  results.  Near  the  receiver  exit,  S would  be  ca.  0.05  because 
of  a lower  Nu  at  that  position.  Also,  the  results  (95)  are  given  for 
a tube  with  L/D  = 50,  whereas  the  experimental  device  had  L/D  = 152 
cm/ca.  1.5  cm  -►  L/D  > 100.  Furthermore,  the  wall  temperature  shown 
(95)  is  approximately  double  the  measured  wall  temperature  here. 

Thus,  the  receiver's  length  and  (relatively)  low  wall  temperature 
serve  to  decrease  the  Stanton  number  influence  on  the  theoretically 
predicted  (95)  radiative  end  losses,  for  the  date  selected.  Since 
the  variation  in  the  Stanton  number  for  all  experiments  is  small, 
it  can  be  assumed  that  the  radiant  loss  over  the  central  60%  of 
the  receiver  tube  can  be  neglected. 

b.  Dimensionless  heat  transfer  coefficient 
The  dimensionless  heat  transfer  coefficient,  H,  is  defined 
as  H = (h/q) (q/o)^” . Here,  q ■ (872  W m ^)(3.67  x 10  ^ m^  cm  ^) 
(0.53)  = 1.7  W cm  ^ length  (the  0.53  is  the  calculated  Solar  Utiliza- 
tion Factor,  or  SUF,  defined  here  as  9jj/9g;  the  above  calculation 

yields  q^^) . h = Nu  • k/D^  = (9.3)  (2.23  x 10“^  W cm"^  °k“^)/2.86  cm 
-4  -2  -1 

* 7.25  X 10  W cm  . H can  now  be  calculated: 
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„ -1  o„-l 


H - [(7.25  X 10  ” W cm  ^ *k“-^)(25.7  cm^  cm"^)/1.7  W cm“-^] 


-1. 


(1.7  W cm"^/25.7  cm^  cm"^  • 5.67  x 10“^^  W cm^  Oj^-ljO.25  _ 

This  high  (95)  convective  heat  transfer  coefficient  reduces  the 
potential  radiative  loss  and  again  limits  it  to  20%  of  the  tube 
length,  at  the  inlet  and  outlet  ((95),  Figure  4(b)).  The  possible 
variation  (experimental)  in  H is  small,  so  this  conclusion  holds  for 
all  experiments. 

c . Dimensionless  gas  inlet  temperature 

The  dimensionless  gas  inlet  temperature,  T . , is  defined 

8»  i 

0 . 25 

Tg  ^ = ((j/q)  * T^,  where  T^  is  the  gas  inlet  temperature.  Here, 


,'12  „ ._-2  o„-l  , .2  „ -1,0.25 


Tg,i  = (5.67  X 10  ^ W cm  ^ ^ • 25.7  cm^  cm"-"/1.7  W cm"-^) 


(120  + 273°K)  =•  1.2 


As  before,  the  effect  of  this  parameter  is  small,  and  it  exerts 

an  influence  only  over  the  first  and  last  20%  of  tube  length  ((95), 

Figure  4(c)).  Because  of  the  exponential,  little  variation  is 

possible  in  T , over  the  range  of  experimental  measurements. 
s>  ^ 

d.  Dimensionless  reservoir  temperature 

T , the  dimensionless  reservoir  temperature  at  the  exit 
r , e 

(assumed  in  (95)  to  be  equal  to  the  gas  exit  temperature) , is  defined 
by 
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not  been  included  in  the  mathematical  formulation.  On  purely 
theoretical  grounds  incorporation  of  the  radiative  interchange  is 
desirable,  since  only  in  this  manner  can  the  peripheral  wall  tempera- 
ture distribution  be  determined.  A single,  average  peripheral  wall 
temperature  will  therefore  be  used  in  the  present  formulation,  to 
determine  the  heat  transfer  both  to  the  gas  and  to  the  ambient 
environment.  This  approximation  slightly  overestimates  the  losses 
(see  Figure  25) . 

4.  Receiver  Shape 

Of  all  the  simple  geometric  shapes,  a circular  cross  section 
in  a tube  yields  the  highest  convective  heat  transfer  coefficient. 

Exact  analytical  solutions  for  the  coefficient  are  extant.  Under 
the  conditions  of  fully  developed  laminar  flow  of  a constant-property 
fluid,  with  constant  peripheral  wall  temperature,  the  Nusselt  number 
for  a constant  axial  heat  input  is  4.36.  Under  the  same  conditions, 
but  with  a constant  axial  wall  temperature  instead  of  a constant 
axial  heat  flux,  Nu  = 3.64  (98).  The  lower  Nu  results  from  the 
decreasing  temperature  difference  between  the  wall  and  fluid,  as  the 
latter  proceeds  through  the  tube. 

A semicircular  cross  section  evinces  poorer  overall  heat 

transfer  because  of  relatively  stagnant  flow  in  the  corners.  With 
t 

conditions  as  before,  and  a constant  axial  heat  flux,  a closed-form 
analytical  solution  for  Nu  is  available.  Sparrow  and  Haji-Sheikh  (99) 
obtained  a solution  of  the  momentum  and  energy  equations  from  analogous 
torsion  and  thin-plate  deflection  solutions,  respectively,  in  solid 
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mechanics  (98) . These  equations  have  also  been  solved  by  an  alternate 
approach,  employing  a computer  routine  that  constructs  orthonormal 
functions  (100).  Both  solutions  gave  identical  results:  Nu  = 4.09. 
Solutions  for  circular  sector  ducts  with  an  opening  angle  < 180° 
are  in  excellent  agreement  with  the  experimental  data  for  air  (101) . 

As  the  circular-sector  opening  angle  approaches  180°  (a  semicircle), 
two  effects  are  noted:  an  increasing  symmetry  of  the  velocity  and 
temperature  profiles  in  the  half-angle  plane  (99,100),  and  a 
decreasing  difference  in  the  Nusselt  numbers  obtained  for  the  cases 
of  peripherally  varying  heat  flux  and  peripherally  varying  wall 
temperature  (100) . 

Before  the  above  results  can  be  used  with  confidence,  the 
three  assumptions  that  were  made  in  the  solutions  must  be  justified. 
These  are: 

a.  Constant  peripheral  wall  temperature 

The  Pyrex  and  graphite  are  relatively  thick,  and  the  graphite 
has  a high  thermal  conductivity.  The  graphite  strongly  reradiates 
the  solar  flux  incident  on  it;  any  portion  of  the  incident  flux 
striking  the  insulation  is  diffusely  reflected.  The  Pyrex  also 
reradiates  strongly.  With  these  high  heat  transfer  rates  in  the 
transverse  plane  of  the  receiver,  any  variations  in  the  peripheral 
temperature  distribution  should  be  small.  The  peripheral  temperature 
variation  was  not  measured  per  se,  but  comparison  of  T7,  T19,  and 
T16  (this  last  approximates  the  graphite  temperature  near  the  axial 
center  of  the  receiver)  under  flow  conditions  gives  some  idea  of  its 
magnitude  (Figure  25  and  Table  8) . 
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Reynolds  (102)  considered  a circular  tube,  as  described 
previously.  The  governing  differential  energy  equation  is  linear. 
This  permits  superposition  of  simpler  solutions  for  the  solution 
of  a complex  circumferential  heat  flux  variation.  For  example,  for 
q ” qQ  (1  + 0.1  cos  a),  describing  a + 10%  sinusoidal  heat  flux 
variation,  the  computed  circumferential  temperature  variation 
Tmay/^  “ 1.22,  with  a local  (i.e.,  circumferentially  varying) 

Nusselt  number.  If,  instead,  the  computation  is  repeated  with 
Nu  = 4.36  (constant  circumferential  wall  temperature  assumed), 

T /T  = 1.10,  fcr  an  error  of  (1.22  - 1.10)/1.22  = 10%.  The 

in3X 

experimental  peripheral  temperature  variation  observed  (again, 
compare  T7,  T19,  and  T16)  is  generally  less  than  the  T^^^/T  = 1.22, 
above,  so  any  error  introduced  by  this  approximation  should  be  small. 
In  any  case,  an  average  wall  temperature  is  used  in  the  flow  calcu- 
lations, reducing  this  potential  error  to  negligible  proportions. 

b.  Constant  axial  heat  flux 

Although  the  solar  flux  entering  the  receiver,  per  unit 
length,  is  constant  axially,  the  actual  flux  input  into  the  flowing 
fluid  is  not.  As  the  fluid  heats,  losses  mount  so  that  eventually  an 
asymptotic  condition  is  reached  with  no  net  heating  taking  place. 

Over  the  central  91-cm  portion  the  temperature  rise  of  the  heat 
transfer  surfaces  is  fairly  linear  (Figure  25);  however,  the  con- 
vective coefficient  varies.  Still,  this  condition  approximates  more 
closely,  over  the  studied  length,  a constant  axial  heat  flux  boundary 
condition  than  a constant  axial  wall  temperature  (the  only  other  case 
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extensively  studied  (103)).  The  error  introduced  by  this  simplifi- 
cation is  difficult  to  quantify.  The  effect  of  a decreasing  heat 
flux  in  the  flow  direction,  as  in  the  experimental  device,  is  to 
yield  a somewhat  lower  heat  transfer  coefficient  (i.e.,  higher 
actual  wall  temperature  ^ greater  losses) . Kays  (104)  shows  that 
the  effect,  qualitatively,  for  a sinusoidally  (0  -*■  maxi.mum  -*•  0) 
varying  axial  heat  flux  is  relatively  minor.  It  will  be  ignored 
here. 


c . Fully  developed  flow 

A fluid  entering  a pipe  at  uniform  temperature  and  velocity 
requires  many  tube  diameters  before  the  fully  developed  velocity 


and  temperature  profiles  become  established;  the  profiles  are 
invariant  thereafter  if  either  the  heat  flux  or  the  wall  tempera- 
ture are  constant  axially.  For  laminar  flow,  the  approximate  length 
necessary  for  a fully  developed  velocity  profile  is  equal  to 
D • (Re/20)  (105).  For  a gas  (Pr  = 0.7),  the  temperature  profile 
develops  somewhat  more  rapidly.  Thus,  a significant  amount  of 
profile  development  takes  place  over  the  receiver  length  and  must 
be  evaluated  separately. 


5.  Development  of  Velocity  and  Temperature  Profiles 


Although  the  classical  analysis  of  temperature  profile  develop- 
ment in  laminar  flow  was  made  by  Graetz  in  1885  (106) , little 
additional  work  on  this  topic  was  performed  until  the  mid-1950's. 
Sellars  et  al.  (107)  extended  the  classical  solution  to  allow 
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consideration  of  more  complex  boundary  conditions.  Kays  (108) 
was  the  first  to  treat  this  problem  in  detail. 

An  exact  analytical  solution  for  a developing  velocity  field 

ft 

in  a circular  tube  has  not  been  derived  (109).  Hence,  solutions 
of  the  temperature  and  velocity  fields  in  forced  convection  are 
derived  from  approximate  analytical  or  numerical  techniques. 

Langhaar  (110)  obtained  a solution  for  the  developing  axial 
velocity  profile  (in  the  radial,  or  transverse,  plane)  by  linear- 
izing the  nonlinear  inertia  terms  in  the  momentum  equation.  Kays 
(108)  Introduced  this  velocity  profile  into  the  energy  equation 
and  solved  it  by  a finite  difference  (i.e.,  numerical)  method. 
Nusselt  numbers  versus  length  are  given  for  constant  axial  wall 
temperature  and  constant  axial  heat  flux,  and  for  both  fully 
developed  and  developing  velocity  profiles,  for  the  case  of  a 
constant-property  fluid  in  a round  tube.  A uniform  circumferential 
wall  temperature  was  assumed.  Sparrow  and  Lin  (111)  refined 
Langhaar 's  (110)  approximations,  but  obtained  a velocity  field 
that  closely  matches  that  of  the  earlier  solution. 

The  boundary  layer  equations  can  be  solved  with  fewer  assump- 
tions than  were  made  above  if  axial  transport  of  vorticity  is 
neglected,  and  the  pressure  is  taken  as  a function  only  of  axial 
length.  Finite  difference  techniques  were  employed  by  Hornbeck 
(112),  Christiansen  and  Lemmon  (113),  and  Schmidt  and  Zeldin  (114) 
to  obtain  the  developing  velocity  profile.  The  results  are  in  good 
agreement  with  each  other  and  with  Langhaar 's  profiles,  establishing 
confidence  in  their  validity. 
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Ulrichson  and  Schmitz  (115)  considered  heat  transfer  with 
simultaneous  development  of  temperature  and  velocity  profiles  of 

i 

an  incompressible  viscous  (Newtonian)  fluid  (Pr  = 0.7,  i.e,  a gas). 

The  boundary  conditions  of  constant  axial  heat  flux  and  constant 
axial  wall  temperature  were  used.  The  solution  represents  a refine-  I 

ment  of  the  work  of  Kays  (108) , who  ignored  the  radial  velocity 
component.  Langhaar's  (110)  axial  velocity  profile  was  employed, 
together  with  the  radial  profile  obtained  form  the  continuity 
equation,  in  a numerical  solution  of  heat  transfer  rates.  The 
results  show  that  the  local  Nusselt  number  near  the  entrance  is 
overestimated  by  about  6%  if  the  radial  velocity  is  neglected,  for 
the  case  of  constant  axial  heat  flux. 

Manohar  (116)  solved  the  non-linear  momentum  equation  by  an 
iterative  technique  to  obtain  the  velocity  profile.  This  velocity 
distribution  was  inserted  into  the  energy  equation,  which  was  solved 
numerically  for  the  temperature  profiles.  The  computed  heat  trans- 
fer closely  matches  that  of  Ulrichson  and  Schmitz  (115)  and 
experimental  data.  Very  recently  Hong  and  Bergles  (117)  treated 
a semicircular  tube  problem  by  the  methods  of  successive  over- 
relaxation and  finite  differences.  However,  their  heat  transfer 
investigation  was  limited  to  the  condition  of  a fully  developed 
velocity  profile. 

The  work  of  Manohar  (116)  is  the  most  recent  and  apparently 
provides  the  most  accurate  results,  and  will  be  used  here.  A least- 
squares-fit  equation  to  Manohar 's  Nusselt  numbers  for  a constant 
axial  heat  flux  condition  is: 
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Nu  = 4,18  + 3.17  X lO'^Cz'*’)"^  - 3.67  x lO'^Cz''')"^ 

+ 1.91  X lO'^Cz"^)"^ 

• _2 
for  a reduced  chi  squared  value  of  1.81  x 10 

An  additional  conunent  can  be  made  now,  on  the  constant  axial 
heat  flux  assumption  discussed  previously.  The  boundary  condition 
of  the  receiver,  although  closer  to  the  constant  axial  flux  case, 
is  actually  delimited  by  this  and  the  constant  axial  wall  tempera- 
ture case.  Over  the  experimental  range,  Kays'  (108)  values  of 
Nusselt  number  for  the  constant  wall  temperature  case  were  only 
ca.  15%  lower  than  for  the  constant  heat  flux  case.  Manohar's 
(116)  constant  wall  temperature  solution  deviates  appreciably  from 
that  of  Kays:  Nu  for  the  constant  wall  temperature  is  ca.  30% 
lower  than  for  constant  heat  flux.  Therefore,  use  of  the  more 
recent  data,  under  the  above  assumption,  carries  a somewhat  higher 
risk  of  error.  At  any  rate,  the  error  from  this  assumption  is 
probably  <<  30%  of  the  Nusselt  number.  The  computed  Nusselt  number 
is  used  to  calculate  the  difference  between  the  gas  and  wall  tempera- 
tures, with  losses  determined  from  the  latter.  Thus,  the  axial 
temperature  profile  is  relatively  insensitive  to  an  error  in  the 
Nusselt  number. 

6 . Combined  Forced-Free  Convection 

Fluid  motion  may  be  on  occasion  strongly  influenced  by  the 
presence  of  gravity,  with  the  magnitude  of  the  induced  body  force 


If 
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field  proportional  to  the  local  fluid  density.  The  classical 
example  is  free  convection,  where  fluid  motion  in  an  otherwise 
quiescent  fluid  occurs  as  a consequence  of  distortions  in  the  body 
force  field  caused  by  density  gradients.  At  the  other  extreme  is 
the  type  of  flow  referred  to  as  (pure)  forced  convection,  whose 
behavior  is  unaltered  by  gravity.  Most  real  forced  flows  are 
influenced  to  some  extent  by  the  gravity  force.  When  the  gravita- 
tional effect  is  significant,  this  flow  is  classified  as  combined 
forced-free  convection.  A measure  of  the  relative  magnitudes  of 
forced  and  free  convection  can  be  obtained  by  comparing  the  Reynolds 
number  with  the  Rayleigh  number,  as  these  are  defined  in  the  nomen- 
clature. When  the  Rayleigh  number  is  of  the  same  order  of  magnitude 
as  the  Reynolds  number,  free  convection  effects  are  significant  and 
cannot  be  ignored  in  a precise  heat  transfer  analysis  (118) . 

Combined  free  and  forced  laminar  convection  has  been  studied 
extensively  in  the  past  only  for  the  two  limiting  cases  of  vertical 
and  horizontal  tubes.  The  paucity  of  inclined-tube  studies  obligates 


a discussion  of  the  vertical  and  horizontal  results. 

a.  Combined  convection  in  vertical  tubes 
The  orientation  of  the  tube  axis  is  important  in  combined 
forced  and  free  convection.  Combined  convection  in  vertical  tubes 
can  be  considered  a special  case  of  combined  convection  in  inclined 
tubes.  In  vertical  tubes  the  buoyancy  forces  are  parallel  to  the 
direction  of  forced  motion  and  some  form  of  axial  symmetry  (rotational 
symmetry,  for  the  round  tubes  generally  studied)  of  the  velocity  and 
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temperature  fields  is  retained,  as  in  pure  convection.  This 
considerably  simplifies  the  heat  transfer  analysis,  and  approximate 
analytical  solutions  are  possible,  just  as  for  the  forced  convection 
condition  (119) . 

Most  of  the  theoretical  studies  posited  constant  axial  heat 
input,  constant  peripheral  wall  temperature,  fully  de^'eloped  flow, 
the  usual  flow  idealizations,  constant  (both  transverse  and  axial) 
fluid  properties,  and  laminar  flow.  Only  the  density,  and  that 
solely  in  the  buoyancy  terms  of  the  momentum  equation,  is  taken  as 
variable.  The  last  is  a common  simplification  of  the  governing 
equations,  and  is  considered  acceptable  when  the  density  is  prima- 
rily a function  of  temperature  (120) . Iqbal  and  Stachiewicz  (121) 
considered  a fully  variable  density  in  the  analytical  expressions 
and  determined  that  this  has  an  insignificant  effect  on  the  resultant 
heat  transfer  coefficients  (Nusselt  numbers) . 

Hallman  (122)  performed  an  analysis  and  compared  it  with 
experimental  data  obtained  by  others;  his  theoretical  Nusselt 
numbers  were  about  20%  higher  than  the  experimental  ones.  Tao  (123) 
formulated  a solution  in  terms  of  Bessel  and  associated  functions, 
and  found  the  heat  transfer  coefficients  for  circular  and  semi- 
circular tubes  to  be  about  the  same.  Iqbal  ^ aJ . (124,125,126, 

127)  obtained  solutions  by  the  methods  of  finite  differences  (124), 
variational  calculus  (124,125,126),  and  a point  matching  method  in 
terms  of  a series  containing  Bessel  functions  (127).  Aggarwala  and 
Iqbal  (128)  utilized  analogous  solutions  from  membrane  vibration 
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problems.  The  duct  shapes  considered  (124-128)  included  triangular, 

rectangular,  rhombic,  and  others;  Nusselt  numbers  were  obtained  for 

0 < Rayleigh  number  < 10^  (based  on  D.  , i.e.,  eight  times  that  used 

n 

in  this  thesis).  It  was  found  that  as  the  Ra^  increases,  the  Nusselt 

number  becomes  relatively  insensitive  to  the  duct  shape.  When 
3 

Ra^j  > 10  , all  solutions  converge  rapidly  (124,127).  Incorporation 
of  a wall  (i.e.,  heat  transfer  surface)  of  finite  thickness  and 
peripheral  (axial  ignored)  conductivity  into  the  analysis  indicated 
that  the  Nu  decreases  as  the  peripheral  wall  heat  conductivity 
decreases  (125) . 

Nayak  and  Cheng  (129)  solved  the  applicable  equations  by  a 
finite  element  method,  with  results  for  triangular  and  square  ducts 
as  above.  Muntjewerf  e£  ajL.  (130)  performed  experiments  with  two 
different  oils  in  rectangular  ducts,  but  only  at  very  low  Reynolds 
numbers  (Re  ^ 100) . 

A more  complex  treatment,  including  developing  flow  and  tempera- 
ture profiles,  was  given  this  problem  by  a number  of  researchers. 
Wors«5e-Schmidt  and  Leppert  (131)  solved  the  boundary-layer  equations 
via  a finite  difference  scheme,  and  found  only  slight  improvement  in 
heat  transfer  due  to  free  convection  in  the  developing  flow  region. 
Lawrence  and  Chato  (132)  used  a finite  difference  method  to  solve 
the  equations,  and  also  performed  experiments  with  water.  The 
theoretical  and  experimental  results  correlate  well;  to  obtain  this 
good  correlation,  both  the  density  and  the  viscosity  had  to  be  consid- 
ered as  nonlinear  functions  of  the  temperature.  From  the  results,  it 


was  concluded  that  the  temperature  and  velocity  profiles  never  become 
fully  developed. 

Experiments  with  nitrogen  gas  at  high  temperatures  (900®K) 

3 3 3 

over  10  < Re  < 3.5  X 10  and  Grashof  number  (based  on  D and  AT  = 

4 5 

T - T,  ) over  3 x 10  to  8 x 10  were  performed  by  Biggs  and 

W u 

Stachiewicz  (133) . A comparison  of  the  experimental  Nusselt  numbers 
with  Kays'  (108)  forced  convection  solution  showed  that  at  values  of 
z"*"  < 0.008  (corresponding  to  a z/D  of  about  18)  there  is  good  agree- 
ment. It  was  therefore  concluded  that,  in  this  inlet  region  of 
developing  flow,  the  effects  of  free  convection  are  largely  suppressed. 
At  values  of  z"*"  > 0.008  the  Nu  increased  sharply  above  that  given  by 
Kays,  which  is  indicative  of  significant  free  convection  effects  in 
this  region.  Maitra  and  Raju  (134)  wrote  the  momentum  and  energy 
equations  for  an  annulus  in  terms  of  Kelvin  functions,  then  approxi- 
mated these  by  polynomial  expressions . Experiments  were  also 

3 

performed  with  water  in  a brass-glass  annulus  over  2 x 10  < Ra  < 3 

X 10^  (based  on  D,  ) and  Re  from  200-1200.  The  experimental  Nusselt 
n 

numbers  were  about  45%  greater  than  theoretical  at  z/D.  >10.  In 

h 

the  entrance  region,  experiments  showed  that  fully  developed  flow 

+ -3 

was  established  at  z = 1.5  x 10  . In  the  same  region,  at  high 

Rayleigh  numbers,  the  experimental  Nusselt  number  was  greater  than 
the  theoretical  Nu  by  a factor  of  approximately  five.  This  runs 
counter  to  the  conclusions  of  Biggs  and  Stachiewicz  (133)  discussed 


above . 
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b . Combined  convection  In  horizontal  tubes 
In  a horizontal  tube  with  appreciable  heating  (of  a low- 
viscosity  fluid)  the  warmer  fluid  near  the  wall  is  less  dense  than 
the  rest  of  the  fluid  and  an  upward  motion  near  the  tube  wall  is 
imposed  on  the  forced  axial  motion.  The  fluid  far  from  the  wall  is 
cooler  and  denser,  so  it  has  a downward  component  of  motion.  For 
the  vertically  oriented  tube,  this  results  only  in  a distortion  of 
the  forced-convection  axial  velocity  profile.  But  for  a horizontal 
tube  a secondary  motion  is  established  which  is  symmetrical  about  a 
vertical  plane  passing  through  the  tube  axis  (Figure  26) . The 
resultant  flow  is  helical,  and  additional  mixing  takes  place  within 
the  fluid.  Heat  transfer  is  enhanced  thereby. 

Ede  (135)  was  probably  the  first  to  gather  experimental  data 
from  many  sources,  including  his  own  work,  for  combined  convection 
of  water  and  also  air,  in  tubes  of  various  diameters.  The  water 
data  is  fairly  consistent,  but  large  errors  are  present  in  the  air 
data  so  that  no  meaningful  results  were  obtained.  Two  effects  were 
noted:  a consistent  relation  between  the  Nusselt  and  Grashof 
numbers,  and  a substantial  variation  in  the  temperature  over  the  tube 
periphery. 

An  early  attempt  at  a theoretical  analysis  of  combined  hori-  ; 

zontal  convection  was  made  by  Morton  (136) , using  a perturbation 

method  with  Nu  • Gr/Re  as  the  expansion  parameter.  Truncation  error  i 

was  large,  ca.  10%  by  Re  • Ra  - 3000,  and  only  low  heating  rates  were  i 

'f 

assumed.  Even  so,  the  velocity  and  temperature  distributions,  and  ' 

I 


jL. 


Figure  26.  Temperature  and  velocity  profiles  (theoretical). 
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the  Nusselt  number,  were  found  to  be  functions  of  the  product  of  the 

Reynolds  number  and  the  Rayleigh  number,  for  a given  Prandtl  number. 

The  inadequacy  of  Che  pure  convection  solution  for  low  laminar  flow 

rates,  even  at  a small  heating  rate,  was  demonstrated;  distortion  of 

velocity  and  temperature  profiles  was  also  observed. 

The  above-noted  aspects  of  combined  convection  were  extensively 

studied  by  Mori  e^  a^.  (137)  in  an  experimental  program  that  included 

laminar  and  turbulent  flows.  A 35.6  mm  I.D.  tube  was  used,  with  air 

as  the  fluid;  its  properties  were  taken  as  variable  in  calculating 

Re  and  Ra.  A Rayleigh  number  of  ca.  40  resulted  from  a heating  rate 

on  the  order  of  0.1®C/(cm  tube  length),  which  is  much  less  than  that 

measured  during  the  solar  collector  experiments  here.  Significant 

distortions  in  the  temperature  and  velocity  profiles  were  noted, 

with  the  temperature  profile  affected  to  a greater  extent.  Turbu- 

4 

lent  flow  at  Re  > 10  overwhelmed  the  Incipient  secondary  flow,  so 

that  the  profiles  and  heat  transfer  were  essentially  unaffected  by 

buoyancy.  In  laminar  flow,  enhanced  heat  transfer  from  free  convec- 

3 

tion  appeared  at  about  Re  • Ra  = 10  , and  increased  substantially 
thereafter.  An  empirical  correlation  formula  for  calculating  the 
combined  convection  Nusselt  number  was  determined. 

Mori  and  Futagami  (138)  continued  this  study,  on  a theoretical 
basis.  A boundary  layer  approximation  was  applied,  yielding  axial 
temperature  and  velocity  profiles  that  are  distorted  by  buoyancy 
forces  but  deviate  somewhat  from  the  experimental  values.  However, 
the  approximation  was  considered  valid  as  Re  • Ra  becomes  very  large; 
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the  calculated  Nusselt  numbers  agree  quite  well  with  the  experimental  I 

data  reported  by  Mori  et  ^1.  (137) . The  calculated  equation  for  the  | 

I 

enhanced  Nusselt  number  is  | 

Nu  = Nug(0.1634) (Re  • Ra)°‘^  , for  Pr  =0.72  . 

The  perturbation  analysis  of  Paris  and  Viskanta  (119)  gave  very 
similar  Nusselt  numbers;  however,  the  analytical  transverse  profiles 
of  the  axial  velocity  and  temperature  components  were  flatter  than 
the  experimental  profiles  of  Mori  eit  (137) . But  the  Nusselt 
numbers  from  experiments  with  water  in  small  (2.5  cm  dia.)  tubes 
obtained  by  Newell  and  Bergles  (120)  are  considerably  lower  than 
those  of  Mori  ^ al.  (137)  for  air. 

The  typical  downward  distortion  and  flattening  of  the  trans- 
verse profiles  as  a result  of  free  convection  were  shown  to  take 
place  in  similar  fashion  in  non-circular  channels.  Cheng  and 
Hwang  (139) , using  a method  of  successive  over-relaxation  of  the 
governing  finite  difference  equations  (this  approach  diverges  at 
low  Re  • Ra) , analyzed  rectangular  and  square  channels.  The  velocity 
and  temperature  profiles  evinced  the  expected  distortion.  Also, 

Hong  and  Bergles  (140)  in  their  analysis  found  that  temperature- 
dependent  viscosity  (i.e.,  Pr)  can  exert  a strong  influence  on 
theoretical  heat  transfer  results. 

The  problem  of  combined  free  and  forced  convection  and 
entrance  effects  (i.e.,  thermally  and  hydrodynamically  developing 
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flow)  has  received  some  attention  in  the  last  decade  or  so.  All 
the  studies  unfortunately  assume  (theoretical) , or  have  (experi- 
mental) , fully  developed  velocity  profiles,  so  that  only  thermal 
development  is  treated.  These  studies  were  severely  limited  in 
scope  and  by  their  assumptions,  and  disagree  widely.  Therefore, 
the  study  results  cannot  be  utilized  here  directly.  Still,  some 
trends  are  evident  and  must  be  taken  into  account,  at  least  quali- 
tatively. 

Shannon  and  Depew  (141,142)  performed  experiments  with  both 

water  and  ethylene  glycol,  in  small  (<  1 cm  dia.)  tubes,  and  at 

3 5 

Grashof  numbers  (based  on  D ) to  2.5  x 10  . Their  heat  transfer 
results,  based  on  locally-measured  fluid  properties,  indicate 
that  the  Nu&selt  number  begins  to  depart  from  the  pure  convection 
case  quite  early.  The  Nu  reaches  a minimum  as  flow  progresses  yet 
remains  above  the  pure  convection  value,  and  then  deviates  sharply 
from  the  latter  (which  becomes  asymptotic)  by  reversing  its  previous 
trend  and  now  increases  with  the  axial  dimension.  The  effect  is 
more  pronounced  at  higher  Gr.  The  Nusselt  number  test  data  scatter 
was  + 25%,  and  the  results  from  water  were  below  the  ethylene  glycol 
values;  backflow  was  postulated  to  account  for  the  difference. 

Very  similar  experiments  with  water  and/or  ethylene  glycol 
have  been  performed  by  Morcos  and  Bergles  (143) , Bergles  and  Simonds 
(144) , and  Hong  £t  a^.  (145) . A large  variation  in  peripheral  wall 
temperature,  to  50®C  between  the  top  and  bottom  of  a glass  tube  (145), 
was  observed  in  these  experiments.  The  results  exhibit  a substantial 
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amount  of  scatter,  but  indicate  that  very  near  the  entrance  the 
Nusselt  number  is  the  same  as  in  pure  forced  convection.  Farther 
along,  the  Nusselt  number  curve  (plotted  versus  axial  distance) 
branches  out  and  approaches  a constant  value  that  is  a function  of 
the  Rayleigh  number.  Development  of  the  temperature  profile  is 
rapid,  requiring  a decreasing  length  of  tube  as  the  Rayleigh  number 
increases . 

Finite  difference  analyses,  which  use  a high-Prandtl-number- 
fluid  assumption  (i.e.,  primary  axial  flow  is  unaffected  by  the 
secondary  flow)  so  that  the  inertia  terms  in  the  momentum  equation 
can  be  neglected,  give  results  in  good  agreement  with  the  experimental 
data  for  water  (145).  Cheng  ^ (92),  and  Cheng  and  Ou  (146) 

compared  their  theoretical  results  with  the  experimental  data 
reported  by  others  and  found  reasonable  agreement,  considering  the 
scatter  in  the  experimental  data.  Cheng  ^ (92)  also  noted  a 

dip,  or  minimum,  in  the  theoretical  Nusselt  number  as  flow  progresses, 
before  the  asymptotic  or  fully  developed  value  is  attained. 

Of  greater  relevance  to  the  present  analysis  are  the  experi- 
mental results  with  air  as  the  heat  transfer  medium.  McComas  and 
Eckert  (147)  investigated  the  range  of  100  < Re  < 900  and  Gr  < 1000 

3 

(based  on  D ) , but  in  a small  (ca.  1 cm  dia.)  tube.  Gas  properties 
were  evaluated  at  the  local  axial  temperature.  At  higher  Grashof 
numbers,  the  local  (axial)  Nusselt  numbers  at  first  followed  the 
forced  convection  results  as  flow  progressed,  then  deviated  and 
became  constant  at  a higher  level  - which  was  a function  of  the  Gr. 


A slight  dip  was  observed  in  the  Nu  versus  axial  position  curve, 
but  only  at  low  Re;  the  peripheral  temperature  variation  was 
practically  nil. 
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A wider  (25.4  mm  I.D.)  tube  was  used  by  Hussain  and  McComas 

(148).  A peripheral  temperature  variation  of  up  to  7°C  for  a wall- 

to-bulk  temperature  difference  of  28°C  occurred  in  the  upper  Grashof 

4 3 

number  range  (ca.  5 x 10  , based  on  D ) . Heat  transfer  results  were 
somewhat  puzzling:  at  Re  < 1200  the  asymptotic  Nusselt  number  was 
lower  than  that  for  pure  forced  convection.  This  may  have  been 
caused  by  exit  end  effects.  For  Re  > 1200,  enhanced  heat  transfer 
was  observed,  and  agreement  with  the  data  of  Mori  ^ al . (137), 
who  studied  a narrower  Rayleigh  number  range,  is  quite  good.  But 
the  empirical  correlation  of  Mori  (137),  discussed  earlier, 

appears  to  overpredict  the  Nu  at  higher  Ra. 

c.  Combined  convection  in  inclined  tubes 

Only  a few  studies  on  inclined  tubes  are  extant;  flow  develop- 
ment with  combined  convection  has  not  been  studied  in  inclined  tubes. 
Iqbal  and  Stachiewicz  (118,121)  calculated  the  temperature  and  veloc- 
ity fields  by  a perturbation  method.  This  method  gives  unrealistically 
high  estimates  of  the  Nusselt  number  at  even  moderate  values  of 

3 

Re  • Ra  (ca.  10  ) (145).  The  results  of  that  study  (118),  to 
3 

Re  ' Ra  < 3 X 10  , not  only  give  erroneous  values  of  Nu,  but  also 
show  the  occurrence  of  a maximum  in  the  Nu  at  an  inclination  angle 
between  20-60° . All  other  studies  indicate  that  the  combined  heat 
transfer  is  a maximum  at  an  angle  of  0°  (i.e.,  horizontal  tube)  where 
interaction  between  primary  and  secondary  flow  is  maximized. 
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Recently  Hwang  and  Cheng  (149)  approached  the  problem  of 
combined  convection  in  horizontal  tubes  by  a combination  of  line 
iterative  and  boundary  vorticity  methods.  The  calculated  trans- 
verse profiles  of  the  axial  temperature  and  velocity  match  the 
experimental  data  of  Mori  e_c  (137)  rather  well,  to  a considerably 
greater  extent  than  the  theoretical  results  of  Mori  and  Futagaml 
(138) , This  numerical  method  has  been  applied  to  inclined  tubes  by 
Cheng  and  Hong  (150) , with  solutions  obtained  by  a combined  boundary 
vorticity/iterative  relaxation  method  of  the  governing  equations, 
which  are  written  in  finite  difference  form.  The  greatest  distor- 
tion of  temperature  and  velocity  profiles  is  shown  to  occur  in  a 
horizontal  tube,  and  it  decreases  with  inclination  angle.  Heat 
transfer  results  show  a similar  trend.  For  a gas  (Pr  = 0.75)  and 
a circular  tube  at  an  inclination  of  45°,  the  Nusselt  number  was 
found  to  depend  on  the  Rayleigh  number  and  the  Re  • Ra  product, 

4 

to  Ra  $ ca.  450  (for  Re  • Ra  < 5 x 10  ) . At  larger  Rayleigh 
numbers,  the  Nu  becomes  independent  of  Re  and  is  a function  of  Ra 
only.  This  numerical  solution  checks  exactly  with  the  known  numer- 
ical results  for  the  horizontal  and  the  vertical  tube. 

The  experiments  of  Sabbagh  e^  £l.  (151)  with  a gas  in  an 
inclined  tube  were  an  attempt  to  corroborate  the  theoretical  results 
discussed  above.  Although  qualitatively  the  experimental  data  are 
similar  to  the  theoretical  results,  they  are  seriously  flawed  since 
even  the  Nusselt  number  obtained  at  0°  inclination  is  much  greater 
than  measured  by  Mori  (137):  7.85  versus  ca.  5.6.  The 
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measured  temperature  and  velocity  profiles  are  also  very  different 
from  those  calculated  by  Cheng  and  Hong  (150) . 

Better  results  were  obtained  by  Futagami  and  Abe  (152)  in 
a simultaneous  experimental/ theoretical  study.  They  were  able  to 
correlate  results  at  various  angles  (a)  of  inclination  using  as  a 
correction  factor  the  cosine  of  a.  For  example,  the  temperature 
and  velocity  profiles,  for  inclination  angles  of  0°  and  45°  and  at 
Re  • Ra  cos  a = 8 x 10  , are  nearly  identical.  A comparison  was 
made  with  the  analytical  heat  transfer  results  of  Hwang  and  Cheng 
(149)  discussed  earlier.  Even  chough  Chose  results  were  for  a 
horizontal  orientation,  reasonable  agreement  was  found  with  the 
values  obtained  by  Futagami  and  Abe  (152) , when  plotted  against 
Re  • Ra  cos  a as  a parameter. 

Inclined  rectangular  channels  of  aspect  ratios  =0.5,  1.0 
(square  cross  section),  and  2.0  were  examined  by  Ou  £t  (153) 

in  a numerical  analysis.  A Prandtl  number  of  5 (i.e.,  water)  was 
used  throughout.  Results  for  various  inclinations  were  correlated 
through  the  use  of  a trigonometric  function  in  the  Ra  and  Re 
calculations.  The  results  clearly  indicate  that  the  effect  of 
cross-sectional  shape,  over  the  range  studied,  was  negligibly 
small.  This  justifies  the  use  here  of  circular-tube  results  for 
the  semicircular  receiver-halves. 

In  summary,  given  the  paucity  and  variability  in  the  above 
studies,  a certain  subjective  judgmental  factor  enters  into  the 
choice  of  the  results  to  be  used  herein.  It  was  decided  to  base 
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the  heat  transfer  coefficient  from  combined  convection  on  the 
theoretical  work  of  Cheng  and  Hong  (150) . Because  that  work  does 
not  adequately  cover  the  full  Re  • Ra  range  encountered  here,  Che 
studies  of  Hwang  and  Cheng  (149)  in  the  low  Re  • Ra  range,  and 
Hussain  and  McCoraas  (148)  in  the  intermediate  Re  • Ra  range  are 
also  used.  Since  the  last  two  were  derived  for  horizontal  tubes, 

Che  correction  cos  a of  Futagami  and  Abe  (152)  is  incorporated 
into  those  results  (Figure  27) . The  complete  set  of  expressions 
defining  enhanced  heat  transfer  from  natural  convection  is: 

(Nu^)^5o/(Nu^)go  = 0.27(Re  • Ra)®'^^ 

for  Re  • Ra  < 2.75  x 10^  , (after  (149)) 

(Nu^)^5<,/(NUf)oo  = 2.0 

for  Re  • Ra  > 2,75  x 10^  , (after  (148)) 


and 


(Nu^)^5<,/(Nu^)qo  = 0.31(Ra)°-^^ 


for  Ra  i 1000  , 


(after  (150)) 


where 


c 45 


is  the  Nusselt  number  in  combined  convection,  at  an 


inclination  angle  of  45°,  and  (Nu^)qo  is  the  Nusselt  number  in  pure 


Figure  27.  Nusselt  number  ratio. 


forced  convection  in  a horizontal  tube.  These  expressions  were 
incorporated  into  the  mathematical  model  in  such  a way  that  there 
was  no  overlap  among  them  (Figure  27) . 
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A number  of  hand  calculations  were  performed  in  order  to 
support  the  choice  of  the  above  expressions.  Data  from  experiments 
with  the  solar  collector  were  used  as  input.  It  was  determined  that: 


(1)  The  cross-sectional  shape  of  the  receiver  can  be  ignored 

(i.e.,  the  Nu  ratio  above  is  the  same,  for  circular  and  semicircular 

tubes) . Two  calculations  were  made  by  integrating  temperature  and 

velocity  profiles,  taken  from  Cheng  and  Hong  (150),  in  the  definition 

of  the  bulk  fluid  temperature,  T,  . The  difference  in  the  calculated 

b 

bulk  temperatures  for  a semicircular  and  circular  cross  section  was 
only  ca.  3%  at  the  Re  and  Ra  values  typically  measured. 


(2)  The  difference  between  calculated  and  measured  f - T , 

w g’ 

for  various  flow  rates  and  at  all  axial  positions  was  quite  small. 

The  ratio  (T  - T.  ) , /(T  - T ^ ..  i • ) spanned  the 

w b calc.  w g meas.  at  centerline 

range  0. 8-1.0,  with  most  values  clustered  near  1.0.  The  experi- 
mental T^  was  found  by  using  the  temperature  distributions  in  the 
Pyrex  and  the  graphite  as  discussed  earlier;  i.e. , it  can  be  shown 
that  T =0.7T  +0.3  T„  . The  good  agreement  between 

calculated  and  measured  results  permits  several  conclusions  to  be 
drawn: 


(a)  The  Nusselt  numbers  are  reasonably  accurate. 


‘ , ir' 
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(b)  The  measured  gas  temperature,  T^,  which  is  measured  at  a 
point  in  the  vertical  center  of  the  upper  half  of  the  semicircular 
receiver,  closely  represents  the  bulk,  or  mean  mixed,  fluid  tempera- 
ture T^.  Such  would  not  normally  be  the  case,  except  when  the 
temperature  and  velocity  profiles  are  completely  flat;  the  ratio  of 

(Tw  - Tb)/(\  - Tg  centerline^  developed  profiles  is 

not  negligible.  Here,  the  downward  distortion  in  the  velocity  and 
temperature  profiles  is  sufficient  to  nearly  offset  the  large  varia- 
tions in  temperature  and  velocity  encountered  in  the  transverse 
plane  of  the  receiver.  This  is  so  because  the  measured  centerline 
fluid  temperature  is  actually  on  the  shoulder  of  the  profile,  rather 


than  at  its  minimum  (Figure  26) . 


Strictly  speaking,  these  conclusions  apply  only  to  the  upper 

half  of  the  semicircular  receiver,  where  all  temperature  measurements 

were  made  (because  of  the  concentrated  solar  flux  in  the  lower 

receiver  half,  no  thermocouples  were  placed  there).  It  will  be 

assumed  here  that  the  heat  transfer  processes  within  the  lower  half 

are  quantitatively,  if  not  qualitatively,  very  similar.  Therefore, 

the  measurements  made  consistently  represent  the  true  thermal  energy 

level  of  the  heat  transfer  fluid,  i.e.,T  ^ ,,  =T,  . This 

g at  centerline  b 

eases  the  computation  enormously. 
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7.  Energy  Balance 


a.  Introduction 


Since  the  gas  used  (Ar)  is  transparent  to  solar  and  thermal 
radiation,  it  heats  by  convection  only.  The  driving  potential  for 
its  heating  is  the  difference  between  the  gas  and  (averaged)  wall 
temperatures;  this  difference  can  be  calculated  if  the  Nusselt 
number  is  known.  The  radial  losses,  consisting  of  conduction 
through  the  insulation  and  radiation  from  the  aperture,  are  purely 
a function  of  the  wall  temperature  if  an  external  (i.e.,  skin) 
temperature  is  known  or  posited.  The  axial  conductive  loss, 
mainly  along  the  absorber,  can  be  determined  from  the  rate  of 
axial  temperature  change. 

b.  Energy  balance  equations 

An  energy  balance  is  now  made  on  a small  (=  differential) 
control  volume  (see  Figure  28) : 


% = * ‘ ’is'  \ - ‘»C2  ’ 


or,  in  finite  difference  form 


p b m b m-1  S w m-1  w m-1  e 


k,  1 [(T  ) , - T ] - (k^  + k^)[2(T,_)  - (T,  ) , - (T,  ) _ ] 

insul.  w m-1  e P G b m b m-1  b nrt-1 


157 


where  the  first  term  on  the  R.H.S.  represents  the  available  concen- 
trated solar  flux  at  the  Pyrex  window;  the  second  is  the  radiative 
loss  from  the  aperture,  both  from  the  Pyrex  and  the  aperture  sides; 
the  third  is  the  loss  through  the  insulation;  the  last  term  repre- 
sents the  axial  conduction  loss,  by  the  Pyrex  and  the  absorber.  For 
clarity,  the  appropriate  view  factors,  areas,  and  radiative  properties 
are  not  shown  in  the  equation.  The  problem  is  a two-dimensional  one, 
with  the  first  two  loss  terms  essentially  consisting  of  (one- 
dimensional) radial  losses.  The  second,  axial,  dimension  enters 
through  the  axial  conduction  loss  term. 

As  described  earlier,  T is  now  the  average  peripheral  wall 

w 

temperature  of  Pyrex  + graphite.  T is  obtained  (except  at  the  first 

w 

step;  see  Appendix  I)  by  solving  the  governing  equation  for  the 
convective  heat  transfer  within  the  collector: 

q„  = h(f  - T.  ) . , 

w b m-1 

where  h = Nu  • k/Dj^.  All  the  parameters  except  are  known  or  can 
be  calculated. 

There  are  two  unknowns,  (T,  ) , which  is  being  solved  for, 

b m 

given  the  known  or  calculated  (T.  ) ,,  and  (T.  ) there  is  only 

b m-1  b m+1 ; 

one  equation.  A second  equation  must  now  be  written  for  an  additional 
step,  i.e. , 


- 
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* - «b>.!  • ’S  - <“<V«  ^ »1«V.  * 

- "insul.KV.  - + V‘2«b).  - 

A typical  node  is  illustrated  on  Figure  28.  Two  additional  approxi- 
mations were  made  in  the  equations  (nodes),  to  make  them  manageable: 

(1)  The  solution  would  not  converge  when  the  wall  tempera- 
tures were  used  in  computing  the  axial  loss.  Thus,  it  was  necessary 

2-2  2 2 

to  replace  these  with  the  gas  temperatures,  i.e.,  (3  T /3z  ) = (3  T, /3z  ), 

W D 

or 


[2(T  ) - (T  ) , - (T  )_] 

w m w m-1  w m+1 


[2(T.  - (T,  ) 

b m hi 


- (T  ) 1 

-1  '•  K-'m+l  ■ 


An  examination  of  the  experimental  data  (Figure  25)  shows  that  this 
is  a reasonable  approximation. 

(2)  The  axial  conduction  loss  term  is  the  same  in  both  equa- 
tions used  in  the  solution,  so  as  not  to  introduce  new  unknowns,  i.e., 
now 


(3^T  /3z^)  = (3V)/3z^)  = constant  , 

W 0 

over  the  two  steps.  Again,  examination  of  the  data  (Figure  25) 
reveals  that  the  slope  of  the  temperature  profile  changes  slowly, 
so  the  approximation  is  valid. 


t 
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Two  further  assumptions  common  to  this  type  of  problem  are 


made : 


(1)  Friction  can  be  ignored  because  of  smooth  surfaces  and 


low  velocities. 

F 

S 

^ (2)  Viscous  dissipation  and  compressibility  effects  can  be 


neglected  (154,155),  again  because  of  the  low  velocities. 


I * 

i 

I 


The  two  equations  considered  can  now  be  readily  solved  by 
an  iterative  approach.  A simultaneous  solution  is  precluded  by 
the  difficult  nonlinear  terms. 

c.  Computational  method 

A computer  program  was  written  in  the  FORTRAN  IV  language, 
and  the  program  was  executed  on  a PDP-9  computer.  The  program 
exceeded  the  computer's  memory  storage  capacity  and  had  to  be 
subdivided  into  three  parts.  The  program  was  executed  rapidly; 
the  relatively  slow  mechanical  printer  was  the  limiting  factor  on 
run  time. 

The  program  performs  an  energy  balance  on  a 1 cm-long  control 
volume  of  the  receiver  and  calculates  the  axial  bulk  temperature 
profile  of  the  receiver  over  the  32-to-123  cm  length.  The  two 
equations  mentioned  are  solved  iteratively  at  each  step,  until  the 
solution  converges  on  a predetermined  value.  The  bulk  temperature 
profile  itself  is  calculated  by  a forward-one  numerical  integration 
of  the  equation,  with  the  solution  marching  on  from  the  initial 


* 

3 


I 


conditions  at  the  32-cm  axial  position.  The  main  features  of  this 
program  are: 

(1)  Use  of  actual  gas  properties  at  the  T,  of  the  control 

b 

element  under  consideration.  The  argon  properties  are  tabulated  by 

Hilsenrath  e^  ^1.  (156) , with  some  properties  also  listed  by 

Raznjevid  (157).  Properties  are  evaluated  at  T,  to  yield  a minimal 

b 

error  in  determining  the  heat  transfer  (37). 

(2)  Inclusion  of  enhanced  heat  transfer  due  to  both  develop- 
ing flow  and  buoyancy.  Each  effect  is  assumed  to  act  independently 
of  the  other  (i.e.,  no  coupling  present). 

(3)  Consideration  of  axial  heat  loss  at  the  inlet  by  conduc- 
tion in  the  graphite  and  glass. 

In  the  actual  computation,  the  program  (see  Figure  29  and  Appendix  I) 

(1)  Determines  the  average  direct  flux,  from  the  two  measure- 
ments described  earlier. 

(2)  Calculates  the  argon  flow  rate. 

(3)  Calculates  gas  properties  and  Nu's  from  the  T of  the 

b 

previous  step. 

(4)  Calculates  thermal  (transverse  only)  losses. 


(5)  Calculates  T 
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(6)  Repeats  (4)  and  (5)  for  an  additional  step. 

(7)  Computes  an  axial  conductive  loss  from  the  temperature 
gradient  difference  between  the  two  steps. 

(8)  Iterates  (5) , (6) , and  (7)  until  the  solution  converges 
within  0.01“C.  The  axial  loss  from  (7)  is  now  included  in  calcu- 
lating the  T,  in  (5) . 

D 

(9)  Goes  on  to  the  next  step  and  begins  at  (3)  again. 

d.  Stability  and  convergence 

It  is  possible  for  some  small  error  (e.g.,  round-off  error) 
to  increase  as  the  computation  progresses.  The  instability  can 
cause  the  error  terms  to  completely  overshadow  the  desired  solution, 
making  it  worthless. 

Karplus  (158)  used  an  electrical  circuit  analog  of  finite 
difference  equations  to  determine  the  criteria  for  stability.  A 
typical  node  in  the  solution  was  considered;  if  the  node  is  stable 
(i.e.,  its  analog  serves  to  dissipate  rather  than  supply  energy  to 
other  nodes,  in  response  to  an  excitation)  then  the  entire  network 
must  be  stable.  When  the  finite  difference  equation  is  put  in  a 
form  with  all  the  terms  on  one  side,  the  stability  criteria  become: 

(1)  If  all  the  coefficients  of  the  finite  difference  expres- 
sions are  positive,  then  the  equation  is  always  stable. 
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(2)  If  some  coefficients  are  negative,  the  Z coefficients 
must  be  < 0 as  a sufficient  condition  for  stability. 

Here,  the  terms  in  the  energy  balance  equation  can  be  rear- 
ranged, and  it  is  apparent  that  all  the  coefficients  of  the  finite 
difference  expressions  are  positive.  The  equation  and  its  solution 
are  therefore  unconditionally  stable.  The  control  element  size 
(i.e.,  the  number  of  steps,  or  subdivisions,  in  the  computation) 
is  thus  detemined  purely  from  considerations  of  accuracy  in  a 
trade-off  against  computer  time  (159) . 

The  normal  control  element  size  selected  was  1 cm,  with  a 
convergence  criterion  of  0.01°C  for  the  iterative  solution  of  the 
axial  conduction  loss.  A fair  number  of  calculations  were  repeated, 
with  an  element  size  of  0.01  cm  and  convergence  on  0.001°C.  Execu- 
tion time  was  greatly  lengthened,  and  the  temperatures  were  only 
slightly  higher  than  those  normally  calculated.  The  difference 
was  between  1-7°C,  and  provided  no  better  agreement  with  the 
experimental  data,  so  the  less  stringent  approach  was  adopted. 

Lax's  (160)  equivalence  theorem  states  that,  given  a properly 
posed  initial  value  problem,  and  a finite  difference  approximation 
to  it  that  satisfies  the  consistency  condition,  stability  is  a 
necessary  and  sufficient  condition  for  a convergent  approximation. 

A problem  is  taken  as  properly  posed  if  a family  of  genuine  solu- 
tions exists,  and  if  the  solutions  depend  uniquely  and  continuously 
on  the  initial  data  (161).  For  consistency,  the  truncation  error  ^ 0 
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as  the  element  size  ->■  0.  Since  both  these  conditions  appear  to  be  i 

satisfied  here,  and  no  problems  with  convergence  beyond  those 
described  earlier  were  encountered,  the  solution  is  assumed 
convergent . 

D.  Results  of  Computation 

1 . Graphite  Absorber 

The  results  of  several  typical  computations,  together  with  the 
actual  experimental  measurements  in  parentheses,  are  shown  in  Appen- 
dix II  and  Figures  30-32.  Agreement  between  the  experimental  and 
calculated  temperatures  is  quite  good  here,  but  becomes  rather  poor 
at  the  higher  flow  rates  as  previously  discussed  in  Chapter  III, 

Section  E (see  also  Figure  33).  The  influence  of  flow  rate,  flux, 
and  inlet  temperature  on  the  temperature  profiles  is  apparent.  The 
efficiency  of  heat  utilization  is,  as  one  might  expect,  high  at 
first  (the  SUF,  a measure  of  thermal  efficiency,  is  above  0.6  for 
two  of  the  runs  illustrated  in  Appendix  II) , but  decreases  rapidly  as 
the  gas  heats  and  losses  mount.  By  the  end  of  the  test  section  it 
is  only  0.04-0.23,  indicating  that  little  additional  heating  can 
take  place  farther  downstream.  The  average  wall  temperature  profiles 
parallel  the  graphite  absorber  temperature  profiles  (measured  graphite 
temperatures  shown  in  parentheses)  but  are  somewhat  lower  because  of 
the  inclusion  of  the  glass  receiver  temperature  in  the  average. 

There  is  some  error  present  in  the  flux  measurements/calculations 
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(compare  FXl  with  FX2  in  Appendix  II)  and  possibly  some  error  in 
an  assumed  5%  loss  because  of  misalingment  between  the  reflected 
flux  and  the  receiver  aperture.  All  the  results  are  summarized 
and  compared  in  Table  9. 

A number  of  calculations  were  also  repeated,  when  there  was 
no  receiver  inlet  shadowing,  starting  right  at  the  first  thermo- 
couple position  2 cm  along  the  aperture.  In  general,  a larger 
error  is  present  in  these  calculated  results  since  the  inlet  end 
effect  becomes  more  severe  and  some  of  the  assumptions  become  less 
tenable . 

Furthermore,  the  use  of  the  heating  tape  was  found  to 
provide  little  new  information  on  the  performance  of  the  collector. 

At  low  flow  rates,  the  exit  temperature  of  the  gas  already  approaches 
its  as)nnptocic  limit  (see,  for  example.  Figures  30  and  31);  preheating 
the  incoming  gas  thus  serves  only  to  move  the  onset  of  near-asymptotic 
conditions  closer  to  the  receiver  inlet.  The  agreement  between  calcu- 
lated and  experimental  temperatures  is  also  somewhat  poorer,  as 
asymptotic  conditions  are  approached  (Table  9 and  Figure  30) . 
Sufficient  variability  in  inlet  temperatures  was  provided  by  the 
changes  in  flow  rate,  ambient  temperature,  solar  flux,  and  aperture 
shadowing;  preheating  was  therefore  a redundant  feature  for  this 
collector . 

2.  Selective  Absorber 

The  standard  program  was  modified  to  account  for  the  different 
dimensions  and  properties  of  the  selective  copper  absorber.  Agreement 
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Table  9.  Comparison  of  measured  and  calculated  temperatures. 


Perp.  direct 


Date 

Time, 

p .m. 

solar  flux, 

W m"^ 

Flow  rate, 

1 • 

1 min 

AT,  ’C  (i. 

Measured 

,e. , T17  - T14) 

Calculated 

15.x. 76 

1:30 

754 

44.7 

135 

147 

2:45 

717 

89.5 

55 

94 

23.x. 76 

1:30^ 

787 

32.8 

172 

184 

2:30 

750 

62.3 

84 

131 

14. XI. 76 

1.00 

788 

98.1 

56 

105 

19. XI. 76 

12:30 

647 

11.9 

201 

189^ 

1:50 

566 

26.6 

158 

158 

14. XII. 76 

12:45 

766 

40.7 

158 

187 

11.11.77 

1:15 

772 

95.2 

63 

103 

2. III. 77 

1:30 

872 

65.2 

110 

148 

9. III. 

1:15^ 

864 

62.6 

98 

141 

11. III. 

1:00 

776 

52.1 

111 

134 

17. III. 

1:00^^ 

852 

38.7 

151 

159 

19. III. 

1:30 

872 

32.3 

155 

161 

23. III. 

1:45® 

872 

29.3 

155 

142 

29. III. 

1:15 

860 

20.9 

139 

121 

30. III. 

12:25 

836 

40.6 

128 

141 

l.IV. 

12:45 

867 

45.1 

134 

151 

12. IV. 

1:45 

721 

31.9 

48 

34 

15. IV. 

2:15 

876 

40.5 

127 

147 

) 


> 


\ > 

I %• 
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Table  9.  Continued. 


Time , 

Perp. 

solar 

direct 

flux. 

Flow 

rate,  AT,  °C 

(i.e.,  T17  - 14) 

Date 

p .m. 

W 

1 min 

^ Measured  Calculated 

17. IV. 77 

1:45 

Very  high  preheat  temperature;  program  stalls. 

8.V. 

1:35^ 

782 

33.0 

96 

99 

14. V. 

1:45 

804 

42.6 

93 

105 

2:45 

793 

42.2 

51 

56 

26. V.^ 

1:30 

721 

20.7 

132 

98 

2:30 

764 

25.8 

124 

110 

3:30^ 

780 

31.6 

111 

117 

3. VI. 

1:30 

795 

36.8 

120 

130 

2:30 

813 

40.5 

117 

138 

3:30 

822 

44.0 

112 

143 

4. VI. 

1:30 

758 

26.0 

119 

108 

2:30 

775 

31.5 

116 

121 

3:30 

773 

35.9 

111 

128 

^See  Figure 

34 

^Calculated 

with  a 

0.1  cm 

control 

volume 

and  0.001  °C 

axial  loss 

convergence.  See  Figure  30 


'See  Figure  33 


See  Figure  32 


'See  Figure  31 


Selective  absorber  now,  and  hereafter 
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between  the  experimental  and  calculated  temperatures  is  not  quite 
as  good  as  for  the  graphite  absorber.  The  buckling  and  rippling 
of  the  selective  absorber  is  the  probable  cause  of  the  larger 
discrepancy,  as  this  undoubtedly  disturbed  the  flow. 

Comparison  of  an  experiment,  where  the  measured  and  theoret- 
ical temperatures  are  in  close  agreement,  with  the  graphite  absorber 
results  for  similar  conditions  (a  composite  of  two  graphite  results 
was  necessary,  to  span  the  temperature  range)  is  shown  in  Figure 
34.  It  is  apparent  that,  within  the  limits  of  experimental  and 
computational  error,  no  gain  results  from  the  selective  surface. 

This  speaks  well  for  the  already  high  efficiency  (i.e.,  low  apparent 
emittance)  of  the  cavity  receiver  design  analyzed  here.  The  higher 
end  loss  from  the  graphite  is  evident  (Figure  34);  its  emittance 
is  of  course  much  higher,  and  it  conducts  more  heat  because  of 
its  thickness. 


E.  Conclusions 


One  phase  in  the  development  sequence  of  a particular  solar 
thermal  collector  design  has  been  described.  The  difficulties  in  a 
realistic  modeling  effort  are  readily  apparent,  yet  do  not  compromise 
the  validity  of  the  results  obtained.  The  experience  gained  with 
this  model  can  prove  useful  for  the  examination  of  a larger  model, 
or  even  the  full-scale  device.  Sufficiently  turbulent  flow  probably 
cannot  be  attained  in  any  experimental  model  short  of  an  extremely 
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large  one,  so  that  some  buoyancy  and  developing  flow  effects  even 
In  a turbulent  flow  regime  would  still  be  present.  However,  their 
Influence  could  be  much  reduced  with  a substantially  larger  collector 
and  higher  flow  rates.  A closer  agreement  between  calculated  and 
measured  temperature  profiles  would  thus  occur  In  a larger  model. 

In  conclusion,  the  basic  concept  studied  here  Is  sound,  and 
reasonably  efficient  utilization  (see  SUf,  or  ASUF,  In  Appendix  II) 
of  solar  energy  at  high  fluid  temperatures  has  been  observed.  This 
utilization  Is  of  course  much  lower  than  that  estimated  for  the 
full-scale  design  (see  Chapter  V),  because  of  end  effects  and 
greater  radial  heat  losses.  Again,  a larger  model  should  evince 
higher  thermal  efficiencies,  approaching  chose  of  the  full-scale. 

An  indication  of  the  potential  power  production  for  the  particular 
device  described  here  Is  Illustrated  by  an  example:  at  Che  32.3 
1 min  ^ flow  rate  (Appendix  II) , energy  is  stored  in  the  gas  as  it 
flows  through  the  91-cm  long  test  section  at  a rate  equal  to  385  W. 
This  power  is  much  too  low  to  drive  any  commercially  available  gas 
turbine. 
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CHAPTER  V 


MATHEMATICAL  MODEL  OF  THE  FULL-SCALE  DESIGN 


A.  Approach 

The  good  agreement  between  experimental  and  calculated 
temperature  profiles  In  the  scale  model  permits  the  formulation 
of  a second  computer  model,  of  the  proposed  full-scale  device  (32), 
with  confidence.  This  mathematical  model  of  the  full-scale  collector 
Is  much  less  Involved  than  the  computer  model  of  the  experimental 
device. 

The  length  (700  m)  and  relatively  flat  axial  temperature 

profile  of  the  full-scale  collector  allow  both  end  effects  and 

axial  conduction  to  be  neglected.  The  highly  turbulent  flow  of 

4 

the  helium  (Re  = ca.  9 x 10  ) in  the  receiver  reduces  free  convec- 
tion effects  to  negligible  values.  Petukhov  ^ al.  (162)  have 
shown  that  at  this  Reynolds  number,  and  at  the  anticipated  maximum 
Grashof  number  in  the  receiver,  free-convection  effects  will  cause 
a deviation  in  the  forced  convection  Nusselt  number  of  less  than  1%. 
The  same  authors  also  show  that  flow  development  is  essentially 
complete  in  ca.  10  to  20  diameters,  so  both  these  factors  will  be 
ignored  here. 
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The  temperature  gradient  between  the  heat  transfer  surfaces 
(glass  + graphite)  and  the  gas  (He  in  the  full-scale)  can  be 
calculated  from  the  relation  given  by  Kays  (163) : 

Nu  ■ 0.022  Pr^‘^  Re^’^  . 

A constant  axial  heat  input  is  assumed  in  the  above;  however,  if  a 
constant  axial  wall  temperature  is  taken  Instead  as  the  more  accurate 
representation,  the  only  change  in  the  relation  occurs  in  the  coef- 
ficient, from  0.022  to  0.021.  At  any  rate,  the  calculated  maximum 
temperature  difference  is  only  ca.  8“C,  a trivial  quantity.  Thus 
the  receiver,  exclusive  of  the  insulation,  can  be  taken  as  isothermal 
radially.  The  computation  therefore  becomes  a straightforward  one- 
dimensional numerical  integration  of  the  energy  balance  on  a small 
control  volume,  and  yields  the  axial  temperature  profile  of  the 
flowing  gas. 

The  main  features  of  the  computation  are  as  follows: 

(1)  A control  volume  (element)  length  of  10  m. 

(2)  Consideration  of  radial  losses  only  (i.e.,  essentially 
1-D,  although  of  course  the  radiant  loss  is  mainly  2-D) , which 
Include  conduction  through  the  Insulation  and  radiation  by  the 
aperture,  both  by  Che  aperture  sides  and  the  glass  window. 
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(3)  Use  of  a variable  view  factor  from  the  aperture,  to 
account  for  examination  of  different  pipe  diameters  and  insulation 
thicknesses.  The  view  factor  is  computed  by  the  crossed-string 
method,  as  described  in  Chapter  IV. 

/ 

The  standard  configuration  of  the  full-scale  device  has  the 

following  dimensions  (from  (32))  and  properties  (which  are  identical, 

except  for  the  aperture  emittance,  with  the  material  properties  of 

the  scale  model):  a receiver  pipe  of  30  cm  inside  diameter,  glass 

pipe  and  graphite  absorber  thicknesses  equal  to  0.05  of  this  diameter, 

10  cm  thick  Flberfrax  insulation  with  temperature-dependent  thermal 

conductivity,  a 6 cm  wide  aperture  with  an  emittance  (both  glass 

and  aperture  sides)  of  0.75,  a mirror  array  width  of  4 m and 

reflectance  ~ 0.6,  and  a convective  cover  transmittance  of  0.85. 

The  more  important  parameters  are  Individually  varied  in  a computer 

program  written  in  FORTRAN  IV  (see  Appendix  III) . The  helium  inlet 

temperature  is  always  taken  as  349°C,  and  the  mass  flow  rate  is  a 

constant  88  kg  a ^ at  ca.  7 atm  pressure  (from  (32)).  A conserva- 

_2 

tlve  direct  perpendicular  solar  flux  of  800  U m , and  a 25°C 
ambient  temperature  (this  is  also  the  external  Insulation  temperature) 
are  assumed  throughout. 

B.  Results  of  Computation 

At  a given  (concentrated)  solar  flux  level,  the  only  signif- 
icant method  for  improving  collector  thermal  efficiency  involves 
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reduction  of  thermal  losses.  This  can  be  achieved  In  any  number  of 
ways;  the  most  Influential  parameters  for  the  radlatlon-cavlty 
collector  design  are  the  following: 


(1)  Increasing/decreasing  the  standard  6 cm  aperture  width 
results  In  approximately  a lO^C/Ccm  of  width  changed)  gas  outlet 
temperature  decrease/increase,  respectively  (see  Figure  35).  The 
aperture  ought  to  be  kept  as  small  as  possible,  but  It  must  be 
recalled  that  a tradeoff  condition  exists  here.  The  pointing  error, 
mirror  surface  error,  tracking  error,  economical  mirror  width,  and 
similar  constraints  limit  the  width  of  mirror  strips  to  a value 
which  Is  probably  considerably  greater  than  the  optimum  suggested 
by  radiative  loss  considerations. 

(2)  Since  at  high  temperatures  radiative  losses  predominate, 
the  greatest  efficiency  Increase  Is  accomplished  by  significantly 
lowering  the  I.R.  emlttance  of  the  aperture  (see  Figure  36).  This 
merits  further  study,  as  will  be  discussed  In  Chapter  VI. 

(3)  Considerable  gains  can  be  made  from  an  Increase  In  mirror 
reflectance  from  0.8  to  0.9  (see  Figure  37;  the  SUF,  as  defined.  Is 
Irrelevant  and  Is  not  shown) . Front-surface  mirrors  with  a reflec- 
tance In  this  higher  range  are  available,  but  degradation  of  the 
reflective  surface  could  pose  difficulties. 

(4)  Some  Improvement  also  results  from  Increasing  the  convec- 
tive cover  transmittance  from  0.85  to  0.95  (see  Figure  38).  In  fact. 


Figure  35.  Effect  of  aperture  width  (theoretical). 


i 

' Figure  36.  Effect  of  aperture  emittance  (theoretical) 
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it  may  prove  advantageous  to  dispense  with  this  cover  altogether. 

The  tradeoff  (not  studied)  with  augmented  free  convection  loss  from  j 

the  aperture  in  the  full-scale  device  would  probably  result  in  a 
net  gain,  given  the  proportionately  small  size  of  the  aperture. 

Again,  the  SUF  is  irrelevant. 

(5)  For  a 30  cm  I.D.  receiver  pipe,  10  cm  of  insulation  is 
about  optimal  (see  Figure  39) . The  reason  behind  the  seeming  insen- 
sitivity of  the  device  to  insulation  thickness  lies  in  the  radiant 
loss  from  the  aperture.  This  is  negatively  affected  by  the  wider 
aperture  sides  (from  thicker  insulation),  although  the  increased 
radiation  is  somewhat  offset  by  a poorer  view  factor  to  the  ambient 
environment.  Thus,  any  gain  from  added  insulation  (beyond  10  cm) 

is  masked  and  overwhelmed  by  associated  radiative  losses.  The  same 
mechanism  operates  in  apparently  making  a 5 cm  thick  insulation 
almost  as  effective  as  the  double  (10  cm)  thickness. 

(6)  Reducing  the  pipe  diameter  to  20  cm,  for  the  same  gas 
mass  flow  rate,  has  a negligible  effect  on  efficiency  (see  Figure  40). 

Use  of  the  smaller  pipe  would  entail  a larger  pressure  drop  (fric- 
tional loss)  which  was  not  included  in  this  calculation.  An 
alternate  conclusion  is  that  an  Increase  in  mirror  width  (i.e., 

flux  concentration  ratio)  is  probably  not  an  effective  means  for 
increasing  the  thermal  efficiency. 
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C.  Conclusions  from  Theoretical  Analysis 

The  good  agreement  between  experimental  and  theoretical 
results  has  allowed  optimization  studies  of  the  full-scale  device 
to  be  performed,  via  a second  computer  model.  These  indicate  that 
average  (over  the  length  of  the  receiver)  heat  utilization  efficien- 
cies (i.e.,  heat  into  fluid/heat  into  receiver,  also  termed  here 
the  Solar  Utilization  Factor,  of  SUF)  on  the  order  of  50%  and  a 
receiver  gas  (He)  exit  temperature  of  ca.  SIO^C  are  realistic 
design  goals.  This  analysis  lends  credence  to  the  concept  as  first 
proposed  (32),  and  confirms  that  the  originally  projected  efficien- 
cies were  not  overly  optimistic. 
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CHAPTER  VI 


SUMMARY  AND  SUGGESTIONS  FOR  FURTHER  STUDY 

A high-temperature  solar  collector  design  has  b^an  described 
herein,  and  the  results  of  experimental  and  theoretical  studies 
have  been  discussed.  In  summary,  these  studies  have  shown  that: 

(1)  High  fluid  exit  temperatures,  at  a reasonable  thermal 
efficiency,  are  achievable  with  this  design. 

(2)  Readily  available  materials  can  be  employed  in  a high- 
temperature  collector  with  little  degradation.  Very  similar 
materials  could  be  used  in  a practical  full-scale  device. 

(3)  The  Archimedes  mirror-concentrator  is  easily  constructed 
and  focused,  requiring  little  maintenance  thereafter. 

(4)  A fixed  glass  receiver  pipe  surrounded  by  rotating 
insulation  and  aperture  is  a simple  and  effective  means  of  radiation- 
cavity  construction. 

(5)  The  use  of  circular  tracks,  chain  drive,  and  a stepping 
motor  is  a method  which  could  be  readily  adapted  to  a full-scale 
tracking  system.  For  example,  earthen/reinforced  concrete  tracks 
could  be  made  at  low  cost  and  to  acceptable  tolerances. 
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of  the  films  are  outstanding  (166) . While  various  methods  exist 
(164)  for  applying  tin  oxide  coatings  to  surfaces,  a straightforward 
technique  such  as  proposed  by  Lytle  and  Junge  (167) , with  boro- 
silicate  glass  as  the  substrate  material,  would  be  worth  investigating. 

A small  pilot  plant,  of  a size  which  could  employ  a commercially 
available  gas  turbine,  would  be  a necessary  element  in  the  progression 
to  full-scale  construction.  At  some  step  along  this  study  sequence, 
attention  should  be  given  to  optimum  receiver  pipe  length,  fluid 
mass  flow  rate,  and  receiver  arrangement  (in  the  last,  tilt  or  hori- 
zontal positioning  tradeoffs  need  to  be  considered).  Last,  but  by 
no  means  least,  a detailed  economic  analysis  of  the  proposed  design 
needs  to  be  made  in  order  to  ascertain  its  viability.  As  the  costs 
of  other  power-production  methods  increase  with  time,  a point  could 
well  be  reached  in  the  near  future  at  which  solar  thermal  power  would 
become  economically  competitive. 
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APPENDIX  I 


COMPUTER  PROGRAM  OF  SCALE  MODEL 


I I 


best  AVMUBIE  COPY 
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>T  TT«-DTl  SOL  4 SRC 


COMMON  PPH,PHG,TAM,T1 ,HFX,  SHD,  ALT,  TB,  TEX,  T I2,DAY ,T 16,ASUF, 
lBSU,B,fiCO  ,ADT 

FORMAT  (IH  , 10X,55HSOLAR  COLLECTOR  TEMPERATURE  PROFILE) 

FORMAT  <IH  , 32HTHE  FOLLOWING  VARIABLES  ARE  USED) 

FORMAT  (F7.2) 

FORMAT  (1H+,13HPPH  = ) 

FORMAT  (1H+,  13HPHG  : ) 

FORMAT  <IH+,13HTAM  = ) 

FORMAT  (IH+, 12HTI  = ) 

FORMAT  (IH+,  13HHFX  = ) 

FORMAT  (1H+,13HSHD  : ) 

FORMAT  (IH+,  13HALT  = ) 

FORMAT  (1H+,12HTB  s ) 

FORMAT  <IH+,13HTEX  s ) 

FORMAT  <1H+, 13HT12  = ) 

FORMAT  (1H+,13HDAY  = ) 

FORMAT  ( 1H+, 13HT16  = ) 

WRITE  (4,1) 

WRITE  (4,2) 

WRITE  (4,4) 

READ  (4,3)  PPH 
WRITE  (4,5) 

READ  (4,3)  PHG 
WRITE  (4,G) 

READ  (4,3)  TAM 
WRITE  (4,7) 

READ  (4,3)  T1 
WRITE  (4,8) 

READ  (4,3)  HFX 
WRITE  (4,9) 

READ  (4,3)  SHD 
WRITE  (4,10) 

READ  (4,3)  ALT 
WRITE  (4,11) 

READ  (4,3)  TB 
WRITE  (4,12) 

READ  (4,3)  TEX 
WRITE  (4,13) 

READ  (4,3)  T12 
WRITE  (4,14) 

READ  (4,3)  DAY 
WRITE  (4,15) 

READ  (4,5)  T16 
ASUF  s 0.0 
6SU  s 0.0 
B = 0.0 
QCO  : 0.0 
ADT  : 0.0 
CALL  S0L5 
CALL  S0L6 
GO  TO  16 
END 
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BESI  AVAIUBIE  COPY 


SUBROUTINE  SOL 5 

COMMON  PPH,PHG,TAM,TI  ,H FX ,SHD .ALT ,TB,TEX ,Tl  2,DAY ,DUM( 2) , 

lOUMl (4),BFX,GM,FL0 

CPT  r PPH*(l .062-0. 0008* TAM) /I .042 

FLO  s EXP(3.41+.4S5*ALOG(CPT)-.062>«EIOG(CPT)**2+.0I61*ALO6(CPT) 

1**3) 

FLO  = FLu*3.59*(72.5+PHG)/(273.0+T12) 

GM  s FLO*0.0297 
VEL  = FLO* 0.93 
FXl  = 36.4*(T1-TAM) 

ANG  z ATAN(91  .5/ALT) 

PRF  r HFX*697.0/S1N<ANG) 

IF  (DAY-183. 0)17,  17,18 

17  ETX  z 1400.0-(0.497*DAY) 

GO  TO  19 

18  ETX  = 1309. 0+0. 497*(DAY-183,0) 

19  RT  z PRF/ETX 

HD  z 3.01-6.85*RT+4.19*RT**2 
FX2  = PRF*(  1.0-HD)*COS(0.86-ANG) 

AFX  z (FXl+FX2)/2.0 
BFX  s AFX*3.67E-3 

WRITE  (4,100)FX1,FX2,AFX  ! 

WRITE  (4,105)  FLO, VEL 
WRITE  (4,110) 

100  FORMAT  (IH  ,SHFX1  = ,F7.1 ,4X,  6HFX2  z ,F7. 1 , 4X  .SHAFX  z ,F7.1)  ; 

105  FORMAT  (1H0,6HFLO  z ,F6.1 ,4X, 10HAND  VEL . z ,F6.1)  L 

110  FORMAT  (1H0,5HPOSN.,4X,2HRE,7X,2HRA,4X,7HT  INCR.,4X, 

16HBULX  T,SX,3HSUF,2X,8HAVG.  SUF,4X,3HTW2) 

RETURN  * 

END  I 


I 

« 


I 

i 


i 


207 


BtSl  AVAIlABLt  COPY 

SUBROUTINE  SOLS 

COnnON  DUn(9),SHD^LTJB,TEX,DUni  (2),TI6,ASUF, 
lBSU,B,QC0,A0T,BFX,6n,FL0 
DT2  s <T16-TB)/A2.0 
QN2  : DT2*Gn*0.921 
DO  40  I : U9t 

VII  s 2.13E-4*(1.0+0.0025*TB) 

RHO  = 0,487/(TB+273.0) 

VI2  : VI I /RHO 

CON  : 1.63E-4*(l.0^0.002S*TB) 

RE  = 4.79E-3*ao/VII 
RA  s 2,75E+3*DT2/(<TB+273.0)*V12**2) 

Z = 1.92=«?E/<32.0+FLOAT(I)-SHD) 

ANU  = 4.18+3.  I7E-20Z-3.67E-9*Z**2+I  .91E-8*Z**3 
IF  (RA>1000.0)24, 21,21 
21  BNU  : 0.31*RA«*0.27 

GO  TO  31 

24  IF  (<RE*RA)-2.75E+9)30,27,27 

27  BNU  = 2.0 

60  TO  31 

30  BNU  : 0.27*(RE*RA)+*0.16 

31  SNU  : 9.74E+34>C0N*ANU»BNU 

TW2  : TB+0.114*ON2/(CON+SNU) 

TAV  = {TV2-TEX)/2.0+273.0+TEX 
QIS2  = 1.04E-3*(TV2-TEX)*<2.87-4.l5E-3+TAV+l.76E-5*TAV**2) 

Q6L2  : 0.90E>lt*(TV2+273.0)**4 

QTR2  s l.79E-ll*TAV**4 

32  QN2  : BFX -GIS2-QGL2-Q00-QTR2 

DT2  : QN2/<Gn*0.52I) 

TV3  s TW240T2 

TAV  s (TW3-TEX)/2. 0+273. 0+TEX 

QIS3  t 1.04E-3*{TW3-TEX)*<2.87-4.l5E-3+TAV+1.76E-5*TAV**2) 

Q6L3  : 0.90E>lt*(TV3+273.0)««4 
QTR3  = 1.79E-I1*TAV««4 
QN3  5BFX-0IS3-Q6L3-QC0-QTR3 
DT3  : QN3/(Gn«0.921) 

QC2  = DT2*3.04 
QC3  : DT34>3.04 
QCO  : (QC2-QC3} 

IF  (ABS(DT2-AOT)-0. 01  >36,36,34 
34  ADT  s DT2 

GO  TO  32 

36  TB  r TB+DT2 

SUF  s QN2/BFX 
BSU  : BSU+SUF 
ASUF  : BSU/FLOATd) 

IF  (I  .LE.IFIXCB))  GO  TO  40 
38  WRITE  (4,119)  I ,RE,RA,DT2,TB,SUF,ASUF,TV2 

B : 6+10.0 
40  CONTINUE 

119  FORMAT  (IH  ,I2,4X  ,F7.  1 ,2X,F7. 1 ,2X,F9.  1 , 9X,F6. 1 , 6X,F4.2,4X ,F4.2, 

I4X,F6.I) 

RHURN 

END 
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Definition  of  Tenns  in  Computer  Program 
(Listed  in  Order  of  Appearance) 


Input  Parameters; 

PPH  measured  AP  of  n-butyl  phthalate,  cm 
PHG  measured  AP  of  mercury,  cm 

TAM  ambient  temperature  (shaded  thermometer  or  thermocouple  T2) , ®C 

T1  measured  temperature  of  Cu  plate  in  flux  bulb  (=  Tl) , "C 

HFX  total  solar  flux  on  a horizontal  surface,  langley  min  ^ 

SHD  measured  length  of  shadow  at  receiver  inlet,  from  beginning 
of  aperture,  cm 

ALT  measured  horizontal  length  of  shadow  cast  by  91.5  cm  high 
vertical  rod,  cm 

TB  bulk  temperature  of  gas  (initially  = T14,  at  m = 0) , °C 

TEX  measured  average  (axial)  temperature  of  insulation  exterior 

(=  (T3  + T8)/2,  = T ) , °C 
e 

T12  measured  temperature  of  gas  at  flow  constriction  (=  T12) , ®C 
DAY  day  of  the  year 

T16  measured  temperature  of  gas  near  axial  center  of  receiver 

(=  T16) , °C,  used  only  at  first  step,  to  start  computation 

Other  Terms  in  Program: 

CPT  AP  of  n-butyl  phthalate  corrected  to  25°C,  cm 
FLO  gas  flow  rate  at  STP  (O^C,  760  torr) , liter  min  ^ 
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gas  mass  flow  rate,  g s 


gas  velocity  at  STP,  cm  s 


calculated  direct  perpendicular  solar  flux,  from  flux  bulb 
AT,  W m-2 


solar  altitude  above  horizon. 


total  perpendicular  solar  flux,  W m“^ 


perpendicular  extraterrestrial  solar  flux,  W m”2 


flux  ratio 


diffuse  fraction  of  total  solar  flux 


calculated  direct  perpendicular  solar  flux,  from  total 
horizontal  reading,  W m~2 


average  of  FXl  and  FX2,  W m”^ 


concentrated  solar  flux  available  at  receiver  Pyrex  aperture, 
W cm~l 


temperature  rise  of  gas  over  step,  (T^)^  - (T^)m-1»  cm" 


net  heat  flux  into  gas  at  step  m,  W cm”^ 


dynamic  viscosity  of  gas  (=  y)  at  g s"-*-  cm 


r'TTl'”! 


density  of  gas  (=  p)  at  g cm"3 


kinematic  viscosity  of  gas  (*  v)  at  cm^  s 


2 0-1 


thermal  conductivity  of  gas  (*  k)  at  (T^,)^.]^,  W cm“l 


Reynolds  number  (*  Re) , dimensionless 


Rayleigh  number  (»  Ra) , dimensionless 
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Z axial  position  from  onset  of  heating  (■  1/z^),  dimensionless 

ANU  ratio  of  Nusselt  numbers,  (developing  velocity  and  temperature 

profiles  Nu)/(fully  established  velocity  and  temperature 
profiles  Nu) , both  Nu’s  at  an  assumed  constant  axial  heat 
input 

BNU  ratio  of  Nusselt  numbers,  (Nu  from  combined  convection  at  41® 
tilt)/(Nu  for  pure  forced  convection,  horizontal  orientation), 
both  Nu's  at  an  assumed  constant  axial  heat  input 

SNU  actual  Nusselt  number,  including  effects  of  developing  flow, 
tilt,  buoyancy,  and  varying  fluid  thermal  conductivity, 
dimensionless 

TW2  average  (peripheral)  temperature  of  heat  transfer  surfaces, 
Pyrex  + graphite,  at  step  m,  ®C 

TAV  mean  (transverse)  insulation  temperature  at  step  m,  ®K 

QIS2  heat  conducted  (lost)  through  insulation  at  step  m,  W cm”^ 

QGL2  heat  radiated  (lost)  by  Pyrex  window  at  step  m,  W cm”^ 

QTR2  heat  radiated  (lost)  by  Transits  sides  of  aperture  at  step  m, 

W cm“l 

TW3  average  (peripheral)  temperature  of  heat  transfer  surfaces, 
Pyrex  + graphite,  at  step  m,  ®C 

QIS3  heat  conducted  (lost)  through  insulation  at  step  m + 1,  W cm~l 

QGL3  heat  radiated  (lost)  by  Pyrex  window  at  step  m + 1,  W cm“^ 

QTR3  heat  radiated  (lost)  by  Transits  sides  of  aperture  at  step 

m + 1,  W cm”^ 

QN3  net  heat  flux  into  gas  at  step  m + 1,  W cm"^ 

DT3  temperature  rise  of  gas  over  step,  (T^j)^j^  - (T^J)^,  ®C  cm”l 

QC2  total  heat  conducted  axially,  Pyrex  + graphite,  at  step  m, 

W cm~^ 
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QC3  total  heat  conducted  axially,  Pyrex  + graphite,  at  step  m + 1, 

W cm”^ 

QCO  net  heat  conducted  axially  (lost),  Pyrex  + graphite,  at  step  m, 
W cm”^ 

ADI  dummy  temperature  used  in  computing  axial  loss 

SUF  Solar  Utilization  Factor,  a measure  of  collector  thermal 

efficiency,  at  step  m 

BSU  axial  sum  of  SUF's 

ASUF  average  SUF  over  axial  length  (■  SUF) 

B dummy  count,  used  to  print  only  every  10th  computation 


<11 


best  available  copy 


PERFORfWCE  ON  29.111.77.  1:15F..^. 


FXl  = 

94S.4 

FX2  s 

774.2 

AFX  = 860 

.3 

n.0  s 

20.9 

AND  VEL. 

. : 19.5 

POSH. 

RE 

RA 

T I NCR. 

BULK  T 

SUF 

A VC.  SUF 

TW2 

1 

323.5 

173.3 

3.5 

180.5  asi) 

0.35 

0.35 

212.1 

(253) 

H 

308.3 

145.7 

2.5 

21T79 

0.25 

0.29 

24^.3 

21 

297.2 

94.2 

2.0 

229.9 

0.20 

0.26 

258.8 

31 

289.0 

63.8 

1.6 

247.5 

0.16 

0.23 

271.0 

41 

282.7 

44.8 

1.3 

261  .5  (272) 

0.13 

0.21 

280.7 

5i 

277.9 

32.3 

1.0 

272.6 

0.10 

0.19 

288  .3 

SI 

274.2 

23.8 

0.8 

281  .4 

0.08 

0.18 

294.3 

71 

271.4 

17.7 

0.6 

288.4 

0.07 

0.16 

299.0 

81 

269.1 

13.4 

0.5 

294.0 

0.05 

0.15 

302.7 

91 

267.4 

10.5 

0.4 

298.4  (316) 

0.04 

0.14 

305.7 

FERFOimCE  ON  19.XI.76.  1:50P.M. 

FXl  s 

582.4 

FX2  s 

550.0 

AFX  s 566, 

.2 

FLO  = 

26.6 

AND  \CL. 

s 24.7 

POSN. 

RE 

RA 

T INCR. 

BULK  T 

SUF 

AVG.  SUF 

TW2 

1 

544.8 

970.7 

3.4 

33.4  (38) 

0.66 

0.66 

79.3 

(90) 

1 1 

510.4 

756.4 

2.7 

67.0 

0.53 

0.58 

119.6 

21 

484.0 

478.5 

2.3 

92.0 

0.46 

0.54 

138.2 

31 

463.2 

517.0 

2.0 

113.7 

0.40 

0.50 

154.8 

41 

446.7 

220.1 

1 .7 

132.3  (131) 

0.35 

0.47 

168.6 

51 

433.3 

153.3 

1 .5 

148.4 

0.30 

0.44 

180.2 

SI 

422.5 

117.1 

1.3 

162.3 

0.26 

0.41 

189.9 

71 

413.6 

88.7 

l.l 

174.1 

0.22 

0.39 

198.2 

81 

406.3 

68.4 

0.9 

134.3 

0.19 

0.36 

205.3 

91 

400.2 

53.5 

0.8 

193.0  (193) 

1 19.111. 77.  1:! 

0.16 

0.3  4 

211.2 

(223) 

PERFORfWiCE  or 

BP.M. 

FXl  s 

855.4 

FX2  : 

388.4 

AFX  : 871, 

,9 

FLO  : 

32.3 

AND  VEL. 

: 30.0 

POSN. 

RE 

RA 

T INCR. 

BULK  T 

SUF 

AVG.  SUF 

TW2 

1 

553.6 

353.7 

3.4 

123.4  (123) 

0.53 

0.53 

175.6 

(life! 

1 1 

523.5 

278,6 

2.7 

152.6 

0.43 

0.46 

210.0 

21 

500.3 

183.3 

2.3 

177.7 

0.37 

0.43 

227.7 

31 

482  .0 

126.2 

2.0 

199.2 

0.31 

0.40 

242.6 

41 

467.4 

90.0 

1.7 

217.5  (22i) 

0.27 

0.37 

255.0 

51 

455.7 

66.1 

1 .4 

233.1 

0.23 

0.35 

265.4 

SI 

446.2 

49.6 

1.2 

246.4 

0.19 

0.32 

274.2 

71 

438.4 

38.0 

1.0 

257.6 

0.16 

0.30 

281.6 

81 

432.1 

29.5 

0.9 

267.1 

0.14 

0.28 

287.7 

91 

426.9 

23.2 

0.7 

275.  1 (281) 

0.12 

0.27 

292.9 

(318) 

PERFORWNCE  ON  15.X.76.  1:30P.M. 

FXl  : 

800.8 

FX2  s 

707.7 

AFX  = 754, 

.3 

ao  = 

44.7 

AND  VEL. 

: 41.6 

POSN. 

RE 

RA 

T INCR. 

BULK  T 

SUF 

AVG.  SUF 

TV2 

1 

853.5 

575.0 

2.6 

69.6  (70) 

0.65 

0.65 

131  .2 

(1801 

II 

812.3 

449.5 

2.2 

93.0 

0.56 

0.59 

159.1 

21 

778.1 

316.1 

2.0 

1 14.3 

0.51 

0.56 

174.2 

31 

749.6 

227.4 

1.8 

133.5 

0.46 

0.54 

188.8 

41 

725.8 

168.6 

1.6 

150.8  (155) 

0.41 

0.51 

201  .3 

51 

705.6 

128.2 

1.5 

i^^I! 

0.37 

0.49 

212.3 

SI 

688.5 

99.4 

1.3 

180.2 

0.33 

0.46 

221.8 

71 

673.9 

78.4 

1.2 

192.7 

0.30 

0.44 

230.3 

81 

661.3 

62.7 

1.1 

203.8 

0.26 

0.42 

237.7 

91 

650.5 

50.8 

0.9 

213.7  (202) 

0.23 

0.40 

244.3 

(240) 
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Defin’.tion  of  Terms  in  Computer  Results 

FXl  calculated  direct  perpendicular  solar  flux,  from  flux  bulb 

AT  (discussed  in  Chapter  IV) , W m"^ 

FX2  calculated  direct  perpendicular  solar  flux,  from  Weather 

Bureau  Station  (discussed  in  Chapter  IV) , W m~2 

AFX  (FXl  + FX2)/2 

FLO  Ar  flow  rate  at  STP,  liter  min"^ 

VEL  Ar  velocity  in  receiver  at  STP,  cm  s~l 

POSN.  axial  position  along  test  section,  cm 

RE  Reynolds  number 

RA  Rayleigh  number 

T INCR.  Ar  bulk  temperature  rise  from  previous  step,  "C  cm”^ 

BULK  T Ar  bulk  temperature  (»  T^j) , ®C,  in  parentheses  are  measured 
temperatures 

SUF  Solar  Utilization  Factor,  a measure  of  collector  thermal 

efficiency  (»  qjj/qs),  dimensionless 

ASUF  average  SUF  over  length,  Z SUF/POSN.,  dimensionless 

TW2  average  (peripheral)  temperature  of  heat  transfer  surfaces, 

Pyrex  + graphite,  ®C,  in  parentheses  are  measured  temperatures 
of  the  graphite  absorber 


APPENDIX  III 


COMPUTER  PROGRAM  OF  FULL-SCALE  COLLECTOR 
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best  AVWlABlt  COPV 

SOLll 

COMMON  DIA  ,TK,WW,EW,ET,RM,TR 

FORMAT  (IW  , I0X;55HSOLAR  COLLECTOR  TEMPERATURE  PROFILE) 
FORMAT  (IH  ,32HTHE  FOLLOWING  VARIABLES  ARE  USED) 


FORMAT  (F7.2) 

FORMAT  (1H+,  I5HDIA  = ) 

FORMAT  (1H+,12HTK  : ) 

FORMAT  (1H+,12HWW  : ) 

FORMAT  (IH  + ,12HEW  = ) 

FORMAT  (IH+,I2HET  = ) 

FORMAT  (IH+,12HRM  : ) 

FORMAT  (1H+,12HTR  : ) 


WRITE  (4,  I) 

WRITE  (4,2) 

WRITE  (4,4) 

READ  (4,3)  DIA 
WRITE  (4,5) 

READ  (4,3)  TK 
WRITE  (4,6) 

READ  (4,3)  WW 
WRITE  (4,7) 

READ  (4,3)  ni 
'<ff?ITE  (4,8) 

READ  (4,3)  ET 
WRITE  (4,9) 

READ  (4,3)  RM 
WRITE  (4,13) 

READ  (4,3)  TR 
WRITE  (4,25) 

FORMAT  (IH  ,5HPOSN.,8X,7HT  I NCR  . , 5X  , 6HB1JLK  T,  9X  ,3HSUF , 7X , 
18HAVG  . SUF) 

CALL  SOL  12 
GO  TO  23 
END 
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BEST  AVAIUBIE  COPY 


SUBROUTINE  S0L12 

COWON  DIA  ,TK,WW,EW,ET,RM,TR 

BSU  : 0.0 

K : 90 

TE  : 349.0 

A : ( TK^^a+l .2*TK*WW+WW**2)**0.5 

Y : 1 .2*TK+WW 

FT  r (TK+Y-A)/(2.0*TK) 

FV'  : (A-TK)/WW 
DO  40  1 : 10,700,10 

3G  : WW*0.1*EW*FV*5.67E-8*(TB+273.0)**^ 

TAV  = (TB-25.0)/2. 0+298.0 

QT  : 0.2*TK*ET+FT*5.67E-8*TAV**4 

B : ALOGC 1 .0+1 .9*TK/DIA) 

QC  = 8.63*(TB-25.0)*(2.87-4.15E-3+TAV+1 .76E-5*TAV+*2)/E 

QS  = 3.2E+4*RM*TR 

QN  = QS-QG-QT-QC 

DT  = QN/4.45E+3 

TB  : TB+DT 

SUF  : QN/QS 

ESU  : BSU+SUF 

ASUF  : 10.0+PSU/FLOAT(I) 

IF  (I  .LE.K ) GO  TO  40 
33  WRITE  (4,35)  1 ,DT ,TE  ,SUF ,A  SUF 

35  FORMAT  (IH  ,I  4, 8X  ,F5 . 1 , 8X  , F S.  I , SX  ,F4 .2, 8X  ,F4 .2) 

K : K+100 
CONTINUE 
RETURN 
END 
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DIA 

TK 

WW 

EW 

ET 

RM 

TR 

POSN. 

T INCR. 

K 


TB 

A.  Y 


FT 

FW 

QG 

TAV 


QT 

B 


Definition  of  Terras  in  Computer  Program 

inside  diameter  of  receiver  pipe,  cm 
insulation  thickness,  cm 

window  width  (i.e.,  aperture  width  at  glass  pipe),  cm 

emittance  of  glass  at  aperture 

emittance  of  aperture  sides  (Transits) 

mirror  reflectance 

convective  cover  transmittance 

axial  position  from  inlet,  m 

temperature  rise  over  element  length,  °C 

counter,  for  printing  every  10th  computation  (i.e.,  every 
100  m of  length) 

bulk  gas  temperature  (=  T^,),  °C 

computed  aperture  dimensions  used  in  determining  the  view 
factors,  cm 

view  factor  of  aperture  sides,  dimensionless 
view  factor  of  glass  window,  dimensionless 
radiant  loss  from  window,  W 

average  (linear)  temperature  of  aperture  sides  at  a given 
axial  position,  °K 

radiant  loss  from  aperture  sides,  W 

accounts  for  mean  insulation  thickness,  dimensionless 
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QC  conduction  loss  through  insulation,  W 

QS  reflected  solar  flux  at  glass  window,  W 

QN  net  heat  into  fluid,  W 

DT  temperature  rise  over  element  length,  “C 

SUF  Solar  Utilization  Factor,  discussed  in  text,  dimensionless 

BSU  ESUF  over  length 

ASUF  average  SUF  over  length  («  SUF) 
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